NASA CR-132494 


EARTH RADIATION BUDGET MEASUREMENT FROM 
A SPINNING SATELLITE: CONCEPTUAL 

DESIGN OF DETECTORS 


by 

L. A. Sromovsky 
K. E. Revercomb 
V. E. Suomi 

V. E. Suomi, Principal Investigator 
L. A. Sromovsky, Co-Investigator 


Prepared under Contract No. NAS-1-13204 by 

University of. Wisconsin 
Space Science and Engineering Center 
1225 West Dayton Street 
Madison, Wisconsin 53706 

for 

NASA Langley Research Center 
Hampton, Virginia 


April 1975 



TABLE OF CONTENTS 


PAGE 

I. INTRODUCTION I_1 

II. SCIENTIFIC REQUIREMENTS II-l 

1. Definition of Planetary Radiation Parameters II-l 

2. Definition of Geographical Radiation Parameters II-4 

3. Required Measurements and Accuracy II-8 

4. Sampling Problems and Requirements 11-10 

5. Deficiencies in Previous Energy Budget Measurements 11-12 


6. Basic Requirements for Earth Energy Budget Observing Systems 11-15 


III. OBSERVING SYSTEM CONCEPT III-l 

1. Plane Flux Measurements, III-2 

2. Measurement of Incident Solar Radiative Flux III-3 

3. In-Flight Calibration HI-3 

4. Sampling Characteristics III- 5 

5. The Satellite Concept III-7 

6. Sensors in Operation III-9 

IV. SOLAR CONSTANT RADIOMETER DESIGN IV- 1 

1. Principle of ACR Operation IV-1 

2. Demonstrated Accuracies of Active Cavity Radiometers IV-2 

3. Design for Present Application IV-3 

4. Servo System Description IV-8 

5. Error Analysis Summary IV-9 

V. PLANE FLUX WIDE FOV SENSOR DESIGN V-l 

1. Basic Construction Geometry V-l 

2. Window Heating Effects V-6 

3. Spectral Response V-8 

4. Angular Response V-12 

5. Required Time Constant V-13 

VI. FAST ACTIVE CAVITY THEORY VI-1 

1. Basic Operating Principles of the Fast ACR VI-1 

2. Cavity Emissivity Calculation VI-5 

3. Analysis of Transient Response VI-7 

4. Linear Servo Analysis VI-10 


5. Angular Response Variations Due to Cavity Non-Isothermalities VI-1* 
VII. FAST ACTIVE CAVITY DESIGN AND ERROR ANALYSIS VII-1 


1. The Design Approach 

2. Design Constraints 

3. The Design 

4. Error Analysis 


VII-1 
VI 1-2 
VII-5 
VII-11 



PAGE 


VIII. 

THERMOPILE DETECTORS AS AN ALTERNATIVE 

VIII-1 


1 . 

Thermopile Characteristics 

VIII-1 


2. 

In-Flight Calibration of the SWS 

VIII-6 


3. 

In-Flight Calibration of the TRS 

VIII-8 


4. 

Deriving Long Wave and Short Wave Fluxes from 




Sensor Outputs 

VIII-11 


5. 

Radiometric Noise Estimates for Thermopile Sensors 

VIII-11 

IX. 

PRE- 

-FLIGHT AND IN-FLIGHT CALIBRATION PLANS 

IX- 1 


1 . 

Calibration Objectives 

IX-1 


2. 

Pre-Flight Calibration 

IX- 3 


3. 

In-Flight Calibration 

IX-5 


4. 

In-Flight Diagnostic Measurements 

IX- 7 

X. 

SYSTEM ERROR SUMMARY 

X-l 

XI. 

INTERFACE REQUIREMENTS 

XI-1 


1 . 

Thermal Interface Considerations 

XI-1 


2. 

Power 

XI-2 


3. 

Weight 

XI-4 


4. 

Signal Interface Requirements 

XI-6 


5. 

Attitude Requirements 

XI-6 


6. 

S/C Spin Requirements 

XI-8 


7. 

Orbital Requirements 

XI-8 


REFERENCES 

APPENDIX A. CALCULATION OF NORMAL AND HEMISPHERICAL CAVITY EMITTANCES 
APPENDIX B. ANALYSIS OF SERVO ELECTRONICS FOR SOLAR CONSTANT RADIOMETER 


APPENDIX C. 
APPENDIX D. 
APPENDIX E. 
APPENDIX F. 
APPENDIX G. 

APPENDIX H. 
APPENDIX I. 

APPENDIX J. 


ERROR ANALYSIS OF SOLAR CONSTANT RADIOMETER 
EFFECT OF WINDOW ON SW SENSOR ANGULAR RESPONSE 
RADIATIVE EFFECTS OF SW SENSOR WINDOWS 

FLUX ERROR PRODUCED BY EMISSIVITY VARIATIONS WITH WAVELENGTH 

FLUX ERRORS RESULTING FROM INCOMPLETE SPECTRAL COVERAGE AND 
IMPERFECT SEPARATION OF SOLAR AND LONG WAVE COMPONENTS 

FLUX ERRORS RESULTING FROM NON-LAMBERTIAN ANGULAR RESPONSE 

ESTIMATION OF FACR ANGULAR RESPONSE ERRORS DUE TO NON-UNIFORM 
RADIATIVE COUPLING BETWEEN CAVITY AND SINK 

PARAMETERIZATION OF THE FACR ANGULAR RESPONSE ERRORS RESULTING 
FROM LEAD WIRE AND SUPPORT CONDUCTION 


APPENDIX K 


THERMAL MODEL OF THERMOPILE DETECTORS 


APPENDIX L. UNCERTAINTIES IN SOLAR INPUT RESULTING FROM ATTITUDE ERRORS 
APPENDIX M. DESCRIPTION OF THE FACR TEST MODEL AND RESULTS 


APPENDIX N 


AC SENSE - DC HEAT CAVITY SERVO DESIGN FOR THE FACR 



ACKNOWLEDGEMENTS 


A large number of people made important contributions to this 

study. Those from the Space Science and Engineering Center were: 

Tom Haig, Bob Dombroski and Bob Wollersheim who provided many 
contributions, especially in the area of spacecraft design, 
interface, and power and weight esitmations; 

Bob Sutton and Bob Herbsleb who fabricated the fast active 
cavity radiometer testi model; 

Mahendra (Mike) Shah who designed the electronic servo for 
the solar constant radiometer; 

Dennis Phillips who performed a preliminary orbit sampling 
study; 

Tony Wendricks and Nai-Kee Mak who drafted many of the figures 
for this report; and 

Nancy Kitzman who did the typing for most of the report. 

The Instrumentation Systems Center, University of Wisconsin 

provided valuable support for the fast active cavity radiometer test 

phase and also performed a study to optimize the electronic servo 

design for the fast active cavity radiometer. We would like to thank 

the following people from the ISC for their contributions: 

Professor Robert Parent and Ken Post who conducted the electronic 
servo study which resulted in the design presented in Appendix N, 
and 

Karl Weitzman who transformed the block diagram of the electronic 
servo for the fast active cavity radiometer test model into an 
operating system. 

We are also indebted to Randy Rasmussen of Electro Optical Systems 
who conducted a preliminary study of the spacecraft power system and pro- 
vided an estimate of the power available from the solar arrays. 

Finally, we enjpyed our interaction with Lou Smith, Langley Research 
Center technical monitor, who made many useful suggestions during the course 
of this study. 



I . INTRODUCTION 


This final report to the NASA Langley Research Center describes the 
results of work performed under NASA Contract NAS-1-13204 initiated on 
6 June 1974 and terminating on 20 April 1975. The objective of this work 
as stated in the Statement of Work is 

"to determine the conceptual design, sensor characteristics, 
sensor performance and accuracy, and spacecraft and orbital 
requirements for a Spinning Wide-Field-of-View Earth Energy 
Budget Detector." 

The guidelines implied by this objective are the scientific requirements 
for measurement of the radiative energy budget of the earth. These 
requirements (stated in Section II) appear to us to demand an observing 
system of at least 6 to 8 satellites in order to obtain adequate space 
time sampling. The economics of multiple satellite systems indicated 
that design simplicity and reliability should also be guiding factors 
in this work. In fact, the idea of using a spinning satellite platform 
for these measurements arose at Wisconsin for these very reasons. This 
system concept, which is described in Section II, was well developed 
prior to the work conducted under this contract and is restated here 
to provide the context within which the present studies were conducted. 
The results of the present study support the original concept, indicating 
that spinning wide f ield-of-view sensors can meet scientific requirements 
and, in fact, have many advantages over earth oriented sensors. 



II-l 


II. SCIENTIFIC REQUIREMENTS 

The purpose of this section is to define the radiation parameters 
which are to be measured, to state the corresponding accuracy require- 
ments for measurement, and to examine the known variabilities of these 
parameters to establish approximate sampling and instrumental require- 
ments . 

1. Definition of Planetary Radiation Parameters 

Figure II. 1 depicts a spherical shell centered over the earth at satellite 
altitude (the exact altitude is not important for the immediate discussion) . 

The coordinate system shown is relative to the earth sun line and the earth's 
orbital plane. Since the x axis points in the direction of the earth's 
orbital motion and the y axis towards the sun, the 0 and (J> coordinates do 
not correspond to earth coordinates of latitude and longitude. Useful 
quantities which can be rigorously measured at the surface of the shell are: 
H(t) = Solar flux incident on the earth at time t 
F^Cr^ , © , c}> ,t) = Plane flux of reflected solar radiation at radius r g , 

coordinates 0, <j> , and time t 

Fg(r g ,0,<t>,t) = Plane flux of earth emitted (long wave) radiation at 

radius r , coordinates 0, d>, and time t. 
s 

At a single instant in time the total reflected solar power is precisely 

equal to the integral of the reflected plane flux F^( r s ,0,<j>,t) integrated 

over the entire surface of the shell of the long wave plane flux F^fr^ ,0,<j>,t) . 

Note that conservation of energy implies that these integrated values do not 

depend on the shell radius r . In order to determine the total solar power 

s 

incident on the earth at time t it is necessary to. define the cross sectional 
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Normal to Earth's 
Orbital Plane 



Figure II. 1 Geometry for defining planetary radiation parameters. Reflected 
radiation passing through an element dS of the spherical shell at satellite 
altitude is F (0,4>)dS. The long wave emitted radiation passing through dS 
is F e (0,4>)dS. r 
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area of the earth normal to the earth - sun line. This is done in terms 
of a reference radius r & corresponding to an altitude 30 km above the 
earth's surface. The expressions for the total radiative power components 
at time t are thus 


rrr^ H(t) = Total incident solar power 


l 

/ 


F (r ,0,<f>,t)ds = Total reflected solar power 
K S 


F (r ' , 0, <f>, t)ds = Total earth emitted power 
-tj s 


From these quantities we can define the more common parameters: the 

planetary albedo A p (t); the average flux of incident solar radiation 
<F s (t)>; the average flux of reflected solar radiation <F R (t)>; the 
average flux of emitted radiation <F„(t)>; and the average net flux 

Ei 

of radiation < ^ NET (t)>. 


A (t) = f F R ds/(irr a 2 H(t)) = <F R (t)>/<F g (t)> 


<F (t) > = irr 2 H(t)/(4Trr 2 ) = -^H(t) 
s a a 4 

<F (t)> = ( F ds/(47rr 2 ) 

K J R a 

s 

<F £ (t)> = $ F E ds/(4irr a 2 ) 
s 

<F NET (t)> = <F £ (t)> + <F R (t)> ~ <F s ( t > > 

= <F £ (t)> - | (l-A p (t))H(t) 

Note that the flux values are defined at the reference radius r . 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


a 

The replacement of the integrals in equations (1), (3), (4), and (5) by 
sums over discrete measurement points is justifiable provided that the 
variations of the plane flux in 0 and <p are adequately sampled by the 
measurement system. At any instant in time 0 and <j> variations due to 
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weather systems on earth would lead to sampling requirements far beyond 
the capabilities of any conceivable measurement system. However, the 
time averaged parameters for a period of the order of a month contain 
much less variability mainly due to the rotation of the earth in the 
0, if) coordinate system. The time averaged parameters are defined as 
follows : 

t+ 2 


<F 


,(t)> T = x f T [fa(t)]dt 


<F E<t> > 


< P E (t) > 


■f. 

■/ 


F r T ds/(4irr a 2 ) 


F T ds/(4Trr 2 ) 
b a 


<F NET (t)> T = <F E (t)> T + <F R (t)> T - <F g (t) > T 
A. p (t) T = <F R (t)> T /<F s (t)> T 


where, if x denotes either E or R, 


F T = 

X 


t/ 


. T 
t+ 2 


T 

t_ 2 


F ( r > t) dt . 
X s 


T 

Note that A p (t) is not equal to the time average of A p (t). 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


( 11 ) 


2 . Definition of Geographical Radiation Parameters 

The previous section dealt with radiation parameters for the planet as 
a whole. Perhaps more significant is the way in which radiative interchange 
is distributed over the earth. The basic character of the general circulation 
of the atmosphere is a result of excess solar energy absorption in the tropics. 
The equator to pole gradient of net flux is the driving force for oceanic as 
well as atmospheric circulation. 
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The desired geographical distribution measurements are measurements 

of local radiation parameters at the reference surface radius r . The 

a 

definition of local parameters is aided by Figure II. 2 which depicts a flat 

surface element receiving an incident radiative flux HcosG from the sun. 

The surface element scatters radiation in all directions in general and 

with highly variable intensity. The integral of all this scattered solar 

radiation is the outgoing reflected flux F at the surface. A similar 

integral of the emitted long wave radiation is the outgoing emitted flux 

F at the surface. The albedo of the surface element is then defined by 
e 

the relation 
F 


A = 


HcosG 


( 12 ) 


and the net radiative loss to space F. T __ is given by 


F..™ = F + F - HcosG = F - (l-A)HcosG. (13) 

iNc-i e r e 

These definitions are applicable regardless of the nonuniformity of the 
surface characteristics and non-Lambertian nature of the scattering. 

However, since both surface character and sun angle can vary with time, 
accurate measurement of these quantities requires instantaneous measure- 
ment of emitted and reflected radiation in all directions. This would 
require so many satellites that they would darken the sky. 

If the surface is uniform (or if it is uniform in a time average 
sense) then it is possible to make an integrated measurement of all 
emitted (or reflected) radiation at all directions by measuring the plane 
flux at (or very close to) the surface in question. This measurement is 
a sampling of the surface and is equal to the outgoing flux averaged over 
the area of the surface. Even if the surface properties are highly 
anisotropic with respect to scattering direction, this measurement yields 



p 


Figure II. 2. Geometry for defining local albedo and local outgoing flux (a), and 
their relationship to plane flux measured at h. The local outgoing flux at 
P is the integral of the intensity at P over all angles (a). The plane flux 
measured at h integrates over intensities and simultaneously averages over 
space (b) . If the region surrounding P has the same scattering (or emitting) 
properties as P then the results are equivalent. 


II-7 


rigorously correct values for A and provided that the assumption of 

uniformity is valid. Measurements of non-planar fluxes do not have this 
property. 

It must be noted that flux measurements at satellite altitude, although 
they provide for rigorous determination of global albedo and energy balance, 
are equivalent neither to the reference level plane fluxes at the subsatellite 
point nor to the average reference level plane fluxes over the sensor field 
of view unless the satellite orbits at the reference level. At altitudes 
required for 5 tq.10 year lifetimes the curvature of the earth can become 
significant. As a result, the sensor field of view, which weights incoming 
intensity as the cosine of the incident angle, does not weight the intensity 
as the cosine of the emitted (or reflected) angle. Although inference of 
local reference level parameters will be more accurate for plane flux 
sensors than would be the case for spherical sensors, errors will be present. 
Errors in this inference are likely to be consistent on the average and affect 
the relative change in reference level radiative parameter estimates far less 
than the absolute values. 

In view of the inability to derive rigorous values of reference level 
fluxed it is better to express the geographical variation of radiative fluxes 
in terms of direct measurements at satellite altitude, i.e., in terms of 
average fluxes at satellite altitude as a function of earth coordinates of 
latitude and longitude. The latter description provides an essential 
boundary condition for climatic modelling and atmospheric and oceanic energy 
transport estimates. This description could also be used to derive approximate 
reference level parameters assuming, for example, that the reference level 


surface is Lambertian. 
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3. Required Measurements and Accuracy 

At present weather and climate models are insufficiently developed to 
provide definite predictive information regarding the effects of changes 
in the earth's radiative boundary conditions. It is in fact a major 
object of earth energy budget measurement programs to provide the long 
term observations to aid in the development and testing of such models. 

In the absence of definitive requirements we must rely on estimates. 

In recent years a number of study groups have dealt with the problem 
of establishing meaningful requirements for earth energy budget measure- 
ments. Reports from these groups (SCEP, 1970; SMIC, 1971; RMOP, 1971) 
indicate an agreement that radiation budget components must be measured 
to an accuracy of at least 1% and perhaps much better. For example, in 
the MIT "Report of the Study of Man's Impact on Climate" (SMIC, 1971), it 
was stated: 

"We recommend 

1. Monitoring the temporal and geographical distribution of the 
Earth-atmosphere albedo and outgoing flux over the entire 
globe, with an accuracy of at least 1%... 

4. Determining the absolute value of the solar constant to better 
than + 0. 5% . . . 

More detailed requirements were outlined very recently in the "Report of 
the GARP Study Conference on the Physical Basis of Climate and Climate 
Modeling" . (1975) . The requirements listed in Table II. 1 make use of this 
later document, as well as previous reports. An attempt was made to 
improve the self consistency of the parameter requirements and to make 
the requirements more specific by stating them in terms of flux errors 
on each component rather than in terms of percentages. It should be 
noted that requirements dealing with local radiation budget measurements 
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have not been included because the spatial resolution required (100 km to 
200 km) is not attainable with plane flux sensors at orbital altitudes. 

4. Sampling Problems and Requirements 

Desired radiation budget parameters are spatial and time averages of 
direct measurements of variable radiation fluxes. Even a perfect radio- 
meter will not yield accurate results unless the variabilities are 

T T 

adequately sampled. In order to estimate the variabilities of F R and F^ 
we shall make use of previous earth energy budget measurements (Vonder 
Haar, 1968). 

Figure II. 3 is a series of meridional plots of outgoing short wave and 
long wave radiation from the earth for three month time averages of non- 
optimally sampled data. Although earth latitude does not accurately 
correspond to 0 in Figure II. 1, the variations shown in Figure II. 3 are 
expected to be comparable in percentages to the variations of F^ with 0. 

The smooth variations shown in Figure II. 3 indicate that polar or near polar 
satellites should provide 0 resolution far in excess of the minimum require- 
ment for seasonal averages; it is expected that monthly and perhaps even 
weekly averages will be adequately sampled as well. 

The variations of F,^ with <j) are expected to be comparable to the 
diurnal variations derived by Vonder Haar (shown in Figure II. 4). In this 
case, although the variation is quite smooth, the difficulty in obtaining 
adequate sampling is much greater than it is for 0 variations if we do 
not want to mix <{> variations with monthly or seasonal variations. Figure 
II. 5 displays estimated outgoing flux errors as a function of diurnal 
sampling intervals based on the data presented in Figure II. 4. 
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MERIDIONAL VARIATIONS 


(a) (b) (c) 





Figure II. 3. Mean meridional profiles of long wave fluxes (a), 
reflected solar fluxes (b) , and net fluxes (c) for each 
season. The solid line in each case indicates the annual 
mean (Vonder Haar, 1968). 


DIURNAL VARIATIONS 



LOCAL TIME (hi) 


Figure II. 4. Diurnal variation of planetary albedo and outgoing 
long wave flux based on TIROS IV measurements (Vonder Haar, 
1968). 
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Since the diurnal modulation of the long wave flux is only + 10% about 
the mean value, a modest diurnal time interval of 2 to 4 hours appears to be 
adequate if errors are required to stay below 0.2%. This time interval 
range roughly corresponds to <J> intervals of 30° to 60° (see Figure II. 1), 
which could be obtained with 3 to 6 polar satellites. 

The 100% diurnal modulation of the reflected solar flux produces much 
more severe sampling constraints based on the estimates shown in Figure 11.5(a) 
It appears, in this case, that even one hour sampling intervals may be 
inadequate. At this point it should be noted that the estimated errors are 
those that would occur if outgoing fluxes were directly averaged over the 
satellite orbit shell without using any external information about the form 
of the diurnal variation. By using such external auxilliary information, or 
by using satellite orbital procession to provide continuous diurnal sampling, 
the actual errors can be reduced substantially below the estimates in Figure II 
5. If we use N near polar satellites providing 2N diurnal samples per day, 
their precessional rates must be large enough to cover the 2N diurnal time 
gaps before seasonal variations become significant. A precession rate of 
360°/2N per month would provide complete diurnal sampling every month and 
probably insure minimal errors due to seasonal effects. 

5 . Deficiencies in Previous Energy Budget Measurements 

The best measurements to date of the radiative energy budget of the 
earth are those of Vonder Haar and Suomi (Vonder Haar and Suomi, 1969, 
and Vonder Haar, 1969) . Their results (probably accurate to + 2%) are 
based on measurements from several satellites, with varying instrumental 
and orbital characteristics. A summary of the instrumental characteristics 


11-13 


(a) Estimated error in average reflected solar flux 
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(b) Estimated error in average long wave flux 
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Figure II. 5. Estimated errors in averages of reflected solar flux (a) 
and long wave flux (b) as a function of diurnal sampling frequency. 
These estimates are for direct averages of discrete samples at fixed 
local times and do not allow for the sampling improvements obtained 
from orbital precession 
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and sampling characteristics of these satellites can be found in a paper 
by Vonder Haar (1969) . Significant shortcomings of these observations 
can be summarized as follows : 

(a) Low resolution spherical sensors (those on Explorer VII, 

TIROS IV, and TIROS VII) suffered from several problems: 

Orbital precession rates required between 65 and 75 
days to obtain maximum diurnal sampling. 

Spherical sensor characteristics produced 3.5 hour 
diurnal gaps. 

Sensors had a crude wavelength descrimination 
capability, all three radiation currents affecting 
each sensor. 

Modeling was required to derive fluxes. 

(b) Medium resolution sensors (MRIR) suffered from very restricted 
angular sampling, incomplete spectral coverage, and post launch 
deterioration, requiring estimate of correction factors and an 
angular distribution model to utilize the data. 

(c) The only plane flux sensors (DISC LRIR's) were in sun-synchronous 
orbits and provided inadequate sampling. 

(d) Latitude coverage was limited in the worst case to + 50° and in 
the best case to + 85°. On the average, sampling of the polar 
regions was inadequate. 

(e) No measurements of the solar output were made thus preventing 
accurate determination of the net radiative interchange. (It 
should be noted that the uncertainty in the solar output was 
not a major error source in this case because of the relatively 
large magnitude of other errors) . 



11-15 


6 . Basic Requirements for Earth Energy Budget Observing Systems 

From the theoretical discussion presented earlier, the work of Bartman 
(1967), Bignell (1961), and others it is possible to list general require- 
ments for measurement systems for determining the radiative energy budget 
of the earth. These should include the following features: 

(a) Measurement of plane flux of radiation from the earth. 

(b) Accurate spectral separation of reflected short wave radiation 
and emitted long wave radiation as well as flat spectral response 
within each range. 

(c) Absolute in-flight calibration capability. 

(d) Monitoring of the absolute radiative flux incident from the sun. 
(This flux normalized to the earth's mean orbital radius is 
termed the solar "constant"). 

(e) Multiple satellites to provide complete spatial sampling and 
rapid diurnal sampling. 

(f) The combined effects of sampling and radiometric errors should 
approach levels indicated in Table II. 1. 
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III. OBSERVING SYSTEM CONCEPT 

The design of detectors for measuring the radiative energy budget of 
the earth is determined not only by radiometric requirements but also by 
the observing system in which they are to be used. The purpose of this 
section is to describe the basic features of the observing system for which 
the spinning wide field of view detectors are designed. These features are 
summarized below: 

(1) Measurement of plane fluxes of reflected solar. radiation and long 
wave emitted radiation at satellite altitude using short time 
constant sensors on board a spinning satellite. 

(2) Periodic measurement of incident solar flux using active cavity 
radiometers on board each satellite. 

/ 

(3) In-flight calibration based either directly on absolute radiometers 
or indirectly with intermediate comparison against the sun. 

(4) Complete spatial and rapid diurnal sampling using 6-8 satellites 
in different orbits launched two at a time. 

(5) Use of simple satellites with 10 year design life which contain 
no moving parts and use spin stabilization and magnetic torquing 
for attitude control. 

(6) Use of a single simple ground station for data readout from all 
satellites and data reduction and archive. 

The observing system concept as outlined above is dictated by the general 
requirements listed in Section II, and by the additional objectives of making 
the system long lived, reliable, and economical both in fabrication and 
operation. These objectives and requirements present a convincing case for 
choosing a simple spin stabilized spacecraft. A spinning spacecraft which 
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permits sensors to view the earth, the sun, and space on each revolution 
provides the capability for accurate in-flight calibration with no moving 
parts. Therefore, the probability for a long lifetime is larger than for 
an earth oriented spacecraft. Also a spinning spacecraft is simpler to 
build than an earth oriented spacecraft making it economically compatible 
with the requirement for a fleet of satellites to properly sample diurnal 
variations. A more detailed system discussion, mainly with respect to 
sensors and sensor operation, is presented in the following subsections. 

1. Plane Flux Measurements 

Measurements of plane flux at satellite altitude will be made in two 
spectral intervals: (1) from 0.2 ym to 4 ym containing approximately 99% 

of the reflected solar radiation; and (2) the complete solar and long wave 
spectrum from less than 0.2 ym to greater than 100 ym. The latter range is 
obtained from a black sensor referred to as the TRS (Total Radiation Sensor) 
and the former using a black sensor, termed the SWS (Short Wave Sensor) , with 
a window probably of fused, water-free, quartz (Suprasil W) . The reflected 
solar flux is determined directly from the TRS and the long wave emitted 
flux from the difference between TRS and SWS outputs. 

Since these sensors are on the rim of a satellite spinning at a rate 
near 6 RPM, they must have moderately short time constants. In order to 
insure response within 0.1% of equilibrium values at sample positions the 
time constant must be -50 milliseconds. The fact that the sensors are 
spinning provides for views of space and incident solar flux which, together 
with absolute solar flux measurements, provide information required for 


in-flight calibration. 
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I 


The sensor field of view, although planar in angular response, will be 
restricted to cover an angular region slightly larger than the angular size 
of the earth at satellite altitude in order to minimize twilight gaps in 
diurnal coverage which result when outgoing fluxes and incident solar 
fluxes both fall within the detector field of view. An illustration of the 
sensor viewing geometry can be found in Figure III.l. The angular size of 
the earth depends on altitude, taking the specific value of 128° for a 
380 nm orbit altitude. The spatial weighting of the plane flux sensors are 
indicated in Figure III. 2 for 700 km and 400 km orbits. 

2. Measurement of Incident Solar Radiative Flux 

An ACR (Active Cavity Radiometer) is used to measure the absolute 
value of the incident solar flux with an estimated error of approximately 
+ 0.1% absolute. This instrument is a primary standard and does not require 
calibration by comparison with any other total radiation instrument. The 
ACR is mounted in the hub of the spinning satellite and is used only when 
the normal to the orbital plane reaches a minimum angular distance from the 
earth-sun line. For a two degree/day orbital precession rate relative to 
the sun the solar flux measurement opportunities occur four times per year 
per satellite. For a six satellite system, this would provide solar flux 
measurements 24 times per year. During each measurement the satellite must 
undergo an attitude change to point the spin axis within 0.5° of the sun. 

3. In-Flight Calibration 

Earth viewing sensors are calibrated by two measurements made of known 
incident flux levels: the zero level is provided by space, and the incident 
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Figure Ill.li Viewing geometry for earth-viewing plane flux sensors. 

Only Earth and Sun sample directions are indicated. 


SATELLITE ALTITUDE = 400 KM SATELLITE ALTITUDE =700 KM 




Figure. II? *. 2. Spatial weighting of plane flux sensors (TRS and SWS) for two 
possible satellite altitudes. Percentage levels describe the 
fraction of received energy coming from within a circle of the 
indicated diameter. 
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solar flux level provides the second point. The value of the solar flux 
level is determined from the solar constant ACR measurement. 

The ACR calibration is based on electrical power measurement standards, 
aperture area measurement, and surface coating emissivity measurement and 
predicted degradation. Since the ACR has a cavity receiver the uncertainty 
in surface emissivity has a greatly reduced effect on the uncertainty of 
the ACR calibration (a reduction by more than a factor of fifty is achieved 
with the proposed design) . 

4. Sampling Characteristics 

The spatial and temporal sampling capabilities of the observing system 
are largely determined by the number of satellites and the configuration of 
their orbits. A fleet of six to eight satellites (shown in Figure III. 3 in 
equally spaced (in longitude) , near polar orbits would provide excellent 
spatial sampling, and with 2°/day orbital precession would provide for rapid 
diurnal sampling (complete diurnal coverage in one month). However, practical 
limitations imposed by booster guidance dispersion and the cost and weight 
required for an orbit plane adjust capability, lead us to believe that a set 
of six to eight satellites in orbit having different inclinations from 20° to 
polar is more practical than a plan based on precise orbit to orbit spacing. 
Preliminary studies indicate that satisfactory coverage may be achieved by 
eight satellites having four different orbit inclinations. These studies 
lead us to believe that system accuracy will not suffer unacceptably with six 
satellites in orbits of different inclinations if the orbits are optimized. 
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Earth Rotates beneath 
Satellites (once per day) 


Terminator 



Satellite 
Rotates 180' 


Satellite orbits 
precess slowly 


Region of Incident 
Solar flux measurement 


North Pole 


Spin Axis 


Spacecraft 



TO SUN 


Figure III. 3. Earth-spacecraft configuration for a six-satellite 
implementation of the radiative energy budget 
observing system. 
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5. The Satellite Concept 

A baseline satellite configuration is shown in Figure III. 4. Many small 
satellites have ’ een designed for, and operated under orbital and other con- 
straints similar to those required for the albedo measuring system. As a 
result, existing designs and proven technology can be taken advantage of. In 
fact, there are no parts of the satellite system which require major develop- 
ment; only adaptions of existing designs to fit together and utilization of 
modern components would be required. 

Because the satellite must operate in all possible positions relative to 
the sun and earth, a shape approaching a sphere is optimum. This same shape 
has, for the same reasons, been used on many satellite programs. The size of 
the satellite is determined principally by the surface area required for 
solar cells. 

The most efficient, most reliable, and longest lived attitude control 
system which will meet program requirements is the spin stabilized magnetic 
torque and magnetic spin control system thoroughly proven in the TIROS and 
military satellite systems. 

The satellites proposed for the albedo and radiative energy balance 
measurements consist of the following basic subsystems: 

(1) Radiometer Subsystem - This consists of earth viewing plane 
flux sensors (TRS and SWS) , the incident solar flux radiometer 
(solar constant ACR) and associated electronics. 

(2) Attitude and Spin Control Subsystem - This consists of magnetic 
coils to adjust spin rate and spin axis by magnetic torquing 
against the earth's field, horizon sensors to provide spin rate 
and attitude error signals, and control logic for adjusting 
torques to match error inputs. 
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(3) Data Collection Subsystem - This provides sample timing, analog to 
digital conversion, data averaging and storage. 

(4) Telemetry and Command Subsystem - This provides for downlink 
transmission of radiometer data and housekeeping information, 
and reception, storage, and processing of commands from the 
ground station. 

(5) Power Subsystem - This consists of solar cells as primary power 
source and batteries to provide continuous operation day and 
night . 

(6) Mechanical Subsystem - Provides support and thermal control for the 
other subsystems. 

A block diagram of the system is shown in Figure III. 5. 

6. Sensors in Operation 

The plane flux sensor operation to measure albedo and earth emission 
is shown in Figure III. 6. The trace labeled "a" shows what is seen by the 
total radiation sensor and the short wave sensor when not seeing the sun. 
When the horizon sensors indicate that these sensors are pointing to the 
nadir, a sample is taken. Shortly thereafter, the sensors are looking at 
empty space and another sample is taken. 

On the sunlit side of the earth, the sensors will also see the sun at 
times when looking away from the earth. The sun pulse is detected and 
sampled at its peak with a photo cell slit trigger system to insure 
proper sampling time. Depending on the position of the sun with respect 
to the earth and the spacecraft, three qualitatively different types 
of signal traces are possible, as shown in Figure III. 6 b-^d. Trace "b" 
illustrates the situation where both an uncontaminated earth and sun 
sample are obtained. If the sun is closer to the horizon than shown in 
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(a) 


(b) 


(c) 


(d) 


Sensor Signal 



(contaminated) (contaminated) 


Figure III. 6. Sensor Operation Diagram for Plane Flux Sensors. In (a) the 

sensor does not see the sun at all during a complete rotation. Traces 
(b) , (c) and (d) show the sun in three different positions to illustrate 
qualitatively different sampling possibilities. A sensor field of view 
restriction of 150° was assumed. Where the signal has contributions 
from both the sun and the earth the individual contributions -rare shown 
as dashed lines. 
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"b", it is possible for the solar sample to be contaminated by earth 
radiation as shown in "c". Such solar data is discarded on the ground 
without loss of the earth sample. Finally, when the sun is very close 
to the horizon it is possible for both the sun and the earth samples 
to be contaminated (see trace "d") . The sampling void caused by such 
cases depends on the sensor field of view restriction. 

It is desirable, and probably necessary, to separate the earth measure- 
ment positively from the sun measurement. To do so, the field of view of 
the flux radiometers will be restricted to about 150° so that the sensor 
sees only slightly more than the earth when the sample is taken. With 
a 150° FOV it will not be possible to obtain isolated "pure" measurements 
of the sun and of space on every satellite spin. However, it is not 
necessary to have such measurements on every spin; it is only necessary 
to collect them frequently enough to prevent data contamination by sensor 
drift or change in solar output, both of which are slow processes. 

With two flux sensors, up to six samples (earth, space, sun) will be 
taken each satellite spin and will be accumulated in six separate registers. 
Every six spins these registers will be averaged and stored in the space- 
craft solid-state memory. This system will provide samples about 216 n. mi. 
apart along the orbit track. Between track spacing with six satellites would 
average about 1800 n miles at the equator and be much closer at higher latitudes 

If data words are 10 bits long, only about 40,000 bits need be stored if 
the satellite is read out twice daily. A single unmanned station located at 
about 45° latitude is sufficient to handle the system even if eight satellites 
are included in randomly drifting orbits since readout can be accomplished 
easily during a small portion of the satellite's pass over the station. 

As the satellite's orbit plane precesses relative to the sun, the sun 
will be seen within a short period of time at all possible angles within 
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the FOV of the plane flux sensors. Thus, even if only one sample is taken 
every six satellite spins, the complete angular response of the sensor will 
be mapped by the sun each six months with very fine resolution. The shape 
of the plane flvx sensor sun pulse will yield information on aging 
characteristics. It should be noted that since the incident solar flux 
radiometer views at 90° to the flux sensors, the observation of the sun 
from one satellite will be used to calibrate the flux sensors in all the 
other satellites — not the flux sensors on board the same satellite. This 
characteristic provides a significant system advantage. With six satellites, 
each reaching calibration position four times each year, each satellite's 
plane flux sensors will be calibrated 20 times each year by five different 


solar reference measurements. 
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IV. SOLAR CONSTANT RADIOMETER DESIGN 

The function of the solar constant radiometer is measurement of the 
incident solar flux with an accuracy approaching .1%. The measurements 
are used to calibrate the solar views of the TR and SW sensors, thereby 
providing accurate measurements of the solar incident flux at daily 
intervals. A series of absolute radiometers known as Active Cavity 
Radiometers (ACR's) has been developed at JPL in recent years. A minor 
modification of one of these designs is especially well suited to this 
task. 

1. Principle of ACR Operation 

The active cavity radiometer consists of a cavity which can be heated 
electrically, a heat sink which surrounds the cavity exterior except for the 
aperture region, and an electronic servo system which controls the cavity 
temperature. The degree of thermal coupling between the cavity and the heat 
sink varies with the specific design. However, all of the designs operate 
using electrical substitution. That is, the power absorbed by the cavity 
from an unknown radiation source is determined by measuring the additional 
electrical power necessary to maintain the cavity temperature when the 
radiation source is removed. The thermal environment of the cavity is 
controlled by the heat sink so that any power losses from the cavity 
(e.g. conduction through the electrical leads, radiative exchange with 
the heat sink, etc.) which occur when the source is present are nearly 
identical to those which occur when it is removed. Therefore, the difference 
of the electrical heating required to maintain the cavity temperature with and 
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without the source present is essentially independent of cavity power losses; 
measurement of the incident source flux depends only on the electrical power 
measurements, the cavity aperture area and the apparent cavity absorptance. 

Since absorptances are notoriously difficult to measure precisely, the ACR 
operation depends heavily on the cavity design which enhances the apparent 
absorptance. The cavity dimensions are chosen such that the apparent 
absorptance is so close to unity that errors in the apparent absorptance 
calculated from the measured surface coating absorptance are much smaller 
than the errors of the measurement . 

2. Demonstrated Accuracies of Active Cavity Radiometers 

A large number of ACR's have been built and operated by JPL since about 
1967. These radiometers fall into two major classes: those in which the cavity 

and the heat sink are thermally isolated (generally referred to as standard 
active cavity radiometers or SACRAD's), and those in which the cavity and heat 
sink are conductively coupled (the first radiometers of this type are 
referred to as primary active cavity radiometers or PACRAD's and the later 
models simply as ACR's). Since the ACR is at present the best radiation 
detector available, it is not possible to measure its absolute accuracy 
directly. However, an important experiment was performed by J. M. 

Kendall, Sr. (1970) which, in conjunction with comparisons of measurements 
made by different types of ACR's, shows that the actual performance of the 
radiometers is consistent with theoretical predictions of their measure- 


ment uncertainties. 
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Kendall used a SACRAD to measure the Stefan-Boltzmann constant which 
can be determined independently to better than .02% using a well established 
relationship involving fundamental constants. The results of measurements at 
four different cavity temperatures show a mean deviation from the theoretical 
value of 0.3%. The theoretical analysis of the SACRAD yields a measurement 
uncertainty of 0.37%. 

A large number of radiometric intercomparisons have been made during 
solar constant measurements from mountain top and balloon. The results are 
summarized below along with the predicted instrument uncertainties. 


THEORETICAL 

UNCERTAINTY 


SACRAD ( In Vacuum) 

0.40% 

(Willson, 

1969) 

PACRAD 

0.22 % 

(Kendall, 

et. al. 1970) 

ACR 

0.22% 

(Willson, 

1973) 


TYPE OF RADIOMETERS COMPARED 
SACRAD'PACRAD 

PACRAD-PACRAD (8 Radiometers) 


MEASUREMENT DIFFERENCES 
<0.5% 

<0.15% 


ACR-PACRAD 


<0.3% 
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The SACRAD-PACRAD comparison was the most difficult to perform accurately 
because the SACRAD had to be operated in vacuum. The results indicate that 
the theoretical calculations of radiometer uncertainties are probably a good 
estimate of the absolute ACR accuracies. 

To further emphasize the capabilities of the ACR it should be noted that 
comparisons have also been made between ACR's and the Angstrom pyrheliometer 
which in conjunction with the Abbott water-flow calorimeter has established 
the International Pyrheliometric Scale (IPS). The theoretical uncertainty of 
the Angstrom pyrheliometer has been calculated to be 2.6% (Willson, 1969). 

The ACR measurements indicate a -2.2% (Willson, 1972) error in the IPS as 
reproduced by Angstrom pyrheliometers. 

3. Design for Present Application 

One of the latest ACR designs, the ACR III, will be used with minor 
modifications for the solar flux measurement. The basic construction of 
the ACR III is shown in Figure IV. 1. The detector cavity is connected to a 
heat sink by a cylindrical thermal impedance. Platinum resistance sensors 
placed at either end of the thermal impedance are used to detect the cavity- 
heat sink temperature difference. A bridge circuit, using the platinum 
resistance sensors as two elements, senses the temperature difference and 
directs power to the cavity heater until the difference is driven to a fixed 
value. Changes in the incident irradiance are thus balanced by changes in 
electrical power. As mentioned earlier, operation of the ACR requires two 
electrical power measurements, one for the unknown source and one for a 
reference source. The cavity thermal environment is maintained constant 
for the two measurements so that the unknown irradiance is proportional to 
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Figure IV. 1. JPL Active Cavity Radiometer (ACR III) Mechanical Configuration 
(from Willson 1971). 
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the difference of the power measurements. 

The ACR III was designed to operate using a shutter for the reference 
source. In this mode of operation the cavity thermal environment is main- 
tained by the heat capacity of the heat sink; that is, the heat sink 
temperature is not actively controlled, but is allowed to drift. To avoid 
the need for a shutter and to control the heating of the heat sink by the 
source Irradiance, a modification will be made to temperature control the 
heat sink to about 10° above the highest spacecraft ambient temperature. 

Platinum heater and sensor wires placed on the heat sink will allow electrical 

_3 

servo control of the heat sink temperature with a reproducibility of +10 °K. 

With this modification, space can be used as the reference source because it 
is not necessary for the solar view and the reference view to be closely 
spaced in time. 

The specific ACR design chosen is illustrated schematically in Figure IV. 2. 
The cavity shape has been altered from that of the ACR III to enhance its 
emittance. A technique for calculating a lower bound to the cavity apparent 
emittance is developed in Appendix A. The results for several surface coating 
emittances are: 

e SURFACE LOWER BOUND TO APPARENT EMITTANCE 
0.88 .9981 

0.90 .9986 

0 .95 .9995 

With a 5% uncertainty or degradation of the surface emittance the apparent 
cavity emittance uncertainty is less than 0.1%. 

The aperture dimensions have also been altered for this application. The 
view limiting aperture has been chosen to block all radiation with an incident 
angle greater than 24°. To restrict the view much further would either permit 




Figure IV. 2. Specific Solar Constant ACR Design Configuration. Dimensions are in 
Centimeters . 


o 
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more solar radiation to be reflected into the cavity causing error, or 
would place more stringent requirements on the pointing accuracy necessary 
to avoid missing part of the sun. For a 380 nm inclined orbit, the 24° 
aperture permits the sun to be viewed without earth radiance contamination 
for at least half of each orbit used for solar flux determinations. The 
pointing accuracy required for this ACR design is about +0.5°. Pointing is 
necessary to provide nearly equivalent electrical and radiative heating of 
the detector cavity and to allow determination of the cavity aperture area 
to +0.05%. 

4. Servo System Description 

The ACR requires very precise control of both the cavity base temperature 
and the temperature difference across the thermal resistance. An electronics 
design study shows that these requirements can be met using state of the art 
technology. (See Appendix B) . The control requirements are as follows: 

Base temperature servo: 

Base temperature setting 300°K + 0.1°K 

Base temperature variation 

between sun and space looks <10 °K 

Base temperature stability 

during 100 second observation <10 °K/100 Sec 

Difference temperature servo: 

Difference temperature (AT) setting 0.3°K + 0.1°K 

Variation in AT between sun and 

space looks <2.5 x 10 °K 

-4 

Stability over 100 seconds <1.0 x 10 °K 
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The difference temperature servo design is shown in Figure IV. 3. The 
important aspects of this design are: • 

(1) The electronics modules which control sensor temperature drift are 
temperature regulated within +0.1°C using an auxiliary servo. 

(2) The base temperatui-e control servo uses the heater winding as a 
temperature sensor. This avoids oscillation problems that would 
arise from heater-sensor time delay. 

(3) The difference temperature control servo uses two platinum wire 
sensors in a bridge circuit. Bridge output (i.e., error voltage) 
is amplified and sign detected. 

Depending on the sign of the bridge output, count in a 16 bit up/down 
counter is updated every time the slow clock pulse occurs. The counter 
output drives a high resolution (16 bit) digital to analog (D/A) converter 
whose output voltage is proportional to the binary value of the 16 bit 
count. D/A output drives a stable constantan wire heater via a unity gain 
power amplifier. Thus, at the end of every slow clock pulse the heater 
power is varied (adjusted) to control the difference temperature without 
an overshoot or oscillation. Moreover, the high resolution 16 bit output 
from the counter provides an accurate high resolution measure of heater 
voltage and hence heater power. 

5 . Error Analysis Summary 

An error analysis of the ACR design shown in Figure IV. 2 is given in 
Appendix C. The basic operating equation of the ACR used there is as 
follows : 

P '-P 

e e 

A c (a c + Pp c } 


H = 




Uj 

. O fy 

Uj ^ p: 
-J ^ Uj 

cq IS 

ssS 

« is 5 


Figure IV. 3. ACR Difference Temperature Control Servo 
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where 

H = Solar incident flux 

P^ = Electrical power input to the cavity while viewing the sun 

P ' = Electrical power input to the cavity while viewing space. 

= Cavity aperture area. 

= Cavity apparent absorptance. 
p = Reflectance of the upper heat sink cavity. 

= Reflectance of the cavity and cavity aperture. 

The analysis shows that a small correction to this equation is necessary to 

account for the cavity heat sink temperature difference increase of 
-4 

2.5 x 10 °K which occurs when the sun is viewed. The correction is 0.13 mW 


or 0.1% for this design. The corrections due to the difference in radiative 
sink-cavity exchange between sun and space views are negligibly small. 

The solar constant measurement uncertainty was calculated by summing 
the standard deviations of the independent cavity parameters (including 
the temperature uncertainties) according to the equation: 



1/2 


where a(£^) is the standard deviation of parameter The largest source 

of error is the cavity absorptance which was assumed known to +0.1%. The 


result for the standard deviation of the solar constant measurement, which 

2 

can probably be slightly improved by optimization of the design, is 0.17 mW/cm 


or 0.13%. 
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V. PLANE FLUX WIDE FOV SENSOR DESIGN 

Two wide field of view (FOV) sensors are required for the earth energy 
budget observations : a short wavelength sensor (SWS) and total radiation 

sensor (TRS) . Their angular and spectral ranges are listed in Table V.l. 

Table V.l. Wide FOV Sensor Identification 


SENSOR 

ANGULAR 

SPECTRAL 

SPECTRAL 

DESIGNATION 

FOV 

FILTER 

RANGE 

SWS 

150° 

Suprasil W 

0.2u- 5p 

TRS 

150° 

None 

<0 . 2p->50u 


The major requirements of these sensors, in addition to spectral response 
are as follows: 

(1) Short time constant (^50 ms) 

(2) Flat spectral response within spectral range 

(3) Lambertian angular response within angular range 

(4) Linearity of the order of 0.1% or better 

(5) In-flight or absolute calibration capability. 

Two approaches are being considered for these sensors. The first uses a short 
time constant active cavity radiometer and the second is based on thermopile 
detectors. The fast ACR appears most attractive at the present. This 
section will discuss design considerations common to both detector 
options. Specific detector characteristics and a comparison will be 
made in the following sections. 

1. Basic Construction Geometry 

Common features of the wide FOV sensors are indicated in Figure V.l. 


These consist of 
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(1) a temperature controlled heat sink, 

(2) a radiation receiver (a defined area within which the detector 
responds to incident radiation) , 

(3) an external aperture defining the angular FOV, and 

(4) a hemispherical window (SWS only) . 

Geometrical parameters identified in Figure V.l are defined below: 
a = radius of an assumed circular radiation receiver 
b = radius of the external aperture 

c = altitude of the external aperture above the plane of the 
radiation receiver 

r = outer radius of the hemispherical window centered at the 
receiver center 
d = thickness of the window. 

The finite detector size has two effects on angular response 
characteristics: (1) the angular FOV is not precisely determined by the 

aperture but has a transition region between zero and full detector 
irradiation; (2) the angular response of the SW sensor is non-Lambertian 
within the region of full detector irradiation because not all parts of 
the radiation receiver are at the center of curvature of the window. 

The first effect can be simply described by the angular size of the 
transition region A0 where A0 is given by 

A0 = tan ■'"[c/Cb-a)] - tan ^[c/Cb+a) ] . (1) 

Note that A0 is the difference between the half angle corresponding to full 
receiver irradiation 0^ and the half angle corresponding to zero irradiation 
© 2 » where 

= x/2 - tan 1 [c/(b-a)] (2) 

© 2 = tt/ 2 - tan 1 [c/(b+a)] 


( 3 ) 
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In order to insure full response to the earth, the half angle 0^ must be 
at least as large as the sum of the half angle subtended by the earth 0 £ 
and the attitude angular tolerance 60. At an altitude of 380 nm we find 
the © E = 64.22°. Allowing for an attitude uncertainty 60 = 0.75°, and an 
angular margin 0 for assuring full view of the earth for 7 detector time 
constants, a reasonable choice for 0^ is 75°. This leads to the condition 
c = (b-a) tan 15°. (4) 

Given a receiver radius a> we can thus determine c and A0 as a function of b. 
Results for a = 0.1 cm (typical of thermopile detectors) and for a = 0.5 cm 
(typical of the fast active cavity detector) are given in Tables V.2 and V.3. 


Table V.2. Aperture Design Parameters for a = 0.1 cm (TP Detectors) 


APERTURE 
RADIUS b 

APERTURE 
ALTITUDE c 

A0 

DIURNAL SAMPLING 
GAP (60 + A0 + 0 ) 

0. 3 cm 

0.054 cm 

7.4° 

1.21 hr 

0. 5 cm 

0.107 cm 

4.9° 

1.04 hr 

1.0 cm 

0.241 cm 

2.6° 

0.89 hr 

1.5 cm 

0.375 cm 

1.8° 

0.84 hr 


Table V.3. 


Aperture Desis 
Detectors) 


Parameters for a = 0.5 cm (Fast Active Cavity 


APERTURE 

APERTURE 

RADIUS b 

ALTITUDE c 

1.0 cm 

0.134 cm 

2.0 cm 

0.402 cm 

3.0 cm 

0.670 cm 

4.0 cm 

0.938 cm 


DIURNAL SAMPLING 


A0 

GAP (60 + A0) 

9.9° 

1.38 

hr 

5.9° 

1.11 

hr 

4.2° 

1.00 

hr 

3.2° 

0.93 

hr 


5.0 cm 


1.206 cm 


2 . 6 ° 


0.89 hr 
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The diurnal sampling gaps tabulated are local time intervals (near sunrise 
and sunset) during which measurements cannot be used because direct solar 
radiation can enter the detector FOV. From the tabulated results an 
aperture radius of 1 cm for TP detectors and 3 cm for cavity detectors 
appears reasonable. 

An important constraint on window size is its effects on angular response. 
The result of detailed calculations presented in Appendix D can be expressed 
as 

R(0) = cos 0 [1-E(0) ] (5) 

where R(0) is the relative power incident on the receiver from an infinitely 

distant point source at angle 0 away from the receiver normal. The cos© factor 
is the ideal response of a plane receiver, and E(0) is the deviation produced 
by the window as a function of angle. R(0) is normalized so that 

R(0) = 1. (6) 

The approximate expression for E(0) is 

E(0) = sin 2 0 a 2 d(n 2 -l)/(2nV) (7) 


where n is the index of refraction of the window. Equation (7) is strictly 
valid only for 

( 8 ) 


4-r « i, i « -St. 

r-d r n-1 


Requiring the maximum value of E(0) to be less than 10 , which insures a 

negligible response deviation, yields the condition 


d < 10 


3 3 

-3 2n r 


2 2 

a (n -1) 


( 9 ) 


where, for fused quartz, n ^ 1.46. For the two a values of interest we 


ft 


have the conditions 
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d < .55 cm r , a = 0. 1 cu (10) 

d £ .022 cm a = 0.5 cm. (11) 

Choosing r = 1 cm for a = 0.1 cm and r = 3 cm for a = 0.5 cm yields 

d < 0.55 cm; a = 0. 1 cm, r = 1 cm (12) 

d £ 0.59 cm; a = 0.5 cm, r = 3 cm (13) 


These conditions leave considerable leeway in choosing window size, 
although the conditions are not exact because the upper bounds , especially 
for (12), do not satisfy (8) very well. 

2. Window Heating Effects 

During its rotation about the S/C spin axis the sensor will be exposed 
to a time varying long wave flux from the earth. The SW sensor window will 
absorb this flux and thus experience radiative heating. Window temperature 
excursion caused by this heating will produce background flux variations 
within the detector FOV. Errors in the SWS flux readings caused by this 
effect can be reduced to negligible levels by invoking appropriate constraints 
on window design. 

The results of the detailed analysis presented in Appendix E can be 
summarized fairly simply. If the time dependent part of the incident long 
wave flux is given by 

F(t) = F q cos cut, (14) 

then the induced background flux variation is given by 

F (t) = 4aT (cucpd) ^F cos (out + (j>) 

15 W O 


( 15 ) 
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where the parameters employed are defined as follows: 
a) = angular rotation frequency of the S/C 
c = specific heat of the window material 
p = density of the window material 
d = window thickness 

F q = half the flux difference between earth and space views 
T w = mean window temperature 

= a phase shift approximately equal to tt / 2 (exactly defined 
in Appendix E) . 

For a fused quartz window and a 6 RPM spin rate we have the following 
specific parameter values: 
ui = 0.628 sec d 

c = 0.753 W-sec/gm-°C 

3 

p = 2.203 gm/cm . 

_2 

Choosing F q = 150 W-M and T g = 300°K we can rewrite (15) in the specific 
form 

F (t) = 0.125 ~ -§S cos (cot + (j>) (16) 

B M 2 d 

where, for a ^ 1 cm, |<f> — tt / 2 | ^ .004 (see Appendix E) . At the Nadir position 
(earth sample) F (t) = 0 since t = 0. The background value for the space 
sample or solar- sample depends on the sampling time. If either is taken at 
anti-nadir position no measurable error is present. For other positions the 
error can be estimated in terms of the angular distance between the sample 
position and the anti-nadir position. If the angular deviation is 0 g then 
the background flux is given by 

= 0.125 sinG . 

B x , 2 J s 


(1.7) 
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The worst case value of 0 is determined by the sensor FOV and has the 

value 0 ~40°. For this case we have 

s 

F = .08 W-M -2 cm/d. (18) 

15 

In order to meet the absolute accuracy target expressed in Table II. 1 
(Section II) for this worst case condition, without correction, the 
window thickness should satisfy 

d >, 0.16 cm. (19) 

A thinner window could be used if window heating effects (which are 
fairly predictable) were corrected for in data analysis, or if samples 
of space or sun were restricted to a smaller angular deviation about the 
anti-nadir point. None of these comments apply to the TR sensor since it 
has no window which can produce such effects. 

3. Spectral Response 

The objective for both the TRS and the SWS is to obtain flat spectral 
response within the spectral passband. In both cases the spectral response 
uniformity is limited by the characteristics of the radiation receiver, which 
in turn is limited more or less by the characteristics of a black painted 
surface. In the case of the fast cavity radiometer this limitation is 
considerably reduced. 

Appendix F contains an analysis of long wave flux errors produced by 
spectral variations in receiver emissivity. These errors can be described 
in terms of the variation of average emissivity with variations in the 
spectral distribution of incident flux. Defining 
e(v) = emissivity at wavenumber v 
F(v) = incident flux at wavenumber v 
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the average effective emissivity for the given incident spectral flux is 


given by 


e(v)F(v) dv/ / F(v) dv 


Average emissivities for spectral fluxes typical of earth emissions are 

compared with the average emissivity for blackbody radiation at 270°K in 

Table V.4. The spectral emissivity values e(v) are for flat painted surfaces 

using Parson's Black in one case, and 3M Black in the other (See Appendix F) . 

The value of emissivity that would be obtained during calibration (at 270°K) 

differs from the average of the e for earth fluxes. This would result in 

r 

a bias error as well as a variable deviation about the mean error. Both 
errors are listed in Table V.5. These values are close to exceeding the 
absolute accuracy target listed in Table II. 1. Allowing for degradation effects 
these may be marginally acceptable. The situation is much better for the 
fast cavity radiometer because of the emissivity enhancement of the cavity 
receiver. For this case errors are estimated to be a factor of ten lower 


than those listed in Table V.5. 


The short wavelength sensor is affected not only by the spectral 
variations in the receiver emissivity but also by the spectral variation 
in window transmission. The spectral variation of transmission for 
Suprasil-W (a water-free fused silica material produced by Amersil, Inc.) 
is displayed in Figure V.2. For most of the solar spectral range the 
transmission is exceptionally flat. The wavelength region which does 
show significant variations is near the long wavelength cutoff (-*-4p) . 
Since window thickness is a major factor in determining SWS cutoff it 
must be chosen to minimize the error involved in spectral separation of 
reflected and long wave radiation components. The results of an analysis 
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Table V.4. A verage Emissivities as a Function of Spectral Flux 
Distribution for Two Black Surfaces 


SPECTRAL FLUX 
CLASSIFICATION 


e FOR e FOR 

PARSON' S BLACK 3M BLACK @373°K 


Blackbody @270°K 
Winter Arctic Atm. 

Western US Summer (Clear) 
Cloudy Tropical Atm. 


0.917 

0.933 

0.918 

0.933 

0.922 

0.926 

0.919 

0.929 


Table V. 5. Flux Errors Produced by Nonuniform Emissivities 


PARSON'S BLACK 3M BLACK @373°K 


Bias error using 270°K 
Calibration 

-0.3% 

+0.4% 

Standard Deviation about 
Mean Error 

0.2% 

0.3% 



TRANSMISSION OF SUPRASIL W 


V-ll 



Figure V.2.. Window transmission as a function of 
wavelength and thickness. 
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of this problem (Appendix G) , which are by no means conclusive, indicate 
an optimum thickness between 1mm and 2mm. This is not in serious conflict 
with other constraints on window thickness. 

We do not have accurate information on SWS flux errors resulting from 
receiver emissivity variations. An analysis of this problem is made 
especially difficult by lack of an adequate means for parameterizing spectral 
variations of reflected flux. At present this is an uncertain area for 
thermopile detectors, and probably very safe for the fast cavity radiometer 
because of the cavity emissivity enhancement (and the resulting factor of 
twelve decrease in sensitivity to emissivity variations). 

4. Angular Response 

For both SW sensor and the TRS sensor angular response is determined by 
the angular response of the radiation receive - '-e ' ave already seen that 
the focusing effects of the SWS window can easily be made negligible) . 

The plane flux errors resulting from non-Lambertian angular response 
are treated in Appendix H for an angular response model of the form 

R(0) = (1-a) cos0 + acos^0, (21) 

where 0 is the angular distance from the receiver normal and ot is an error 
parameter. Equation (21) is equivalent to an angular emissivity variation 
of the form 

e(0) = e [l-a(l-cos0) ] . (22) 

o 

The flux errors resulting from such a variation depend on the angular 
distribution of incident fluxes as well as u. Table V.6 summarizes errors 
for limb darkening and brightening (100% modulation from nadir to limb). 
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Table V.6.. Percentage Flux Errors Resulting as a Function of the 
Angular Nonuniformity Parameter a . 


a 

LIMB DARKENING 

LIMB BRIGHTENING 

.01 

+0.2% 

-0.2% 

.02 

+0.4% 

-0.3% 

.05 

+1.1% 

-0.8% 

.10 

+2.1% 

-1.7% 


Measured data for 3M Black (40.1-C10) are presented in Figure V.3. Although 
they do not fit the model (Equation 22) very well, an a value of 0.04 seems 
to be the best characterization. According to Table V.6 the measured response 
could lead to errors of +1% in the visible flux (where large nadir to limb 
modulations are frequent) and probably much less than this in the IR where 
only (relatively small) modulations are typical. However, if the SW and TRS 
sensors have different angular variations in emissivity, then both reflected 
solar and long wave fluxes could be significantly in error. 

Present information suggests that thermopile detectors may not be able to 
obtain adequate angular response unless surface coatings better than 3M (401-C10) 
can be employed. This is not the case for the fast cavity radiometer which 
again benefits from the cavity receiver. For this detector angular response 
errors can be reduced below 0.2% even using 3M(401-C10) for the interior 
coating of the cavity. 

5. Required Time Constant 

Earth viewing sensors rotate at 6 RPM (nominal) about an axis normal to 
the orbit plane. During each rotation three samples of detector output- are 
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FIGURE V.3. Angular emissivlty variation oi 3M 40I-G1U. 
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recorded (if viewing requirements are satisfied). These samples are listed 
below: 

(1) Earth-only (taken at nadir) 

(2) Sun-only (taken at the rotation angle corresponding to the 
minimum angular distance between the sun and the sensor normal) 

(3) Space-only (sampling angle not critical except for the constraint 
that the edges of the detector FOV be at least 5° away from both 
sun and earth nearest points). 

It can be shown that near nadir and solar sample angles the sensor 
radiation input is approximately proportional to the cosine of the angular 
deviation of the sensor normal to the nominal sample angle. For the purpose 
of estimating the required time constant we shall assume that the sensor 
signal varies as the cosine of the rotation angle for large angular deviations 
as well as small ones. 

If the sensor time constant is x then the sensor output at time t in 
response to a step function radiation input AF at time t' is given by the 
following expression 

F (t) = AF [l-e' (t_t ' )/T ]. (23) 

m 

The response to a time dependent function F(t) is just the integral of the 
responses to the differential step changes in F(t) , i.e 

[l-e- (t - t,)/T ldf (24) 

o 

where we assume F(t') = 0 for t' ^ 0. Since F(t’) has been assumed to have 
the form 

F(t') = F sin (mt') , 
o 


( 25 ) 
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where w is the spacecraft angular rotation frequency, the sensor response 
is found to be 


ir/2u) 


F m (t) = to F q J cos (cot ' ) [1-e (t t,)/T ]dt'. 
o 


(26) 


Integrating (26) and writing the solution in terms of the fractional 
deviation from full response yields 


F -F (T) „ , 

— - Z (1 + (tor) ) 


(27) 


In order to keep the response within 0.1% of F q , t must satisfy the 
condition 


t < 0.03co 


-1 


(28) 


For to = 0.628 sec ^ (corresponding to 6 RPM) we find the requirement 

t <48 ms. (29) 

It should be noted, however, that (29) is an accurate condition only for a 
plane sensor with a hemispheric field of view (in. fact, for this case, it is 
a conservative condition). For a plane sensor with an aperture limited FOV 
the assumed cosine behavior only applies over an angular region equal to the 
difference between the angle subtended by the earth and the angle subtended 
by the aperture. In order for the cosine assumption to be valid starting 
seven time constants prior to the nadir position, the aperture clearance 
should be approximately 12° in half angle. This would lead to a twilight 
diurnal sampling gap of approximately one hour in local time. This should 
not result in significant diurnal sampling errors because the level of 
reflected fluxes during twilight is relatively low. 
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VI. FAST ACTIVE CAVITY THEORY 

The advantages of cavity radiometers have been referred to several 
times in Section IV and in Section V in reference to the solar constant 
radiometer. Among these advantages the following should be emphasized. 

(1) Flat spectral response 

(2) Lambertian angular response 

(3) Absolute calibration 

For the wide FOV earth viewing sensors we must also obtain a relatively 
small time constant in the range from 40 to 60 milliseconds. The purpose 

of this section is to present the theoretical basis for design of a fast active 
cavity radiometer (FACR) without sacrificing any of the major ACR advantages. 
Laboratory measurements verifying the theoretical predictions have been made 
with a FACR test model described in Appendix M. 


1. Basic Operating Principles of the Fast ACR 

The basic design features of the fast ACR cavity and viewing aperture 
are illustrated in Figure VI. 1. (The specific design details are discussed 
in the following section). The nearly spherical cavity is a thin shell formed 
of a thermally sensitive resistance wire coated with a high emissivity paint. 

The resistance wire acts both as the cavity heater and as the cavity temperature 
sensor. The cavity is conductively insulated from its heat sink such that 
the primary coupling between the cavity and the heat sink is radiative. The 
operating principle of this ACR is fundamentally the same as that for the 
ACR discussed for the solar constant measurement. Electrical power supplied 
to the cavity heater is servo controlled to maintain a fixed temperature 
difference between the cavity and the heat sink. Two measurements of the 
electrical power supplied to the cavity are again required, one obtained 




FIGURE VI. 1 Basic design concept for a "high speed" ACR. 

The wire wound cavity provides heating and 
sensing functions. 
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while viewing the unknown source and one obtained while viewing space. 

The difference of these measurements is proportional to the unknown 
irradiance with a proportionality factor depending on the cavity area 
and absorptance. 

Two different operating modes have been considered for the fast 
ACR during its development. For one, the cavity and the heat sink are 
operated at the same elevated temperature somewhere near 400°K. This 
forces the net radiation from the cavity to always be outward such that a 
constant temperature can be maintained by electrical control. For this 
mode of operation the radiation coupling between the cavity and the heat 
sink should be minimized to minimize errors. The second mode relies on 
a large radiative coupling between the cavity and the heat sink to transfer 
the power input to the cavity to the heat sink. Since the radiation resistance 
is necessarily large, a temperature difference of about 20° between the cavity 
and the heat sink is required. The second mode was chosen for two reasons: 

(1) cavity operating temperatures could be reduced and thus the thermal 
requirements on high emissivity coatings could also be reduced; (2) the 
presence of an IR blocking window did hot allow the first method to be 
used for the short wavelength sensor. 

The ACR considered for earth viewing is different from the solar flux 
ACR in two major aspects. First, because the earth viewing sensors view 
the earth, the sun, and space for each revolution of the spacecraft, the 
time constant of the ACR must be about 50 ms while the time constant for 
the solar output ACR is essentially unconstrained. (The time constant of 
the fast ACR is such an important aspect of its operation that it will be dis- 
cussed in considerable detail in the following subsections). Second, while 
the solar constant ACR is always measuring normally incident radiation this 
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is not true for the earth viewing radiometer; thus its design must insure 
that the response is independent of the portion of the cavity heated by 
the incident flux. This cannot be achieved by coupling the cavity to the 
heat sink conductively unless the conducted power is uniformly distributed 
over the cavity surface. However, using radiative coupling to the sink 
and a nearly spherical cavity shape will essentially guarantee that the 
response will be the same for all angles of incidence because all por- 
tions of the cavity have nearly identical views out the aperture and 
towards the heat sink. The importance of the geometry for radiation 
exchange between the cavity and the sink in obtaining uniformity in 
angular response is discussed in a separate subsection. 

The short cavity time constant is obtained using a high gain servo 
system to enhance the natural time constant. The natural time constant 
for the establishment of a steady state following a discontinuous jump 
in the input power to a cavity connected by a thermal resistance R to a 
heat sink is given by t = RC, where C is the cavity heat capacity. With- 
out servo control the cavity temperature changes exponentially until a 
temperature change equal to the product of the thermal resistance and the 
input power change has been attained. For the earth viewing ACR the 
resistance R is a radiative resistance which is given by 

R = ( T "T )/eA0(T 4 -T ' 4 ) = [ eAo (T 2 + T 2 ) (T + T ) ] _1 , (1) 

cs cscs 

where e is the effective emittance for radiative transfer between the 
cavity and the heat sink, a is the Stefan Boltzmann constant, T £ and T^ 
are the cavity and sink temperatures respectively, and A is the exterior 
cavity area. For the cavities under consideration this time constant is 
on the order of 30 sec. The desired time constant of 40 to 60 ms represents 
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an enhancement factor of approximately 600. With sufficient gain the 
servo system can detect very small changes in cavity temperatures and 
rapidly change the electrical power to the cavity. When the net power 
to the cavity reaches equilibrium no further temperature change can occur. 
Although the initial change in cavity temperature as a result of the change 
in irradiance proceeds at a rate determined by the natural time constant, 
the servo can produce power equilibrium long before the full temperature 
excursion can occur. If T is the servo gain in units of electrical power 
input to the cavity per degree change in cavity temperature, the radiometer 


time constant can be expressed as 


RC 

T 1+fR • 


( 2 ) 


Since TR is much greater than 1, we see that the radiative resistance R 
becomes unimportant and the time constant is given by C/T. Detailed 
analysis of the fast active cavity transient response can be found in 
the following sections. 


2. Cavity Emissivity Calculation 

The variations of the apparent emittance of a cavity with wavelength and 
angle can be made much smaller than those for a flat plate. This point was 
mentioned in regards to the solar constant radiometer, but is discussed in 
additional detail here because of its relevance to the angular response of 
the earth viewing sensors. 

The angular response of a cavity of course depends on its shape. We 
will restrict attention to spherical or nearly spherical cavities because 
isotropic emission and absorption is desired. An isothermal spherical 
cavity with a diffusely emitting and reflecting surface coating displays 
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an isotropic angular response because all portions of the cavity have 
identical view factors to the aperture. Although a flat plate has this 
response by definition for a diffusely reflecting and emitting surface 
coating, deviations of the coating from this ideal are suppressed by a 
cavity design. The apparent hemispherical emittance of an isothermal 
spherical cavity with a diffusely emitting and reflecting surface coating 
is given by 

£ a = e + (l-e)(l-A /A ) (3) 

c s 

where A^ is the cavity area and A^ is the area of the whole sphere. The 
derivation of this result is presented in Appendix A. The apparent 
hemispherical emittance for two values of the surface emittance e are 
given in the Table VI. 1 as a function of the ratio of the aperture radius 
to the cavity radius (^) . The third column labeled Ae /At is the ratio 
of the apparent cavity emissivity change to the change in surface emittance. 
The advantages of the cavity radiation receiver are apparent from the Table. 
If the surface emittance varies with angle, wavelength, or as a result of 
degradation, the apparent emissivity of the cavity varies a great deal less. 
For example, for a radii ratio of 0.5 the cavity emissivity changes by only 
0.4% if the surface emissivity changes by 5%. This same radii ratio with 
the 3M Black (401-C10) coating discussed in Section V would insure that 
flux errors due to angular response variations would stay well below +0.2%. 
Flux errors due to spectral variations of emissivity would be even smaller 
(probably below +0.05%). 
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r/R 

e (e = 0.9) 
a 

e (e * 0.95) 
a 

Ae /Ae 
a 

.2 

.9989 

.9995 

0.012 

.3 

.9975 

.9988 

0.026 

.4 

.9954 

.9978 

0.048 

.5 

.9926 

.9965 

0.078 

.6 

.9890 

.9948 

0.116 

.7 

.9844 

.9925 

0.162 



The transient behavior of the fast cavity radiometer depends on 
both the servo characteristics and the thermal properties of the cavity 
receiver. The latter properties require the following parameters for 
description : 

T^Ct) = Cavity temperature at time t 

T 0 = Equilibrium cavity temperature with zero incident flux 

T g = Sink temperature 

= Cavity aperture area 
= Cavity exterior surface area 
e = Effective cavity emissivity within its aperture 

SL 

e s ' = Cavity effective emissivity for radiative exchange with the sink 

= Conductance between cavity and sink via leads, supports, etc. 
P g (t) = Electrical power dissipated in the cavity at time t 
P r (t) = Incident radiative power absorbed by the cavity at time t 
C = Cavity heat capacity 
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The rate of change of cavity temperature depends on the difference between 
incoming and outgoing powers. If we define the incoming power at time t by 


P in (t) = P r (t) + P e (t) ’ (4) 

then the power leaving the cavity at time t is given by 

P ft) = e A CT (t) 4 + e A a [T (t) 4 - T 4 ] + K. (T (t) - T ). (5) 

out a a c ss c s 1 c s 

Since the cavity temperature will never differ from T q by more than a 

small fraction of a degree, a linearization of equation (5) introduces 

completely insignificant error. The result of the linearization is 


P „(t) = P + K(T (t) - T ) 
out v ' out c o 


( 6 ) 


where we have defined 
P 


out 


= oT 4 [e A +e A ] - e A oT 4 + L (T - T ) , 
oaa ss sss lo s 


K = 4oT [e A + e A ] + K. 
oaa s s 1 


(7) 

( 8 ) 


Conservation of energy requires that the rate of energy stored in the 

cavity satisfy 
dT (t) 


dt 


= P. (t) - P _(t) 
xn out 


which can be rewritten, using equations (4) and (6) , as 
dT (t) 

= P r (t) + P e (t) - P° ut - K(T c (t) - T o ). 


dt 


(9) 


( 10 ) 


Proceeding any further in the analysis requires characterization of the 

servo, i.e. a defined relationship between P (t) and T (t). The desired 
• e c 


relationship is given by 


P e (t) = P e ° - r [T c (t) - T q ] (11) 

where P° is the electrical power required to keep the cavity at T q while 
viewing space. The parameter T characterizes the servo gain in units of 
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power per unit temperature change. The ideal nature of the stated relation- 
ship will become clear when the final solution for P (t) is obtained. The 

e 

electronics required to obtain the relationship defined by equation (11) are 
described in the following subsection. 

Defining T(t) = T c (t) _ T , and substituting (11) in equation (10) yields 
the following differential equation for T(t), i.e. 


C d ^ t -- + (K + T) T(t) = P (t) + P° - P° . 
dt r e out 


( 12 ) 


Since T, dT/dt, and P^(t) are all zero when the radiometer response to space 


has reached equilibrium, we conclude that 

P° = P° . 
e out 

As a result, equation (12) takes on the simpler form 

C + (K + T) T (t) = P (t). 

dt r 

The general solution to the equation (14) is given by the following 
expression 


(13) 


(14) 


T(t) = T (o) + tT ( o) (1 

where the time constant t is defined by 
C 


dP (f) 


-t/T. . T f “V' 

)+ W dt 7 


[1 - e - (t - t,)/T ]dt' 


T = 


K+r 


(15) 


(16) 


The first term of equation (15) represents the initial condition, the second 
term the decay from the initial rate of change at T at t = 0, and the third term 
the response to the time dependent radiative power incident from 0 to t. 

For the special case of a step function change in P^(t) at t = 0 from 
P to P the temperature response is given by 


P <2)_ p (1) 

T(t) = T (o) + r ~ [1 - e _t/ * ] . 


( 17 ) 
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The corresponding time dependence of the electrical power is just 


P e (t) ' P e (1) ‘ST t p r (2> -r r C1) l [l-e' t/T l. 


(18) 


At equilibrium we find that 


P <« - p <» -- [p (2 >-p (1 >] [i-s,. 
e e ~ r r r 


(19) 


For K/F << 1 the change in electrical power is equal and opposite to the 

change in absorbed radiative power. Since typical values of K/T are of 
-3 

the order of 10 (see next section), the electrical power difference is 
only about 0.1% less than the radiative power difference. Since this is a 
correctable effect, the actual error introduced is very much less than 0.1%. 

Thus we see that the transient response of the radiometer can be simply 
characterized by a single time constant t. This is an important feature 
resulting from the servo relationship stated in equation (11) . Servo 
systems producing power changes which are not linear in the cavity 
temperature difference would lead to non-linear differential equations, 
the solutions of which would in general contain a mixture of short and long 
time constant effects and, in some cases, resonances. The servo system 
design concept which satisfies equation (11) is described in the following 
subsection. 


4. Linear Servo Analysis 

The electronic servo system which produces the desired linear relation- 
ship between electric power and cavity temperature, as stated in equation 
(11), is based on a standard Wheatstone bridge circuit. A block diagram of 
the electronics is presented in Figure VI. 2. The operation of the servo 



V 

a 


HIGH GAIN 
PREAMP 



FIGURE VI. 2 


■ Servo Electronics 
Diagram 
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can be roughly explained by considering the response of the system to a 
decrease in irradiance of the cavity. The sequence of events which follow 
are listed below: 

(1) The cavity begins to cool, lowering its resistance R^ . 

(2) V increases due to the R decrease, thus increasing the imbalance 

cl X 

voltage V -V^* 

(3) The high gain preamp amplifies the imbalance voltage to the level 

g< VV • 

(4) The amplified imbalance voltage undergoes several analog operations 
which ultimately result in an increase in the bridge supply voltage 
to the level V = z[gy (V -V,)/V]^^. 

3 D 

(5) The increased supply voltage causes increased electrical dissipation 
in the cavity. 

(6) The cavity begins to heat as a result of the increased power 
dissipation, thus tending to increase the cavity resistance toward 
its initial value. 

(7) As the cavity resistance increases the voltage imbalance decreases 
eventually stabilizing the electrical power at its new equilibrium 
value. 

The function of the divider-multiplier circuit is to prevent the positive 
feedback effect which results from the imbalance voltage being proportional 
to the bridge supply voltage. Without this circuit it is possible to have 
resonance or very strong negative feedback, depending on the size of the 
preamplifier input offset voltage. The square rooter circuit is required 
to make the power linearly related to the cavity temperature. 

It is easy to mathematically verify that the servo satisfies equation 
(11) and thus behaves as described above. Since the bridge imbalance voltage 



is given by 


V ab (t) 


v a (t) - v b (t> - V(t) 1^5- 


l' w "2 

and since the bridge supply voltage satisfies 


V R 4 


1, 


( 20 ) 


V(t) = z[gy V ab (t)/V(t)] 


1/2 


( 21 ) 


the supply voltage is directly related to the resistance imbalance, i.e. 


V(t) = z[gy] lyf2 [ 


R. 


R 1 (t)+R 2 R 3 +R 4 


R 1/2 

4 3 • 


( 22 ) 


Note that R 2> R^> and R^ have been assumed constant, a result insured by 
using resistors with low temperature coefficients and proper coupling to 
a temperature controlled sink. 

In order to proceed further it i$ necessary to linearize the resistance 
imbalance relationship. Defining the cavity temperature coefficient of 
resistance a by the relationship 


R l (t) = R l [1 + aT(t)] 


(23) 


we can make the extremely accurate approximation 


^2 . R 2 m ° R l T(t> , (24) 

R 1 (t)+R 2 R° + R 2 1 R“ + R 2 J ’ 

Since the product aT(t) is typically smaller than 0.0001 errors in the 
linearization are completely negligible. 

Inserting (24) into equation (22) and squaring yields 


2 2 ^ R 2 R 3 R i R 4^® 


(25) 


where we define 
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Since the electrical power dissipation in the cavity is given by 

P e (t) = [I a (t)] 2 R 1 (t) = V 2 (t)R“ (R° + R 2 )- 2 , (27) 

where the approximation error is again negligible, we can obtain the 
desired result from equation (25), i.e. 


P (t) = p° - r T(t) , 
e e 

where we have used the notation 


(28) 


P° 

e 


2 < R 2 R 3 - R °lV R o 
82 Y 1 


(R° + R 2 ) 


r 


„ 2 
G z y 


R, 


(R° + 


R 2 } 


-2 


(29) 


(30) 


Note that equation (28) is identical to equation (11). 

In addition to the basic power temperature equation several of the 
servo relationships derived above, as well as several additional corollary 
relationships, are significant in optimizing the servo design. These are 
summarized in Table VI. 2. 


5 . Angular Response Variations Due to Cavity Non-Isothermalities 

Angular response deviations for a nearly spherical isothermal cavity 
have already been treated in Section VI. 2 and shown to be well within 
required limits for typical design parameters. However, the radiometer 
cavity is not strictly isothermal as a result of three distinct effects: 

(1) Variation in the angular distribution of incident flux; 

(2) Variation over the cavity surface of radiative transfer between 
the cavity and the heat sink; and 
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Table VI. 2. Summary of Servo Relationships 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


0 = r°/r 2 

r = 

Cl+0) *1 

3 -1 

v ab (t) = (R“P e (t)) 3/2 (l + |) (gz 2 y) 

P° = gz 2 y (R 3 -eR 4 )(l+g)“ 3 (R 3 +R 4 ) _1 (R“)' 1 S 2 

P,j,(t) = Total electrical power dissipated in cavity leg of the 

bridge 

P T (t) = (l+|)P £ (t) 

= (gz 2 yV ab (t)) 1/3 = (1 + -|) (R“P e (t)) 1/2 


(7) 


V(t) 
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(3) Non-uniform heat conduction resulting from support and lead 
wire contact with the cavity. 

If (1) were not true, the effects of (2) and (3) would be insignificant, 
as it is in the case of the ACR used for measuring the solar constant. 

The way in which (1), which does apply to the wide FOV ACR, can make (2) 
and (3) significant is best illustrated by an extreme example. Suppose 
the cavity is in contact with the heat sink at one point and thus has a 
fixed temperature there. Radiation striking the cavity a large distance 
away will result in a temperature gradient between the point of incidence 
and the contact point, thus raising the mean cavity temperature and causing 
the servo to reduce the electrical power input. However, if the same 
radiative power input occurs at the contact point, the mean cavity 
temperature will not be raised and the servo power will not change. Since 
different electrical power levels are produced for the same irradiance 
level, the measured powers are in error. 

In order to treat these effects numerically it is necessary to simplify 
the geometry of the cavity. Effects of wire conduction will be treated 
using a disc shaped "cavity", and non-uniform radiative transfer effects 
will be treated using a one dimensional "cavity". Although the results 
obtained with these simplified models are not exact, they do establish 
parameterizations which show what factors influence the error and also 
yield order of magnitude estimates for the size of the errors. The only 

significant uncertainty in using these model results is in transforming 
a geometrical coefficient to the case of the spherical cavity. (Accurate 

values for the coefficients can be obtained by measurement) . 
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The effects of non-uniform radiative transfer between the cavity and 
the sink are parameterized in Appendix I. The resulting errors can be 
kept well under control by appropriate cavity design. The effects are 
significant enough, however, to result in an important design constraint. 
How this constraint is used in the cavity design is discussed in Section 
VII. 

The parameterization of errors produced by wire and support conduction 
is developed in Appendix J. To first order these effects produce errors 
which are small and independent of cavity size and do not seriously 
constrain the cavity design. 
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VII. FAST ACTIVE CAVITY DESIGN AND ERROR ANALYSIS 
1. The Design Approach 

The specification of the FACR design requires the following cavity 
physical properties, electronics parameters, and operating conditions: 

(a.) Properties of the wire wound cavity shown in Figure VI. 1. 

i. Wire material 

p - resistivity 

a - temperature coefficient of resistivity 
p - density 

k - thermal conductivity 
c - specific heat 

ii. Wire size 

w - wire diameter 

iii. Wire wound cavity size parameters 

r - cavity maximum external radius 

A^ - cavity exterior area for unit cavity maximum radius 

r /r - ratio of the cavity aperture radius to the cavity maximum radius 

iv. Radiative coupling to the sink 

e T- effective emittance for radiative exchange 
s ' 

(b.) Servo parameters 

g - bridge resistance ratio as defined in Table VI. 2 

2 

gz y - servo gain, square rooter, and divider parameter product 
as defined in Section VI. 4. 

(c.) Operating conditions 

T - cavity temperature 
c 

T g - heat sink temperature 

P° - electrical power supplied to the cavity when viewing space 
(a super ° will be used throughout this section to denote 
the space view) . 
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The design approach was to first select a wire material and size for the 
basic cavity geometry shown in Figure VI. 1, to then apply a series of 
constraints imposed by the theory of the FACR operation and by the system 
measurement concept, and finally to select the design which was the most 
practical with respect to servo construction and power consumption. The 
wire material and size selected was 3 mil nickel. This choice was based 
on the desire to make the bridge signal voltage large without making the 
time constant too large. To do this a material with a large resistivity 
and temperature coefficient of resistivity was chosen since the bridge 
signal voltage is proportional to the cavity resistance to the 3/2 power 
and the time constant is inversely proportional to the temperature co- 
efficient of resistivity (See Table VI. 2). The values of the pertinent 
properties of nickel which were used for design considerations are the 
following: 

= 6.55 x 10 ^ft-cm a t o°C 
a = .0064 at 0°C 

3 

p = 8.9 g/cm 
k = .578 W/cm-°K 
c = .481 J/g-°K 

2 . Design Constraints 

Constraints on the values of the FACR parameters listed in subsection 1 

i 

are imposed by the theoretical expressions for the time constant, the bridge 
signal voltage, the radiative power loss of the cavity, the non-uniformity of 
the angular response, and also by the desired cavity emittance enhancement. 
Using the same notation as Section VI, they are as follows: 



The FACR time constant t = C/T can be expressed in terms of the 
parameters listed in subsection 1. The required time constant 
was found to be 48 ms in Section V. 

A bridge signal voltage AV which is much larger than the 

aD 

amplifier noise is required. The signal voltage is defined by 

V, = V , - V where V ° is the total imbalance voltage for 
ab ab ab ab 6 

the space view and V ^ is the imbalance voltage for the source 
view. A signal voltage AV ^ of 250 yV was selected for an 
irradiance at the cavity aperture corresponding to 1 solar 
constant. The resulting signal to noise ratio will be 
discussed later in this section. 

The electrical power supplied to the cavity is related to 

the cavity and heat sink temperatures by the expression for 

the radiative losses of the cavity to the heat sink and to 

space. This relationship is given by 

4 4 4 

P = Aeo(T -T ) + A e oT . 
e ssc s aac 

The first term represents cavity-heat sink exchange of radiation 

and the second the radiation to space. The conductive transfer 

\ 

of power to the sink was assumed to be negligible. 

The non-uniformity of the angular response due to the variation 
over the cavity surface of the effective emittance for cavity- 
sink radiative exchange e g should not limit the accuracy of the 
measurements. The fractional error f in the measurement of the 
source power absorbed by the cavity caused by an emittance 
variation of Ae can be approximated using the results of 


Appendix I to yield. 
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f = 


A Ae aT 
s s c 

6kw 


( 2 ) 


where A is the cavity exterior area (A = A, r ) . This fraction 
s si 

-3 

f was chosen to be 5 x 10 which according to Appendix H corresponds 
to a .1% error for a scene with extreme limb darkening or limb 
brightening. 

The variation of £ g over the surface of the cavity occurs 
because of the variation of the relative cavity-heat sink geometry. 

A detailed analysis of the variation of £ g over the cavity has not 
been performed. However, because of the nearly spherical shape of 
the cavity and heat sink a reasonable estimate of e g and its 
variations can be obtained from a simple result for concentric 
spheres. The effective emittance for radiative exchange between 
concentric spheres with diffusely reflecting and emitting surface 
coatings is given by 


e l £ 2 


'eff e 2 + e^Cl-e^) a^/ 


(3) 


where anc * are t ^ ie surface emittance of the inner 

(outer) sphere and the area of the inner (outer) sphere respectively. 

For e^ = e 2 = *93 and an area ratio of .70, which is reasonable for 
the FACR, the effective emittance is 0.89. The major deviation of 
from this result will occur at the apex of the FACR cavity where 
an area element on the sink has less of its field of view filled by 
the cavity. As a result the effective emittance is larger there be- 
cause less of the cavity radiation can be returned by diffuse reflection 
Equation (3) shows that e e ££ is always less than c for any finite area 

ratio so a reasonable estimate of the variation of £ is less than 

s 

0.04. For design purposes we have used a Ae of 0.03 which is 
considered a conservative value because the area of the cavity 
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where is large is small compared to the area irradiated by 

radiation from a given direction. A value of 0.90 was assumed 

for e . 
s 

(e) The ratio of the cavity aperture radius to the maximum radius, 

r^/r should be sufficiently small to allow accurate determination 
of the apparent cavity emittance from measurements of its surface 
coating emittance and also to make the cavity emittance insensitive 
to variations of the coating emittance due to its spectral response 
characteristics or degradation. A ratio of 0.5 was selected based 
on the results shown in Table VI. 1 which show that for this choice 
the apparent emittance uncertainty is a factor of 12 smaller than 
the uncertainty in the coating emittance. 

3. The Design 

The final design was determined from the theory of the FACR operation 
by applying the constraints discussed in the last subsection for a cavity 
constructed from 3 mil nickel wire. The geometrical design of the cavity 
as shown in Figure VI. 1 was represented by choosing the cavity exterior 
area A^ for unit maximum radius to be 10.97, the effective emittance 
for radiative exchange between the cavity and sink to be 0.90, and 
the ratio of the viewing aperture radius to the cavity aperture radius 
to be 0.80. The resulting FACR design parameters were determined for a 
number of different cavity and sink operating temperatures. 

It was found that the power required to operate the cavity decreases 

as the sink temperature T decreases and also as the cavity-sink temperature 

s 


i 
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difference T - T decreases. However, there are reasons for not making 
c s 

T and T - T too small. First the power required by the sink temperature 
s c s 

control servo increases as T g approaches the maximum spacecraft temperature. 

Second, it is clear that the operation of the cavity control servo depends 

on T -T being large enough that the electrical power supplied to the cavity 
c s 

to maintain the cavity-sink temperature difference while viewing space is 

larger than the power absorbed by the cavity for any source irradiance. 

The choice of T and T was not optimized, but values of T = 345°K and 
c s c 

T g = 320°K were selected as reasonable with respect to the resulting power 
requirements and the above considerations. The resulting FACR design 
parameter values are shown in Table VII. 1. The corresponding operating 
characteristics are summarized in Table VII. 2. 

The noise characteristics of the FACR were analyzed assuming use of 
the AC sense - DC heat servo discussed in Appendix N. The major advantage 
of this servo design is the noise reduction obtained by placing a trans- 
former between the bridge output and the preamplifier. In the DC servo 
case the dominant source of noise is the preamplifier. However, through 
the use of a transformer in the AC sense case, both the signal and the 
noise are amplified before reaching the preamplifier and the Johnson 
noise of the bridge dominates. The RMS noise voltage V n at the bridge 
is about 20 nV (see Appendix N) . 

The equations for the AC sense - DC heat servo corresponding to those 

summarized in Table VI. 2 for the DC servo are given in Table VII. 3. Except 

for equations (3) and (7) these relationships are the same as those for the 
2 2 2 

DC servo with gz y replaced by z A^ nmA c V^,. The signal to noise ratio was 
calculated as follows using the DC servo design value of 3.4 x 10 ^ for 
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Table VII. 1. Fast ACR Design Parameter 

Values 

Cavity Maximum Radius (r) 

1.31 cm 

2 

Viewing Aperture Area (ir(.8r) /4) 

2 

0.865 cm 

2 

Cavity Aperture Area (ffr /4) 

1.352 cm 2 

2 

Cavity Exterior Area (A^r ) 

2 

18.88 cm 

Cavity Apparent Emittance 
(assumed surface coating emittance of 
.90) 

0.993 

Cavity Heat Capacity 

3 mil Ni wire 0.48 w-s/°K 
paint 0.48 w-s/°K 

0.96 W-s/°K 

Cavity Resistance at Operating 
Temperature 

520 

Cavity Temperature 

345 °K 

Sink Temperature 

320°K 

Ratio of Cavity Space View 
Resistance to the other resistance 
in the cavity leg of the bridge (8) 

1.131 

Servo Power Gain Parameter (T) 

20.2 W/°K 

2 

Servo Parameter z gy 

/r c -i i \ 

3.41 x 10 7 


(See Section VI) 
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Table VII. 2. 


Fast ACR Operating Characteristics 


Time Constant 


48 ms 


Total Bridge Power for the Space 
View (P°) 

Electrical Power to the Cavity for 
the Space View (P°) 

Minimum Electrical Power to the 
Cavity (Normal Incident Solar 
View) 


0.80 W 
0.42 W 


0.31 W 


Per Cent Variation in Angular 

Response due to Variations of the 
Effective Emittance for Radiative 0.5% 
Exchange between the Cavity and 
the Sink 


Bridge Signal Voltage AV^ for 

1 Solar Constant Incident Flux 
using the DC Servo 


250 mV 
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Table VII. 3. Summary of AC Sense - DC Heat Servo Relationships 


( 1 ) 0 = 


R ! /R 2 


(2) r = 


(1+3) 


a[z 2 A p 2 nm A c V AC ] 

r: 


where, 

V ^ = AC sense voltage applied to the bridge (constant) 
n = Transformer turns ratio 

m = Demodulator Gain 


A^ = Preamplifier voltage gain 

Ap = Power amplifier voltage gain 

z = Gain of analog square rooter 

(1+1/0 ) 2 R° P° C (t) V 
(3) \T^(t) = 1 e AC 


(A) P 


,DC 


[z 2 A 2 nm A V.„] 

1 p c AC 

= [z 2 A p 2 nm A c V AC ] ^-gR^ (R^) - 1 R° - 1 (l+ 0 ) _ V 


DC 


(5) P T (t) = Total DC electrical power dissipated in cavity leg 


DC 


of the bridge 


DC . 


( 6 ) P T (t) = (1+1/S) P^(t) 


r 2. 2 ... , AC 

[z A nm A V. _JV , 1/9 „„ 

(7) V DC (t) = [ P- --- C -. AC , ab jl/2 = (1+1/g)[R opDC (t)] 

V AC X 


1/2 
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2 2 

z A nmA V and an AC sense voltage V of 6V. 
P C AL AL 


AV 


Signal to noise = 


AC 

ab 


n 


where AV ab is the bridge signal voltage defined as the change of the 
bridge AC imbalance voltage for a given ir radiance input to the cavity. 


INCIDENCE IRRADIANCE mW/cm 


& 

ab 


SIGNAL TO NOISE RATIO 


137 (1 solar constant) 

35 (typical irradiance for 
the total radiation sensor) 

21 (typical irradiance for 
the shor.t wave sensor) 


38.4 yV 
10.2 yV 

6.0 yV 


1920 

510 

300 


The noise equivalent flux NEF is probably a more useful quantity than 
the signal to noise ratio because it is essentially independent of the 
input irradiance value. It can be calculated from the equation. 


NEF = 


[(6P DC ) 2 + (SP°) 2 ] 1/2 


A e 
a a 


„DC 


where the cavity heating power variation P g caused by a noise voltage 
V n on the bridge imbalance voltage is given by: 

3P DC P DC 

= ( — — ) v = — V 
e n yAC n 


ab 


ab 


(4) 


(5) 


Equation (3) of Table VII. 3 was used to evaluate the partial derivative. The 

2 

noise equivalent flux for V^ = 20 nV is approximately .1 mW/cm . Although 

2 

this NEF is .5% of the typical short wave flux of 21 mW/cm , many 
samples of the same region of the earth obtained by the 6 to 8 satellite 
system will be averaged either before transmission to the ground or 
during data processing. The averaging process, which is inherent to 
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this mission, will make the effects of noise unimportant. 


4. Error Analysis 


An error analysis of the FACR was performed using the same technique 
explained in Appendix C for the solar constant ACR. The basic operating 
equation for the FACR can be written as: 


F = 


(P° - P ) [1 + K/r] 
e e 

A e 
a a 


where A is the aperture area, e is the apparent absorptance of the 

Si SL 

aperture and K is defined in Equation (8) of Section VI. The ratio 

-4 

K/r is very small (8.4 x 10 ) and is known to better than 5% so the 

factor [1 + K/r] introduces no significant error. The electrical 
power P^ supplied to the cavity is determined from two voltage measure- 
ments as follows: 


P = 
e 


(v-v 2 )v 2 


( 6 ) 


(7) 


where V is the total voltage across the bridge circuit and is the 
voltage across the resistance R 2 in series with the cavity resistance 
in the current carrying leg of the bridge. 

The measurement uncertainty of F was calculated for three different 
input fluxes by root sum squaring the weighted standard deviations of 
the parameters involved in equations (6) and (7). The standard deviations 
assumed for each parameter and their contribution to the total uncertainty 
are shown in Tables VII. 4, 5 and 6. As was the case for the solar constant 
ACR, the uncertainty in the cavity absorptance is the largest source of 
error. Therefore, most of the uncertainty is an absolute not a random 
error. The effects of cavity temperature variations were found to be 



VII- 12 


negligible. The effects of sink temperature variations are also 

negligible when a reference view of space is available on each rotation 

of the spacecraft. However, for certain orbits a space view may not be 

available for 10 minute periods where the sun is near the anti-nadir 

position. For these cases sink temperature stability relies on the 

sink control servo which allows temperature variations of + 1 m°K (the 

same sink control servo is used for the FACR as for the solar constant 

radiometer). This temperature uncertainty translates into a flux 

2 

uncertainty of .02 mW/cm which should be root sum squared with the 
results of Tables VII. 4 , 5 and 6 to get the results when a space view is 
not available for each spacecraft spin. 

The measurement uncertainty calculations presented in Tables VII. 4, 5 
and 6 assume no use of the solar constant radiometer for calibration 
of the FACR. Incorporating the solar constant radiometer measurements 
into the calibration of both the TR and SW sensors will probably reduce 
the measurement uncertainties from those presented here. 
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Table VII. 4. 

FACE Uncertainty with 1 Solar Constant InDut 

Parameter 

Nominal 

Value 

Standard 
Deviation (a.) 

l 

|o. (mW/cm ) 
1 

e 

a 

0.993 

0.003 

0 . 41 

A 

a 

0.865 cm 2 

0.0007 cm 2 

0.11^ 

R 2 

460Q 

0. 1ft 

0.03 

V 

23.8 V 

1 mV 

0. 03 

v° 

28.0 V 

1 mV 

0. 03 

V 2 

11. 2 V 

1 mV 

0.01 

V 2 

13.1 V 

1 mV 

0.01 



2 2 1 1/2 
[f (ff) ai ] 

1 9f i 

0.43 mW /cm 



Percent Uncertainty 

0. 31% 

Table VII. 5. 

FACR Uncertainty for Typical Total Wave 

Sensor 


2 

Input (35 mW/cm ) 



Parameter (£.) 

i 

Nominal 

Value 

Standard 1 

Deviation (o^) ■ 

|a i (mW/cm ) 
^i 

e a 

0.993 

0.003 

0.106 

A 

a 

0.865 cm 2 

0.0007 

0.028 

R 2 

460ft 

0. 1ft 

0 . 008 

V 

27.0 V 

1 mV 

0. 032 

v° 

28.0 V 

1 mV 

0.033 

V 2 

12.7V 

1 mV 

0. 004 

V 2 

13.1V 

1 mV 

0.004 



3F 2 2 1/2 

[T (§ )o , 2 ] 

X 1 

2 

0.11 mW/cm 


Percent Uncertainty 


.32% 



Table VII. 6. 


FACR Uncertainty for Typical Short Wave Sensor 
Input (21 mW/cnr) 


Parameter 

Nominal 

Value 

Standard 
Deviation (o^) 

| g r | o (mW / cm ) 
’i 

e 

a 

0.993 


0.003 

0.063 

A 

a 

0.865 cm 2 


0.0007 

0.017 

R 2 

460H 


0 . IQ 

0.008 

V 

27.4V 


1 mV 

0.033 

v° 

28.0V 


1 mV 

0.033 

CM 

> 

12 . 9V 


1 mV 

0.004 

v° 

2 

13 . IV 

[Z 

i 

1 mV 

8F 2 1/2 
Ovj 

0. 004 

= 0.08 mW / cm 2 


Percent Uncertainty .38% 



VIII-1 


VIII. THERMOPILE DETECTORS AS AN ALTERNATIVE 

Except for apparent deficiencies in spectral response, angular 
response, and possibly time constant, thermopile detectors offer a 
number of advantages over the fast active cavity for use in the wide 
FOV sensors. Among these are: 

(1) high signal to noise; 

(2) relatively simple associated electronics; 

(3) wide range of operating temperature; and 

(4) small and rugged construction. 

Since most of the limitations of thermopiles depend on absorbing 
characteristics of the black coating on its radiation receiver, they 
must be considered for possible application pending more complete 
analysis of the coating problem. The purpose of the present section 
is to describe the characteristics of thermopiles, their special 
requirements, and their overall capabilities compared to the fast 
active cavity. 

1. Thermopile Characteristics 

A wide variety of thermopile detectors are available commerically . 
Among these we are mainly concerned with thin-film devices because of 
the time constant constraints for the spinning sensors. (Most wire 
wound thermopiles have time constants measured in seconds instead of 
milliseconds) . 

Table VIII. 1 summarizes the significant performance characteristics 
of devices which have been considered because of their relatively short 
time constants. Note that there are two time constants listed for each 
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Table VIII. 1. Thln-Film Thermopile Characteristics (SI = Sensors Inc. ; 
BEC = Barnes Engineering Corporation ) 



SI Cl 

BEC 2-214 

BEC 2-215 

Target Size 

1 mm diam 

1 m x 8 mm 

4.7 x 5.6 m 

No. of Junctions 

12 

120 

308 

Operating Temp. 

-65°C - 125°C 

-60°C - 125°C 

-60°C - 125° 

D* (500°K, DC) 
(cm Hz 1/2/W) 

5 x 10 8 

0.5 x 10 8 

0.5 x 10 8 

Responsivity 
(500°K, DC) 

15 Volts /W 

1 Volt/W 

1.1 Volt/W 

Time Constant* 
with vapor 
deposited coating 

~20 ms 

25 ms 

30 ms 

Time Constant with 
3_M Black Coating** 

~100 ms 

~50 ms 



* All time constants are stated for vacuum conditions 

** Time constant estimates for 3M-Black coatings are based on telephone 
conversation with manufacturer representatives 
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detector. The short time constant is obtained using a thin vapor 
deposited black coating. As indicated in Figures VIII. 1 and VIII. 2 
the spectral response of the detectors using this thin coating is 
rather poor compared to 3M-Black or Parson’s black which are themselves 
marginal for accurate integrated flux measurements. The second set of 
time constants in Table VIII. 1 are obtained for the same sensor using 
a coating of 3M-Black instead of the thin film coating. The penalty 
for the improved spectral flatness is seen to be an increased time 
constant. The added heat capacity produced by the more massive 3M 
coating can probably be compensated for by reducing the thermal resistance 
between the radiation receiver and the heat sink (See Appendix K). The 
only penalty in this case would be a drop in D* and responsivity, both of 
which are well beyond requirements anyway. Thus we expect that thermopile 
detectors with short time constants (~20 to 40 ms) and relatively flat 
spectral response can be obtained, although we have not yet seen any as 
standard products. 

The output characteristics of thermopile detectors when mounted in 
the wide FOV sensors is derived in Appendix K . The basic equations for 
the SW sensor and the TR sensor are somewhat different due to the emission 


from the window used for the SW sensor. Using the parameters 


= window transmission at wavelength X 
li'sw ^ = incident short wave flux at wavelength X 
F^ = long wave flux emitted by the aperture 

F„ = long wave flux received from window emission and reflection 
W 



( 1 ) 



RELATIVE RESPONSE (PERCENT) RELATIVE RESPONSE 
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FREQUENCY (Hz) 


FIGURE VIII. 1. Performance Characteristics of Sensors 

Incorporated series Cl Thermopile detectors 
(Technical Bulletin 10, Sensors Inc.) 
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FIGURE VIII. 2. Relative Spectral Response of Barnes Thin - Film Thermopile Detectors 
(Bulletin 2-212, Barnes Engineering Corporation) 
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where e is the average long wave emissivity defined in Appendix K and 

J_iW 

where 


T ■ f\ F sv,x din sw 

J o 

oo 

e SW = /^A T A F SW,X dA/TF SW 
J Q 


( 2 ) 

(3) 


F sw ■ / F sw,i di - <4 

O 

Since both the aperture and the window which are in contact with the heat 
sink will be at very nearly the same temperature as the sink, the term in 
(1) involving e will be very small. However, in view of the possibility 

-LiW 

of slight deviations of the window temperature from we will retain this 
term in the analysis. 

The voltage output of the total radiation sensor is given by the 
simpler expression 


V TR k 2 [£ SW F SW + E LW [F LW + F A ° T s ]] 


(5) 


where k^, £^ w and e^ w are not identical to the same parameters in equation 
(1), because slight variations in construction result in differences in k^, 
and e from one sensor to the other. 

Figure VIII. 3 displays the general character of the output as a function 
of incident flux for both SW and TR sensors. For illustrative purposes we 
consider two cases for the TR sensor: (1) incident flux is all short wave; 

and (2) incident flux is all long wave. 


2. In-Flight Calibration of the SWS 

In condensed form, we can write the equation relating the SWS thermopile 
voltage output to sensor input as follows: 
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Figure VIII. 3. Voltage output characteristics for SW and TRS sensors. 
The voltage intercept of the SW sensor could be positive or 
negative, depending on the relative temperature difference between 
the window and the heat sink. 
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V = R F + V 
SW SW o,SW 


( 6 ) 


where the calculated value of the responsivity R and the voltage intercepts 

V „ T , can be found by comparing (6) with (1). It is the object of in-flight 
o , SW 

calibration to determine these parameters on a regular basis so that deviations 
from ground measurements and theoretical calculations can be tracked and 
accounted for as the sensor changes its characteristics in orbit. 

To determine R and V we propose to use voltage outputs measured 

O 9 o W 

when the sensor views space and when it views the sun. Assuming that the 
flux input from space is zero (a good assumption) and that the solar flux 
is known (it will be periodically measured by on-board absolute radiometers 
to an estimated accuracy of better than 0.2%) we obtain the following system 
of equations for the desired parameters: 

(7) 


v sw (sm) - R F(sm0 + Vsw 

V SW (SPACE) - V o _ SH 

Subtracting (8) from (7) we obtain the responsivity value 
R - [V SW (SUN) - V sh (SPACE) ] . 


( 8 ) 


( 8 ) 


The same technique would be used for the fast active cavitv to account for 
changes in window transmission. 


3. In-Flight Calibration of the TRS 

The total radiation sensor output as indicated in equation (5) depends 
on both Fg W and F . The condensed expression for its voltage output 

V TR = R SW F SW + R LW F LW + V o,TR 

shows three constants to be determined by calibration: the two responsivities 

R , R ; and the offset V . It should be noted that for detector coatings 
bW LW O y 1 K 
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with perfectly flat spectral response, which is not the case here, 

R„„ = R,.,,. The short wave responsivity R„„ and the offset V _ 

SW LW SW o , TR 

can be determined from space and sun measurements the same way R and 

V were determined, 

o , SW 

At first sight it would appear that determining R^ requires an 
on-board blackbody reference to provide the second point in a two point 
calibration process. The usual way to implement this second calibration 
point is to rotate the sensor so that it views an on-board blackbody 
cavity instead of the earth. Since this involves the complexities of 
moving parts, changes in thermal gradients, and changes in frequency 
content of the flux input, it conflicts with the system objectives of 
simplicity, reliability, and long life. An alternative method, which 
uses a temperature variable aperture and no moving parts, appears to 
be feasible. 

Rewriting the TRS voltage output response to a long wave input 
flux F as 

V - V - B, (10) 

the specific form of the background B is given by 

B - \w ^s 4 - r a oT a 4 - (11) 

where T = sensor temperature and T. = aperture temperature. The 
S A 

responsivity factor R^ is dependent on the long wave emissivity, the 
number of thermopile junctions, dV/dT per junction, and the conductivity 
of the thermal resistor (See Appendix K). The responsivity factor R^ 
differs from R^ by a geometrical factor accounting for the cosine-weighted 
solid angle of the aperture and an emissivity connection factor accounting 
for the fall-off near grazing incidence. 
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Two space views at two different aperture temperatures are used to 
determine R^ and R^. Voltage outputs for the two space views are given 


by 

V 1 * V T aIi - hw 0T S 4 

V 2 ' V T a! 2 - *U* 0I S 4 - 

It has been assumed that an accurate temperature control servo keeps 
constant for both views. The solution for R^ and R^ in terms of the 
measured voltages and temperatures are given by 

r a - <w wt a * 2 - t a *p r 1 . 

\w " <T A,1 V 2~ T A,2 V l )[aT s <T A,2 “ T A,1» ' 


( 12 ) 

(13) 


(14) 

(13) 


In order to estimate the effects of measurement uncertainties in 
the calibration on the derived flux values we shall first invert equation 
(10) to express the incident long wave flux F in terms of the calibration 
parameters, i.e. 


F -\w [V+o( \ W T s 4 - R A T A 4)! - 

where in this case T and T will be very nearly equal. The standard 

A S 


( 16 ) 


error in the flux (a ) is related to the errors in the measured parameters 

F 


x^, k = 1,7 by the equation 

• , I ,3F , 2 2,1/2 

°F { k=l ( 3x k ) ° x k } 


( 17 ) 


where x , ...x = V ,V_,T ,>T A „,T ,T ,V. The derivatives are determined 

J. / X Z A| X A)Z S A 

from equations (14), (15) and (16). 

For the typical parameter values listed below, the value of o can 

r 


be estimated 
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*lw = 15 V/VH 


M 


-2 


r a - °- 2 

T = T. = T = 270°K 
s A A,1 

T A = 320°K 
A, 2 


aT = 0. 05°K 
s 


a T. = 0. 10°K 
A 


°v = lyV 


Inserting these values into equation (17) yields a standard flux err*or of 
2 

0.3 W/M which is within the required limits presented in Section II. 


4. Deriving Long Wave and Short Wave Fluxes from Thermopile Sensor Outputs 


For a single measurement the short wave flux from the earth is obtained 


by inverting equation (6), i.e. 

F = — (V - V ) . 

SW R v SW o,SW' 


(18) 


The long wave flux emitted by the earth is determined by inverting equation 
(9) and inserting (18), i.e. 


1 R SW 

F = (V - V ) — F 

LW R lw v TR o,TR' R lw SW 


(19) 


Similar equations would apply to the fast active cavity option, except that the 
short wave and long wave responsivities of the TRS will be known absolutely. 

5. Radiometric Noise Estimates for Thermopile Sensors 


The noise equivalent flux (NEF) can be determined from the following 
relationship 

NEF = y(Af. T ) 1/2 (A 1 ^ 2 D*)~ 1 , where 
N d 


( 20 ) 
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Af = the system noise equivalent, bandwidth 
= the detector area, 

D* = the specific detectivity of the thermopile, and 
y = the amplifier noise factor. 

In the case of thermopiles the noise power spectrum is flat; thus, assuming 
a similar noise spectrum for the electronics, the noise equivalent bandwidth 
is equal to the information bandwidth. In order to prevent reduction of the 
inherent time constant of the detector, we will assume an electronic band- 
pass of 0 - 500 Hz. Assuming a conservative amplifier noise factor of y = 2, 
the NEF's for the three TP detectors described in Table VIII. 1 are found to 
be 

NEF = 0.005 W M -2 (SI Cl) 

NEF = 0.03 W M -2 (BEC 2-214) 

NEF = 0.02 W M -2 (BEC 2-215) 

l 

all of which are well below the required limits (by at least an order of 
magnitude) . 
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IX. PRE-FLIGHT AND IN-FLIGHT CALIBRATION PLANS 


Although sensor calibration procedures have been discussed previously 
in Sections III, IV, VI and VIII, the importance of this topic warrants a 
separate and more detailed description. The system calibration discussed 
in the following subsections applies specifically to the preferred option 
which employs a slow ACR as the solar constant radiometer and two FACR' s 
for the wide field of view sensors. The special calibration problems of 
the thermopile alternative are discussed in Section VIII. 

1. Calibration Objectives 

The basic objective of any instrument calibration is to define the 
relationship between instrument measurements and fundamental laboratory 
standards. Successfully meeting this objective requires thorough under- 
standing of instrument behavior, carefully designed specific comparison 
measurements referenced to laboratory standards, and instrument stability 
between calibration and utilization. In the case of radiometers based on 
thermopile sensors, instrument behavior is not understood sufficiently to 
obviate the need for direct comparison with laboratory radiometric stan- 
dards. Given a linear instrument, the comparison with laboratory radio- 
metric standards is required to determine the constant of proportionality. 
The ERB calibration plan (Hickey, 1974) provides a complex mechanism for 
making this determination. The ERB long wave calibrations are based on 
blackbody cavity sources ultimately traceable to the IPTS 1968 temperature 
scale. The short wave calibrations are derived from comparison with 
reference radiometers using controlled radiation sources. The ERB flight 
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instruments are calibrated by comparison with a reference sensor model 
which is itself calibrated by comparison with normal incidence pyrheli- 
ometers which are in turn transfer standards referenced to a Kendall- 
type Primary Active Cavity Radiometer (PACRAD) which serves as the 
primary standard. 

Since the sensor system proposed in this document contains an 
absolute cavity radiometer standard (the solar constant ACR) as part 
of the flight package, the calibration plan in this case does not 
require laboratory comparisons with radiometric standards and, in 
general, is much simpler and more direct than that required by ERB. 

Both the Wisconsin Solar Constant ACR and the PACRAD referred to by 
Hickey are absolute radiometers in the sense that they derive their 
calibration mainly from absolute electrical standards and not from 
other radiometric standards. The theory of these types of electrical 
substitution cavity radiometers is well understood and experimentally 
verified (See Section IV); and their stability is exceptionally good, 
mainly limited by the stability of the power measurement electronics. 

The wide FOV sensors proposed for the Wisconsin system (the Short 
Wave FACR and the Total Radiation FACR) are also electrical substitution 
absolute cavity radiometers and will be calibrated in much the same way 
as the solar constant ACR. However, the cavity emissivity enhancement 
factor for these instruments (~12) is not as large as it is for the 
solar constant ACR (~50). Since this might lead to small changes in 
calibration resulting from surface coating changes, the solar constant 
ACR, which is four times less sensitive to this effect, will be used 
in-flight as a reference radiometer for updating the wide FOV radiometer 
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calibration using the sun as a source. 

2. Pre-Flight Calibration 

The pre-flight calibration of the solar constant ACR is established 
by a combination of laboratory measurements and theory. The several 
factors involved in this process are indicated in the block diagram of 
Figure IX. 1. Surface emittance measurements are used in conjunction 
with the theory of cavity emittance to establish the emittance of the 
cavity. Since the normal incidence emittance enhancement for this 
cavity is approximately 50, absolute surface emittance uncertainties 
of <5% are required. Aperture area measurement, power measurement 
electronics calibration, and small correction factors (see Section IV) 
complete the absolute calibration. The small correction factors are 
derived from the ACR theory and some auxiliary measurements (e.g. the 
approximate thermal resistance between the cavity and its associated 
heat sink). The two comparisons indicated in Figure IX. 1 serve two 
distinct functions. The PACRAD comparison serves as a quality control 
test. The Solar Constant ACR and the PACRAD should agree to within 
predicted uncertainties provided both function properly. The comparison 
with Angstrom Normal Incidence Pyrheliometers embodying the 1956 
International Pyrheliometric Scale is planned merely to establish the 
relative difference between that scale and the Solar Constant ACR. 

The FACR' s used for the wide FOV sensors will undergo pre-flight 
calibration which is structurally similar to that described for the 
Solar Constant ACR (Table IX. 1) . The major difference occurs in the 
case of the short wave sensor, which uses a Suprasil-W window. For 
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Figure IX. 1. Absolute Calibration Chart for the Solar Constant 
ACR. The comparisons are made for reference, not 
as part of the calibration. 
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this sensor an additional step, i.e. determination of the window trans- 
mission, is required to establish pre-flight calibration. Since both 
wide FOV sensors are also calibrated in-flight using the solar constant 
ACR their pre-flight calibrations are, to a degree, redundant. Such 
redundancy is, however, a desirable feature which can aid in rapid 
identification of possible degradations resulting from launch. 

3. In-Flight Calibration 

As indicated previously, in-flight calibration of the wide FOV FACR's 
is accomplished using the sun as a source and the solar constant radio- 
meters (ACR' s) as reference radiometers. Because of the constraints of 
the spinning spacecraft assumed in this study it is not possible to per- 
form these in-flight calibrations with a 'single spacecraft; both wide 
FOV sensors and the narrow FOV solar constant ACR on the same spacecraft 
cannot simultaneously view the sun. In the multiple spacecraft system 
proposed here the wide FOV sensors on S/C #1 use the reference ACR's on 
spacecrafts #2, #3, ..., n. This procedure is illustrated in Figure IX. 2. 
Assuming 2°/day orbital precession relative to the earth-sun line, each 
S/C will be in favorable sun viewing position for calibration of wide FOV 
sensors at approximately 90 day intervals. At each such opportunity 
(which lasts for the order of a week or two) any of the Solar Constant 
ACR's on the other spacecrafts of the system can be used to measure the 
solar flux which, combined with the simultaneous wide FOV response to 
the same solar flux, establishes the in-flight calibration of the wide 
FOV sensors. Generally only one or possibly two other ACR's will be 
used at any given opportunity since only this number of the other S/C's 
in the system will have attitudes near optimum for pointing their spin 



Calibration (Electrical Standards) 



Figure IX. 2. Wide FOV Sensor In-Flight Calibration using Solar Constant ACR’s on Multiple Spacecrafts. 

In-flight Calibrations could occur every few months for each set of wide FOV Sensors. 

For simplicity only the calibration sequence for sensors on S/C //I are shown in the diagram, 
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axis at the sun. Using all other Solar Constant ACR's is not only 
unnecessary but also inconvenient and inefficient since, at any 
given time, several of these S/C’s would require large attitude 
changes before and after each calibration event. 

Both the Solar Constant ACR and the wide FOV sensors actually do 
not measure absolute flux directly. Instead they .measure absolute 
flux differences. Consequently, an external flux reference is required 
to establish the offset. In both cases a view of space provides this 
reference, i.e. zero incident flux. Figure IX. 3 illustrates the role 
of the space reference in transferring calibration from ACR to FACR. 

4. In-Flight Diagnostic Measurements 

As already indicated, an important aspect of the process of 
calibration is thorough understanding of instrument behavior. The 
purpose of the in-flight diagnostic measurements is to monitor 
instrument characteristics so that changes from the theoretical and 
preflight behavior which can affect performance, yet not be observed 
easily in the transfer calibration methods, can be identified. Among 
these diagnostics the following list is selected for illustration of 
their significance and utilization: 

(1) Space power offset variation with time - this measurement 
of FACR behavior uses the normal sampling mode, reveals 
out of specification electronics drifts. 

(2) Intercomparison of wide FOV sensors (SWS and TRS FACR's) 
against the sun uses normal sampling mode, allows rapid 
discovery of sudden degradation of either of the wide 


FOV sensors. 


Incident solar flux magnitude 
determined from slope and power 
difference between space and 
sun views . 


Solar Constant ACR 
Slope determined by 
Absolute Calibration. 


Electrical power 
provided to cavity 


Solar View 
Power 

Measurement 


Space View 
Measurement 


(a) SOLAR CONSTANT ACR 


These two points determine FACR 
-slope providing in-flight calibration 


Electrical Power 
provided to Cavi 


Space View 
Power 

Measurement 


Space View 
Measurement 


(b) WIDE FOV FACR's 


Figure IX. 3. Calibration transfer between Solar Constant ACR and 
Wide FOV FACR's. 


Wide FOV sensor response to incident solar flux as a 
function of rotation angle - uses rapid sampling mode, 
establishes in-flight angular response characteristics, 
reveals non-uniform degradation of windows, degradation 
of black surfaces, and misalignments between cavity and 
aperture. 

Wide FOV sensor response to space (night measurement 
using SWS) - uses rapid sampling mode, yields a direct 
measure of window heating amplitude and phase, reveals . 
changes in window thermal coupling to sink. 
Intercomparison of ACR measurements - verifies relative 
error predictions, identifies out-of-specification ACR's. 



X. SYSTEM ERROR SUMMARY 


A large number of system error sources are discussed both 
in the body and -'.n the appendices of this report. In order to put the 
errors in a system perspective, representative values for the important 
errors are shown in Table X.l. The errors are simply compiled for easy 
reference: the table is not meant to represent a complete system error 
analysis which cannot be performed at this time. Such an analysis would 
require an optimized orbit set for prediction of the sampling errors 
which are crucial to the system performance (see Section II). Also, 
better data is required on black surface characteristics than is 
generally available. The error totals shown in the table are worst 
case estimates, since the errors were simply summed. 

The single sample errors are divided into those which depend on the 
scene viewed and those which do not. This was done to separate errors 
having different impacts on the system. The scene independent errors 
shown for the FACR are substantially larger than those for the TP sensor. 
However, the FACR single sample noise and the measurement uncertainty are 
errors which do not impair the system performance significantly. The 
single sample noise will be reduced to a negligible level by the required 
averaging of samples from designated earth areas, both on-board each 
spacecraft, and for data processing. The FACR measurement uncertainty is 
dominated by the small cavity absorptance uncertainty which is an absolute 
percent error. Since data from FACR' s on different spacecrafts can be 
intercompared this error constitutes a constant system percent error. Also 
the absolute error introduced by the cavity absorptance uncertainty can be 
reduced by intercomparisons of solar irradiance measurements made by the 
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solar constant ACR and by the total radiation and short wavelength FACR's. 
The relative error portion of the FACR measurement uncertainty and 
calibration is comparable to the calibration error given for the thermopile 
sensors . 

The scene dependent errors of the FACR indicates its substantial 
advantage over any instrument with a flat detecting surface. The errors 
introduced by non-uniformity of the detector response to angular and 
spectral variations of the scene irradiance are reduced from those of a 
flat plate by a factor 12.5 for the FACR design specified in this report. 
However, the expected thermopile errors indicate that a system of 8 
satellites employing TP ' s would have capabilities adequate to meet the 
requirements of this mission as stated in Table X.l and would be clearly 
superior to any past system used to measure the radiative energy budget 
of the earth. Which of the two wide FOV sensor options should be employed 
ultimately depends on surface coating emissivity characteristics and 
their in-orbit stability. Uncertainties in this area would heavily weight 
the balance in favor of the FACR option. 

The solar constant radiometer accuracy is fully discussed in 
Section IV and Appendix C. 
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XI. INTERFACE REQUIREMENTS 

The sensor design presented in this document has been developed for 
use in a fleet of spinning spacecrafts dedicated entirely to measurement 
of the earth's radiative energy budget. The mission objectives and 
corollary sensor requirements together determine spacecraft and orbital 
requirements. The purpose of this section is to describe the general 
requirements thus implied. Since the sensor system has not been designed 
as a bolt-on and plug-in package for riding on operational satellites, 
but is instead meant for integration in a dedicated spacecraft as the 
only observing instrument, specific configurations of subsystem components 
which would be determined mainly by dynamical considerations are not 
described. 

1. Thermal Interface Considerations 

All sensor subsystems are thermally controlled by servo systems and 
insulated from the S/C structure. However, there must be sufficient 
conductive coupling to maintain positive power flow from the sensor to 
the S/C under minimum temperature differences. If we define 

= the hottest and coldest S/C temperatures experienced 
during a full year of precessional cycles, 

= the sensor sink temperature, 

= the minimum power which must be conducted from sink 
to S/C, 

= the maximum power conducted from sink to S/C, and 
= the thermal conductance of the insulation between 


.00 (c) 

S/C’ S/C 


T 

s 

P .. 
mm 


P 

max 

K 


sink and S/C 
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then the thermal conductance must satisfy 
K ' P „in (T a - T s/C>' 1 - 


( 1 ) 


and the resulting maximum power is thus found to be 


max 


,(cK _ 

‘S/w min ' s 


P = K(T_ - T^') = P_,_ (T_ - T^) (T„ - 


, (h) s -1 

s l S/C ; 


( 2 ) 


For an integrated sensor S/C design T s> T g/c* wou -*-^ be system 


optimized as a group considering all constraints. From a bolt-on package 

point of view we will choose a fixed sink temperature T = 320 °K and 

s 


require 


,(h) 

L s/c 


T - T v “'_ > 40°C 


T Srt - T s/c s 20 ° c 


(3) 

(4) 


in which we obtain 


P <1.50 P . 
max ~ min 


(5) 


which will be used in estimating total power requirements. 


2 . Power 

Lacking detailed electronics subsystem design, power estimates are 
only approximate. Values listed below are based on comparison with 
results for in-house spaceflight hardware of similar function and 
complexity. 
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Power (watts) 


Subsystem 

W/0 CMOS 

W/CMOS 

Command & Control 

0.8 

0.4 

Data Processor & Memory 

1.0 

0.5 

Power Conversion 
(Converter Loss ~25%) 

2.0 

1.5 

Attitude & Spin Control 

0.65 

0.65 

Transmitter 

0.07(4.0) 

0.07(4.0) 

Receiver 

0.015 (1.0) 

0.015(1.0) 

Housekeeping 

0.1 

0.1 

Sensor Servo Control Circuitry 

1.0 

1.0 

FACR Sensor (2) 

1.6 

1.6 

FACR Sink Control Power (2) 

0.6 (0.8) 

0.6 (0.8) 

ACR Sensor 

.003 (0.25) 

.003 (0.25) 

ACR Sink Control Power 

0.004(0.4) 

0.004 (0.4) 

TOTALS 

7.84 

6.44 


The maximum powers for those subsystems with variable loads are 
shown in parentheses. The transmitter and receiver powers assume a 
duty cycle of 2 and 5 minutes every 12 hours respectively. The 
solar constant ACR will only be used four times per year. 

The power available from the solar panel array was estimated by 
R. Rassmussen of Electro Optical Systems to be 14.1 watts at the end 
of a 10 year mission. The following conditions were assumed: 

1. Spherical spacecraft with 24" diameter. 

2. Of the spacecraft surface, 90% is available for solar 
panels and 72% for solar cells. 


3. Circular orbit with an altitude of 380 n.mi. 
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4. Worst case illumination of 54.6 min day light and 35.6 min night. 

The resulting power margin assuming a battery efficiency of 0.8 is 4.84 watts 
if CMOS are used and 3.44 watts otherwise. The satellite power profile is 
illustrated in Figure XI. 1. 


3. Weight 


The weight estimate for a complete satellite is given below. The values 
are approximate but conservative and are based on comparisons with in-house 
spaceflight hardware of a similar function and complexity. 

Subsystem Estimated Maximum Weight (lbs) 

Sensors 


FACR (2) 

Solar Constant ACR 
Sensor Electronics 

Power 

Communication 
Cable Harness 
Structure 

Command and Control 

Attitude Control 

Data Processor 
TOTAL 

Two satellites of this weight 


6.0 

3.0 

2.0 

18.0 

5.0 

4.0 

15.5 

5.0 

5.0 

6.0 

69.5 lbs ~ 32 Kg. 

easily be launched together on a 


single scout launch vehicle. 
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Figure XI. 1. Satellite Power Profile 
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4. Signal Interface Requirements 

The data processor subsystem (see block diagram Figure XI. 2) requires 
the following timing signals: 

(1) Spin clock (frequency of S/C rotation x 1024) 

(2) SWS pointing at earth nadir 

(3) SWS pointing at min. relative sun angle 

(4) Read out clock frequency (to be determined) 

(5) Read envelope 

Timing signal errors (for (1) through (3)) are assumed to be directly 
dependent on attitude errors. Requirements for the latter are outlined 
in a separate subsection. The read clock and read envelope are used to 
shift serial digital data from processor memory to the S/C communication 
system for transmission to the ground station. Data storage of approximately 
40,000 bits requires ground read out twice per day. Additional commands 
required are 

(1) Read Command 

(2) Mode (normal, solar constant measurement, other) 

(3) Power on/off 

Detailed signal characteristics (pulse widths, rise times, logic levels, 
etc) are to be determined. Housekeeping data for sensor systems are encoded 
digitally in the 40 K bit data block and do not require analog signal lines. 

5. Attitude Requirements 

In order to keep attitude errors from affecting accuracy of reflected 
solar and emitted fluxes we must have the nadir pointing accuracy better 
than + 1.8°. The ACR used for measurement of the solar constant requires 



Command 
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Figure XI. 2. Sensor System Block Diagram 










XI- 8 


sun pointing accuracy of + 0.5°. It is also necessary to have relatively 
high attitude accuracy relative to the sun for calibration of the SWS and 
the short wave responsivity of the TRS. The latter requirement is deter- 
mined as a function of the minimum pointing angle between the sun and the 
sensors in Appendix L. A reasonable pointing accuracy in this case is 
+ 0.5° which would allow for wide FOV shortwave response determination 
error less than 0.2% for sensor plane sun line angles in the range + 30°. 

6. S/C Spin Requirements 

Spacecraft spin axis must be parallel to the orbit plane normal 
within 1° and known accurately enough to meet the attitude requirements 
outlined in the previous subsection. The chosen spin direction results 
in the S/C "rolling" along the orbit. The spin rate should be 
approximately 6 RPM. Spin rate accuracy is not a separate requirement 
provided timing signals provided are sufficient to meet attitude 
requirements. 

7. Orbital Requirements 

Optimum orbit sets have not been determined. The general requirements 
are for 6-8 satellites in approximately circular orbits 380 n.mi. in altitude. 
The altitude should be as low as possible within the constraint of the 
desired systems lifetime. A 380 n.mi. circular orbit has a 98% probability 
for a 6 year lifetime. Multiple satellites with precessing orbits are 
necessary to obtain sufficient diurnal sampling (see Section II) . Although 
three equally spaced satellites with precession rates of 2 ° /day relative to 
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the sun would probably provide adequate monthly averaged sampling for most 
longitudes, more satellites are necessary because no orbit adjust capability 
is employed to maintain the satellite spacings. Also, the inclusion of a 
polar orbit and a near equatorial orbit would improve the sampling. Sun- 
synchronous orbits are not desired because they severly limit diurnal 
sampling capabilities. 
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APPENDIX A. CALCULATION OF NORMAL AND HEMISPHERICAL CAVITY EMITTANCES 


The accuracy of the solar constant measured with an active cavity 
radiometer depends heavily on knowledge of the cavity absorptance. Methods 
of approximating and bounding cavity emittance (and therefore absorptance) 
are developed in this appendix. The results are sufficient to justify the 
uncertainty assumed for the ACR absorptance in Appendix C and Section VII. 4. 


!• Apparent Emittance of the Cavity Surface 

Assuming that the cavity surface is opaque and grey (p = 1-a = 1-e) 
that it emits and reflects diffusely, and that no radiant flux enters 
the cavity from the outside, the flux balance for a point x on the 
surface is given by the following integral equation (Sparrow et al. , 1963) 

B(x) = e a T 4 (x) + p jf B(y) dF x _ y (1) 

where 


e - surface emittance 
T(x) - surface temperature 

P - surface reflectance, p = 1 - a , a is the absorptance. 

B(x) - total emitted flux. 

4 

B(x) = e(x) a T (x) , e(x) is the apparent surface emittance. 
dF^_ y - view factor from position x to position y. It represents 
the fraction of the energy emitted from a small area da(x) 
which arrives at da(y) . 

The integral is performed over all y on the surface of the cavity. If 
the temperature of the cavity is uniform. Equation 1 reduces to: 

e(x) = e + p JC e(y) dF x _y* (2) 

To determine the apparent emittance e(x) for a given cavity directly 
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from Equation 2 requires a numerical computer solution. However, 
without obtaining such a solution for each possible cavity design, 
we can learn a great deal about the dependence of emittance on cavity 
design by examining Equation 2 more closely. 

First, it is clear from Equation 2 that e(x) is always greater 
than e since the integral portion of the equation is positive. Therefore, 
the crudest possible bounds on the apparent surface emittance are: 

e $ e(x) $ 1. 

Substituting these bounds into the RHS of Equation 2 gives an 
improved set of bounds : 


e + ep JdF $ e(x) $ e + p JdF 

j x-y K x-y 


k 


or 


e + ep(l-F )* e(x) $ e + p (1-F ) 

X U A U 


(3) 


where F = 1 - JdF is the view factor from x to the cavity aperture, 
x-o x-y J r 

Second, by writing the formal solution of Equation 2 as an infinite 
series in increasing powers of the reflectance, the bounds of Equation 3 
can be further improved. The formal solution is as follows: 

e(x) = e + pe JdF + p 2 e jfdF dF + ... 

x-y y-z x-y 

= e + pe(l-F ) + p 2 ef(l-F ) dF + ... 

x-o y-o x-y 

Now define F* to be the maximum of F for all x on the cavity. We 

x-o 


get a lower bound to e(x) as follows: 


e(x) £ e + p e (1-F*) = p^e (1-F*) fdF _ + ... 

x y 

= c + p e (1-F*) + P 2 e (1-F*) (1-F X _ 0 ) + ••• 


oo r\ n i 

>, e l P" ( 1-F*) n = 7 i i — ^ 

n=o 1-p ( 1-F* ) 



A-3 


An upper bound of the same form can easily be found. Then using 

Equation 2 the following bounds result: 

defining F F £ F„ 

1 x-o 2 


e. = - i=l,2 

1 l-p(l-F i ) 

e + pe r ( 1_F X _ 0 ) S e(x) $ e + p e 2 (1-F x _ q ). 

Third, for a spherical cavity with a uniform surface coating, the 
apparent surface emittance can easily be calculated. The view factor 
from one area element da(x) of a sphere to another da(y) is: 


dF 


x-y 


_ da(y) 
4 ttR 2 


for a sphere of radius R. This result is interesting because dF x _ y is 
not dependent on x. Therefore, Equation 2 shows that e(x) is a constant 
and is given by: 

e(x) = e + e(x)p f dF_ 


x-y 


or 


e(x) = 


1-p /dF 


x-y 


1-p (1-F ) 

K x-o 


where in terms of the area of the cavity A , 

c 


/dF _ .. = 1-F 


x-y 


x-o 


4itR 


2 ' 


Based on the result for the sperical cavity, we suspect that for nearly 
spherical cavities the apparent surface emittance will not vary a great 
deal over the surface. Therefore, Equation 2 can be approximated by 

e(x) = e + pi fdF 

J x-y 

= e + pe (1-F ) 

■ x-o 


(4) 


(5) 


( 6 ) 


(7) 


( 8 ) 



where e(x) inside the integral is approximated by its area average over 

the surface: 

~ _ Je(x) da(x) 

/ da(x) 

Now, we have developed bounds, Equation 3 and 4, and an approximation, 
Equation 8, from which we can proceed to obtain relations for the apparent 
hemispherical and normal emittance of the cavity as a whole. 


2. Cavity Hemisp herical Emittance 

The apparent hemispherical emittance is the ratio of the emitted 
flux to that which would be emitted by a perfect blackbody. That is: 


H 

e 


a 


i 


B(x) F (x) da(x) 
x-o 


/ a T (x) F x _ q (x) 


da(x) 


For a uniform temperature distribution, 

H / (x)F x- Q (X) da(x) 

E a / < F (x) da(x) 

J x-o 


■ V^‘ 


e(x) F x _ q (x) da(x) 


where A is the aperture area. 

Substituting e (x) from Equation 4 and Equation 8 gives the following 
results : 


(9) 


e + pe 1 (1-f) $ $ e + pc 2 (1-f) (10) 

= e + p e(l-f) (11) 

/F 2 da(x) 

J x-o 

TV da(x) 

J x-o 


where f = 
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A more useful approximation than Equation 11 can be obtained by letting 

U 

£ = e which is a weighted average of e over the surface. Then: 
a 

H ^ e 
£ a l-p(l-f)' 

Substituting 1-e for p, 

H & £ 
e a e +(l-£)f 

For a spherical cavity of radius R, this result is exact. We can see 
from Equation 7 and the relation ^^ x _ 0 da(x) = A that: 


f 


F 

x-o 





( 12 ) 


(13) 


(14) 


where A is the aperture plane area, A is the area of the cavity and r 
ap c 

is the aperture radius. The spherical approximation (Equation 13 with 
Equation 14) was compared with the emittance of three nearly spherical 
cavities formed from cyclinders and cones calculated numerically (Syndor , 1970) . 
The approximation was accurate to within . 2 % when a sphere of the same area 
as the nearly spherical cavity was used for the approximation. 

Finally, it should be noted that the lower bound of Equation 10 is 
quite useful. While it may appear difficult to determine e^, for most 
cavities is the emittance of a spherical cavity with the same aperture 
as the cavity of interest and with a radius such that it just fits inside. 


3. Cavity Normal Emittance 

The equation for apparent normal emittance is formally very similar 
to that for hemispherical emittance. However, instead of dealing with all 
of the flux emitted by the cavity, we consider only the flux emitted to 
an area equal to the aperture area and essentially an infinite distance 
away in the normal direction. The solid angle from the emitting surface 
to the distance area is independent of position on the cavity so the normal 



(15) 
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emittance is given by: 


N 

£ =, 
a 


Jb(x) n • d a n (x) 


I 


o T (x) n *d a (x) 
n 


where n is a unit vector in the normal direction, the area vector 
d a^(x) is normal to the surface at x and the integration ranges 
over the inner cavity surface a^ bounded by the normal projection of 
the cavity aperture boundary. It is clear that for a spherical cavity 
with a constant surface temperature (B(x) and e(x) are constants), the 
normal emittance is equal to the hemispherical emittance. 

For an isothermal cavity with n • d a^(x) independent of x (e.g. a cone) 


Equation 15 reduces to: 

N / e < x ) d a n < x ) 


£ = — - 
a 


fd a n (x) ' 

Again we can use the bounds of Equation 4 and the approximation of Equation 
8 to yield: 


(16) 


e a + pe l ^ 1-g ^ * E a * £ + pe 2 ^ 1_8 ^ (17) 

and 

= e + p e (1-g) (18) 

where g = K o (x) d a n (x)/J*d a n (x). 


The approximation £ used in Equation 18 has a form very similar to that 
for g. However, e is the average of e(x) over the whole cavity area, not 
just over a^. 

Now consider the following special cases: 

CASE 1. Right Circular Cone. 

Since for the cone the area a is the whole cavity area, 

n 


IF da (x) = A 
J x-o n v 

y d a (x) = A 
n 


cone 


ap 



and 
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Therefore, g = sin 0, where 0 is the half angle of the cone, and 

Equation 17 and 18 reduce to: 

N 

£ + pe..(l-sin0) < e $ e + pe„ (1-sin 0) 

1 cone 2 

N . e 

cone £ + (1 -e) sin 0 

CASE 2. Cone Plus Upper Cavity: 

For simplicity we assume that the height of 

the upper cavity is much larger than the radius 

of the aperture. Then: 
r 2 

g = (-) sin 0, 

and e N = e + p e (1 - Cr) 2 sin 0). 
a h 

Also, the normal emittance is bounded by 

r 2 jg r 2 

e + p£. (1 - Or) sin 0) £ e £ e + pe (1 - Or) sin 0). 
l h a 2 n 

A useful lower bound can usually be found, as mentioned in Section 2, by 

letting e^ be the emittance of a sphere which just fits into the upper 

cavity. 

4. Application to the Solar Constant ACR 

The cavity absorptance of the solar constant ACR was assumed to be 

.999 + .001 in Appendix C . We can get a lower bound to a £ for an isothermal 

cavity using the result of the last section, case 2. 

r 2 

a e + (1-e) e , (1 - Or) sin0) 

c sphere h 

where £ , is the emittance of a sphere which fits into the cavity and can 

sphere 





be calculated using Equation 13 and 14. The appropriate dimensions are: 


r = 

. 564 cm 


h 

3.4 cm 

£ = - 166 
n 

0 

15° 


R . 
sphere 

1.128 cm 


The results for 

a few values of 

surface emittance are: 

e 

£ sphere 

lower bound to a 

c 

00 
| CO 

.9909 

.9981 

.90 

.9926 

.9986 

.95 

.9965 

.9995 

These bounds on 

the absorptance 

are significantly large 

given in some JPL reports on the 

ACR. Those estimates < 


been made using the hemispherical emittance instead of the normal 
emittance. 

We see that an uncertainty of + .1% for the cavity absorptance 
is reasonable even if significant degradation of the black coating 
occurs . 

5. Application to the FACR 

The values of the apparent emittance of the FACR given in Section VI 
were calculated from equations (13) and (14). Since the cavity design is 
more nearly spherical than the cavities used for comparison with numerical 
results, it is felt that they are accurate to better than .1%. 



APPENDIX B. ANALYSIS OF SERVO ELECTRONICS FOR SOLAR CONSTANT RADIOMETER 

This apper lix briefly describes the two temperature control servos 
for the thermal resistance difference temperature and the heat sink 
temperature. It also presents the control requirements and design and 
error details. 

1. The Difference Temperature Servo 

Referring to Figure 1, two platinum wire sensors R1 and R2 along with 

ultra stable resistors R3 and R4 measure the temperature difference AT. 

R1 and R2 are wound on the upper and lower (heat sink) end of the thermal 

resistance respectively. The bridge output is proportional to the temperature 

difference. This output is voltage amplified by the high gain amplifier. The 

amplifier output is filtered by a simple low pass R-C filter to reduce the 

noise (mostly contributed by the high gain amplifier) to an acceptable level. 

The filtered output is then fed to a zero crossing detector with small 

hysteresis. if AT is lower than the set value (set by R3, R4) , output 

of this detector controls the up /down counter to count up and vice versa. 

The 16 bit counter output along with the 16 bit digital to analog converter 

and the unity gain power amplifier provides a DC voltage that is 

controllable with very high resolution. Adding a count (count up) to 

the counter increases V, to increase cone heater power and vice versa. 

h 

At the time of every slow clock pulse a decision is made whether AT 

L 

is lower or higher than the set point and a corrective action is taken by 
counting up or down by 1 bit - slowly tweaking to regulate AT to its 
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Figure 1. ACR Difference Temperature Control Servo 
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proper value. The clock period is chosen to be long enough to allow 
temperature settling for each correction. This is done to avoid thermal 
oscillations in the servo which would easily set in because of very high 
loop gain. Note that R3, R4 and the op-amp will be placed in a temperature 
controlled environment to reduce their thermal drift errors. An auxiliary 
servo will keep these key components within +0.05°K over the measurement 
periods of 100 seconds regardless of ambient temperature variations. 

2. Difference Temperature Control Requirements 

Difference temperature setting (AT) 

Variation in AT between sun and space looks 

Variation in AT over 100 second observation 
period 

3. Design and Error Details for the Difference Servo 

Cone heater power P^ = 150 milliwatts 

Desired accuracy in P^ = 150 microwatts 

Thermal resistance of cavity support Pedestal = 500 milliwatts/°K 

The difference temperature sensors are 100 turns each of 1 mil 

diameter platinum wire wound on the 2 cm diameter pedestal. Each 

sensor resistance is 1240J2. The bridge will be excited by a stable 

6 volt source - stable within +0.05%. 

Using Vishay Technology's Ultra Precision resistors for R_ and R 

3 4 

and a MONO OP-7 as the op-amp long term drift components can be summarized 


0. 3°K + 0. 1°K 
< 2 x 10 _4o K 

< 1 x 10 _4o K 


as follows: 
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Long term drift due to R3, R4 

Long term drift due to op-amp 

Long term drift - total - due to 
R3, R4 and op-amo 

Note: The Platinum sensors may drift 

by 0.05°K over many years 
(3 years to 5 years) 

Total of above four terms 


1.39 x 10 ^°K/month 
2.22 x 10 ^°K/month 

3.61 x 10 ^°K/month 

8 x 10 ^°K to 1.3 x 10 3 °K/month 
1.522 x 10~ 3o K to 2.022 x 10 -3 °K/month 


4. The Heat Sink (Base) Servo 

The heat sink servo maintains the heat sink or the base at a temperature 
higher than the highest expected ambient temperature during observations. 

Referring to Figure 2 the heat sink heater R1 is used as a sensor-heater. 

Along with resistors R2, R3 and R4 it forms a bridge whose output voltage 
imbalance is proportional to the temperature deviation from the set point 
value (300°K) . The imbalance voltage is amplified by a high gain voltage 
amplifier and used to drive the bridge via a power amplifier. If the tempera- 
ture falls below the set point, R1 decreases, making the + input of the voltage 
amplifier more positive. This increases the bridge voltage V^; more current 
flows through R1 and R2; and R1 heats up. Rl, having a positive temperature 
coefficient, increases to the point where the bridge is balanced at a new 
equilibrium point. Due to the high gain, the bridge is very closely main- 
tained at null, and the heater winding is thus maintained at a precisely 
constant temperature. 

Note that since the heater and the sensor are the same, heater-sensor time 
delay is extremely small; hence a large loop gain is possible without oscillations. 
R2, R3, R4 and the high gain voltage amplifier will be placed in temperature 
controlled environment to minimize their thermal drift errors. The same 




Figure 2. Solar Constant ACR Heat Sink Temperature Control Servo 
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controlled environment used for the difference servo will be shared by 

/ 

these components. The environment is controlled to +0.05°K over a 
measurement period of 100 seconds regardless of the ambient temperature 
variations . 

5. Heat Sink Temperature Control Requirements 

Heat Sink temperature setting 

Heat sink temperature variation between 
sun and space looks 

Heat sink temperature, instability during 
100 second observation 

Heat sink temperature instability over a 
month 

6. Design and Error Details for the Heat Sink Servo 

The heat sink heater is made up of no. 40 Copper wire having a temperature 
coefficient of 0.005/°C. The heater resistance is taken to be 50ft and 
resistor to be 10ft. 

_3 

Bridge output sensitivity = 3.14 V/10 °K 

R2, R3, R4 with 0.001% tolerance will allow sink temperature setting tolerance 

of 3 x 0.001% x 300°K = 9 x 10 _3 °K. 

For a 6 volt +0.05% excitation, power dissipated in each sensor will 

be constant within + 36 microwatts. This is much smaller than the desired 

accuracy in P , i.e. 150 microwatts, 
h 

R3 and R4 will be Vishay Technology's Ultra Precision resistors having 


300°K + 0. 1°K 
10 _3o K 

10 _3o K/ 100 sec. 
10 3 °K/month 


following specifications: 
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Resistance tolerance 0.01% 

Temperature Coefficient of resistance - TCR 1 PPM/°C 
Self Heating 2°C/0.1watt 

Maximum power dissipation 1 watt 

Drift 5 PPM/year 

Thermocouple effect error due to 

temperature gradient across the 1 microvolt/°C 

resistor 

_3 

Bridge Output voltage sensitivity 4.5 microvolt/ 10 °K 

The Op-amp will be a MONO-OP 7 having the following specifications: 
Op-amp drift. Maximum (0.6pV/°C) 

Op-amp noise. 0.1 Hz to 10 Hz, (typical = 0.35pV P-P) 

Op-amp power supply drift coupling through power 
supply rejection ratio (PSRR = 100 DB) (power 
supply drift of 100 PPM/°C is assumed) 

Op-amp Common mode error (CMRR) = 120 DB 

Short term drift, noise error, and long term drift were calculated based 

on the above specifications. Results are summarized in the following 

table: 
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BRIDGE OUTPUT EQUIVALENT ERROR IN 

ERROR VOLTAGE HEAT SINK TEMPERATURE 


Drift due to R3 and R4 
(maintained at +0.05°K) 

0.0787yV 

2.5 

x 10“ 5 

Op-amp drift (max. 0.6yV/°C) 

0.03yV 

9.5 

x 10~ 6 

Op-amp PSRR error 

0.06yV 

1.9 

x 10~ 5 

Op-amp CMRR error (for power change 
from 1 watt to 2 watts) 

0. 322yV 

1 X 

io' 4o k 

R2 thermal EMF error (Assumes gradient = 
1/5 temperature rise) 

0. 2yV 

3.8 

x 10~ 5 

R2 self heating error (for a 1 watt 
to 2 watt power change) 


1 X 

10“ 4o K 


TOTAL drift during measurement 

cycle due to R2, R3, R4 and the , 

op-amp 2.915 x 10 °K 


TOTAL noise error will be 
essentially that due to op-amp 
noise alone; allowing small 
increase due to noise in Rl, 

R2, R3, R4 it will be about 
0.4 V PP for a 0.1 Hz to 10 Hz 

bandwidth 1.27 x 10 *°K P-P 


Long term drift due to R2, R3, R4 


3. 1 x 


10 4 °K/month 


Long term drift due to op-amp 


3.2 x 10 4 °K/month 


NOTE: Long term drift of heater wire needs to 

be investigated. (Platinum wire heater 
would give better long term stability 
but will have 25% less temperature 
sensitivity and hence drift and noise 
errors will be 25% larger - which is 
tolerable) . 



APPENDIX C. 'ERROR ANALYSIS OF SOLAR' CONSTANT RADIOMETER 


The solar constant measurement uses an active cavity radiometer 
(See Section IV) . This appendix investigates the accuracy obtainable 
from the specific ACR design illustrated in Figure 1. Some of the 
important design parameters are given in Table 1. 

The ACR operating equation is derived and studied in Section 1. 
Using this equation, the solar constant measurement uncertainty is 
calculated in Section 2. The major sources of error are also summarized 
there. 

1. ACR Operating Equation and Correction Terms 

Assuming nearly steady state operation, the power balance equations 
for viewing the sun and space are as follows 

Solar View: 

Power in = Power out + Power stored. 

. • 

A (a + pp )H + P = (P - P ln ) + P + P + CT 
c c r c e r r tr £ 


( 1 ) 
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TABLE 1 


Cavity aperture 
Area 
Radius 


1 cm 
.564 cm 


View limiting aperture 

2 

Area 1.265 cm 

Radius .634 cm 

Limiting angle 24° 

Cavity (silver) 

Area 24.6 cm^ 

Thickness .0254 cm 

Heat capacity (neglecting resistance wires) 1.53 J/°K 


Heat sink (Al) 

Volume 

Mass 

Heat Capacity 
Thermal resistance 


432 cm^ 
1.17 kg 
1044 J/°K 


Conductance 


500 mW/°K 
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where 


P 

H 

Pe 


p°ut _ p in 


tr 

P * 

CT 


- Cavity aperture area. 

- Absorptance 

- Reflectance of the cavity and cavity aperture. 

- Reflectance of upper heat sink cavity. 

- Solar constant 

- Electrical power input. 

- Radiative power exchange between the cavity and the 
heat sink and apertures. 

- Power conducted through the thermal resistance. 

- Power conducted through the electrical leads. 

- Power stored in the cavity to increase its temperature 

at a rate T. C is the cavity heat capacity. This 
term represents a deviation from steady state operation. 


Space View: 
Pe’ = (P° Ut 


• I • 

- P ) + P' + P' + CT' 

r ' tr l 


( 2 ) 


where the same notation has been used with primes indicating the space 


view. 


The form and magnitude of the terms on the RHS of equations (1) and 
(2) are discussed in the next three subsections - the radiative exchange 
terms in 1.1, the conductive terms in 1.2 and the power stored term in 1.3. 

The principle of operation of the ACR depends on using the solar 
view in conjunction with the space view. Subtracting equation (2) from 
equation (1) gives: 

A (a + pp )H + (Pe - Pe’) = [(P° Ut - P in )- (P° Ut - P in ) ] + (P - P' ) 
c c KK c r r r r tr tr 

+ (p* - V> + C < T - T ’> 

The terms on the RHS will be shown to be small in subsection 1.4. The 


uncertainty of solar constant measurements, the topic of Section 2, is b;ised 
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on equation (3) . However, the operating principle for the ACR is expressed 
by the equation: 


H . Pe ' - P e 

A (a +pp ) 
c c c 


(4) 


The cavity absorptance is assumed to be known. The error in 
can be made small because the cavity design makes a almost unity, and 


also makes a insensitive to degradations of its surface coating, 
c 


1.1 Radiative Power Exchange ((P r - P^. ) in Equation (1)) 

The power exchange between the cavity aperture and the upper heat 
sink cavity and the exchange between the cavity exterior and the lower 
heat sink cavity will be considered separately. 


1.1.1 Exchange with the Upper Heat Sink Cavity 

If the view limiting aperture' were shuttered, the power exchange 

AP U would be the same as that between two infinite planes with effective 
r 

emittances e, e c and temperatures T, T^, i.e. 

AP U = — — 9 A a (T 4 - T 4 ) 

r £+£ ~££ C C 

C C 

were A^ is the cavity aperture area. Since e and £ c are nearly equal, 

AP U = - A a (T 4 - T 4 ) 

r 2-e c c 

Now, since the view limiting aperture is not shuttered, the power received 
by the cavity is decreased by: 

A a T 4 

TT C C 

where ftap is approximately the solid angle of the view limiting aperture 
from the cavity aperture. Because the cavity emittance e is nearly one, 
the power output of the cavity is essentially unchanged. Therefore, the 



total power exchange can be approximated by: 


AP 


2-e 


A a (T 


- T 4 ) 


+ £ 


flap 


A o T 


To estimate the size of AP U let: 

r 

e =1 

A , 2 

A = 1 cm 

c 

flap/* = (^rlff ) 2 = - 055 

T = 300 °K 

AT = . 35°K, 

The result is: AP u = .21 + 2.55 mW. 

r 


1.1.2 Exchange with the Lower Heat 

• £ 

Assume that the exchange AP^ 
concentric cylinders: 


AP = 


e 2 G l 


£ 2 +e i (A 1 /A 2 ) (l-e 2 ) 1 


A. a 


Sink Cavity 

can be modeled by that for infinite 



where subscript 1 (2) refers to the inner (outer) cyclinder, i.e. 


o 4 4 

AP = £ cc A a (T - T ). 
r eff ce c 


where A is the area of the cavity exterior and e is less than the 
ce eff 


emittance e of the cavity exterior, 
ce 


Letting e e ££ 


= .3 


A = 24.6 

ce 

T = 300°K 

AT = .3°, 

o 

we get: AP =1.36 mW. 

r 


1.2 Conductive Power Losses (P + P in equation (1)): 

cr x> 

Most of the power input to the cavity flows directly to the heat 



sink through the thermal resistance. The power flow is given by: 


P = K (T - T) 
tr tr c 

where K is the conductance and - T is the temperature drop from the 
cavity to the heat sink. is about 150 mW (AT = .3°K, K = 500 mW/°K), 

There is a small conductive loss to the platinum heating and sensing 
wires, which can be written as: 


P„ = K. (T - T ) 
Z Z c cb 


where K is the lead conductance and T 


cb 


For 


is the circuit board temperature. 


= .1 mW/°K 

T - T = 1° K, 
c cb 

P z = .1 mW. 


1.3 Power Stored in the Detector Cavity (CT in Equation 1): 

The cavity deviates from steady state operation when small changes in 
its temperature occur. Although the net temperature change during a 
measurement can be made very small, this term is significant because it is 
uncompensated. That is, the temperature change occurring while viewing the 
sun is not correlated to that occurring while viewing space . 

P = CT 

where C is the cavity heat capacity. With heat sink control accurate to 

_3 

10 °K, the rate of change of temperature for a measurement time of 100 sec 

is less than 10 3 °K/sec, Then 

P = (1.53 x 10 3 ) (10~ 5 ) = .015 mW. 


1.4 Corrections to the Basic Operating Equations 

Using the results of subsections 1.1-1. 3, the solar view - space view 
difference terms of Equation (3) can be evaluated. We assume that the 



C-8 


TABLE II - CORRECTIONS 


Radiative Exchange: 


Upper heat sink cavity 

+.0061 mW 

Lower heat sink cavity 

+.0005 mW 

Conduction : 


Thermal Resistance 

+.13 mW 

Electrical Leads (circuit board 
temperature variation assumed 
. 1°K) 

+.01 mW 

• • 

Temperature Variations C(T-T') 

+.03 mW 
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temperatures of the upper heat sink cavity, the lower heat sink cavity, and 
the detector cavity are uniform. Therefore, define: 

= Upper heat sink cavity temperature 
T = Lower heat sink cavity temperature 

JL 

= Detector cavity temperature. 

The temperature differences which will appear are: 


AT 

= 

T - 

T = . 35 °K 

u 


c 

u 

AT 

= 

T - 

T = .3°K 

L 


c 

L 

6T 

= 

T - 

T ' = ±10~ 3o K 

c 


c 

c. 

6T 

- 

T - 

T t ' = ±io~ 3o k 

L 


L 

L 

ST 

- 

T - 

T ' = +10 _2o K 

u 


u 

u 

<5(AT L ) 

E 

ATj - 

- AT l ' = 2.5 x 10 _Ao K 

5 (AT) 

— 

AT - 

- AT ' = -10 _2o K 

u 


u 

u 

The 

cavity heat sink temperature difference AT is determined by the 


size of the thermal resistance. The reproducibility of AT is represented 

L 

by 6 (AT ) and is controlled by the cavity heater-sensor servo. Likewise, the 
L 

reproducibility of the heat sink temperature (ST ) is controlled by the sink 

JLi 

servo. The temperature differences for the upper cavity are assumed to be 
somewhat different from those for the lower cavity because of the different 
power exchanges experienced by the two. 

Table II gives the corrections obtained using the above temperature 
differences and Equation (5) — (9) . It is interesting that the errors caused 
by radiative exchange are essentially negligible. The only term which 
requires correction is that due to conduction through the thermal resistance. 
This term arises because the cavity servo allows a larger sink-cavity 
temperature difference during the solar view than during the space view 
(larger by 6(AT^)) . 
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2. Solar Constant Measurement Uncertainties 

The expression for the solar constant from Equation (3) is: 

H = [A (a + pp )] -1 {(Ve ' 2 - Ve ? ')/R + [ (P ° Ut - P in ) ' + (P - P’ ) 
c c r J v r r tr tr' 

+ (P £ - P £ ) + C (T - T') } 

2 

where we have used Pe = Ve /R with Ve = voltage across the cavity heater 
resistance R. The standard deviation of H is: 


a(H) = [Z (||) 2 o 2 (€.)] 1/2 

t H 1 

where are independent parameters with standard deviations 0 ( 5 .). 

Table III gives the parameters and their assumed uncertainties. For 

i 3H I 

those parameters making a significant contribution to o(H) , j— — | and 


r-r— 0 (£.) are shown in Table IV, 
o£, . 1 


The resulting value for the solar constant uncertainty is: 


4 . 

o(H) = .17 mW/em which corresponds to a per cent uncertainty of .13%. 
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TABLE III: PARAMETER VALUES AND UNCERTAINTIES 





a 

c 

.999 

+ 

.001 




A 

c 

1.000 

+ 

.0005 cm 2 




p c 

.010 


.005 




p 

.02 

+ 

.02 




R 

1000.0 

+ 

.1ft 




« 

Ve 

12,500. 

± 

2 mV 




Ve 

819. 

± 

2 mV 



1 

t =£ 
c 

.95 

± 

.05 




e ef£ 

.2 

± 

.2 




A 

ce 

24.6 

± 

c 2 

.5 cm 



ftap/ir 

.06 

± 

.01 




K tr 

500 

± 

100 mW/°K 




K £ 

.10 

± 

.03 mW/°K 




C 

1.5 

± 

.3 x 10 +3 mW- 




T 

c 

300 

± 

. 1°K 

• • » 

T, T 

(100 sec. 

measurement 

time) 

0. 

± 

10“ 5 °K/sec 

in 

H 

<1 

CT c - T l ) 

= AT 5 (T 
u c 

- T ) 
u 

.3 

± 

. 1°K 

ill 

O 

H 

(T - t') 
c c 

- ST l 5 (I L 

- V 

.3 x 10~ 3 

± 

io _3 °k 

6 T = 

u 

(T - T 1 ) 
u u 

- -«(4T ) - 
u 

-(AT - AT') 
u u 

.01 

± 

. 01°K 

6<i V 

= (4T l - 

4 V 


2.5 x 10” 4 

± 

— A 

10 °K 


5I cb 

5 T cb - 


0. 

± 

. 1°K 


2 

H ~137 mW/cm 
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TABLE IV 


9 H 


+PP ) 337 m W/cm 2 


R 


f 

Ve 


Ve 

K tr 
6(AT ) 

ij 

«(T , ) 
cb' 

T 


T' 


6 (AT ) 
u 


H 

A ~ T37 mW/cm 
c 


f P 

%+pT )H = 2 -74 mW/cm" 


c c 
P 

_ c 

C ■ c 


A c 

a c +Pp c' H = 1,37 m W/cm 2 


H __ 

R * 137 raW/fi- cm 2 
2 Ve 

A^~r = -025 mW/mV-cm 2 
2 V e _ _ _ 

A^r - *0016 mW/mV-cm 2 

S (^ )/A ( a+pp ) _ . 

L c P V - 2.5 x 10 %/crn 

K tr /A c (%+PP c ) = 500 mW/°K_ cm 2 

VVV>V , ^ mW/ . K _ ce2 

C/ VV*V - 1.5 , 10 3 n „_ sac/ . K _ 
C ^ (a ^°= ) • x io 3 *. sec/ . K . c 

O T 

•61 mW/°K-cm 2 



/ 3H | • 


I5TV__ 

(mW/cm 2 ) 


.137 


.069 

.014 

.027 

.014 

.050 

-003 

.025 
.050 
.010 
• 015 
.015 


.006 



APPENDIX D. EFFECT OF WINDOW ON SW SENSOR ANGULAR RESPONSE 


The short wave sensor requires a window to achieve the required 
spectral separation. A window can affect the angular response of the 
sensor in two ways : 

1. The reflectivity of the window varies with the angle of 
incidence with respect to the window surface. 

2. The window has an imaging effect on the transmitted 
radiation which can vary with the angle of incidence. 

The angular response for a flat window is not impaired by imaging 
effects. However, the variation of the transmission for angles of 
incidence from 0° to 65° can be as large as 30% and for large angles 
depends heavily on the polarization of the incident flux. To minimize 
the effect of the window reflectivity variations, a spherical window 
was chosen. For a spherical window with radius large compared to the 
detector size, radiation striking the detector is always incident 
nearly normal to the window surface. The imaging effects of a spherical 
window were studied with a ray tracing technique. It was assumed that the 
problem could be modelled by the two dimensional problem for which the detector 
is linear with a length 2& as illustrated in Figure 1. The response of 
the detector to radiation with incident angle y measured in a plane containing 
the detector is determined by the size, w + w + , of the incident beam which 
strikes the detector. The normalized sensor angular response can be written 
in terms of w + and w as : 

R(0) = cos0(l-E(0) ) 
w + (0) + w'(0) 
w + (0) + w (0) 


( 1 ) 
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Figure 1. Geometry of the 2 dimensional model for determining the 
imaging effects of a spherical window on the detector 
angular response. 
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where 0 is the angle of incidence and E(0) is a measure of the deviation 
of the response from the cos0 dependence of an ideal flat plate detector. 

The procedure used to find w"*" and w will be described next and is followed 
by the result for E(0) . 

First, the portion of the inside of the window irradiated by flux which 
strikes the detector is calculated. The area of the irradiated portion is 
less than that which would be irradiated at the same location if no window 
were present because the window is a diverging lens. Its focal point lies 
in front of the window and the focal length (measured from the center of 

curvature) is given by the equation, 
r(r-d) n 


f = 


( 2 ) 


d (n-1) 

where r is the window outer radius, d is the window thickness and n is the 
index of refraction. Since the detector was centered at the center of 
curvature, for normal incidence it lies in the principle plane of the window. 
For other than normal incidence, a portion of the principle plane smaller than 
the detector is irradiated. This portion is denoted by z + + z (See 
Figure 1) where z + and z are given by, 


1 „ n . fcsinO/f n 

2 - *cose [l t 


(3) 


By projecting the extreme rays striking the detector back to the focal 
+ - + + 

point from z and z the intersections (x , y ) with the inner surface of 

+ + 

the window can be determined. The results for x and y are, 

2 


x = 


£ 

f + 


2 (—) ^ 2 
(1 + V ) ~ z ] 


1/2 


(4) 


[1 + (f) ] 


y = z (l-x/'f) 


( 5 ) 


where the + superscripts have been dropped. 
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Now, from the geometry illustrated in Figure 2, w which give the 


angular response can be calculated using the following relationships: 


i = <}> + tan 1 (z/f) 


<j> = sin 1 (-^) 

- — -■ y = - (Snell's Law) 
sm i n 


( 6 ) 

(7) 

( 8 ) 


[(\t-x) 2 + (w-y) 2 ] = [r 2 + (r-d) 2 (l-2sin 2 t) ] 

2 2 2 2 2 2 ^ 

- {[r + (r-d) (1-2 sin t)] - [r -(r-d) ] } 

2 2 

w = y + [(v-x) + (w-y) ] sin (<j>-t) 


(9) 


( 10 ) 


The result for E(0) defined in Equation (1) which is valid when d/r<<l 
and l/(r-d)<<l is, 

E(0) = — y 1 - ^ - -^y- sinG. 

2n r 

It is evident that, as expected, the deviation from a purely cosine 
angular response can be made small by making the window thin and much 
larger than the detector. This relationship is evaluated in Section V for 
values appropriate to this application. 
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Y 



Figure 2. Geometry and definition of symbols used to find the width 
of the incident beam which strikes the detector. 



APPENDIX E. RADIATIVE EFFECTS OF SW SENSOR WINDOWS 


1. Introduction 

The shortwave (SW) sensor window transmits radiation in the hand 
from 0.2pm to 4ym and absorbs and emits radiation outside this band. 

Since the detector is sensitive to radiation emitted by the window, it 
is necessary to estimate the variability of this emission as the sensor 
is exposed to different long wave fluxes from the earth, sun, 
and space. Although hemispherical windows are used in the sensors, flat 
windows will be used in the thermal analysis for the sake of simplicity. 
Because of this simplification the results of the analysis will be 
approximate rather than exact. Nevertheless, they will be useful in 
identifying significant parameters affecting window emission variability. 

2 . Thermal Model of SW Windows 

Figure 1 displays the SW window configuration and the thermal model 
used in the following analysis. In this model we use the following notation: 
x = window thickness 

a = window radius 

r = radial distance from the center of the window 
T g = window sink temperature 

F e = average external flux incident on the window 

F^ = average flux incident on the window from the 

detector and detector sink 

F = average flux emitted by the window (either 

w outward or toward the detector) 




window 


x 


window 


FIGURE 



(a) ACTUAL WINDOW CONFIGURATION 



(b) WINDOW THERMAL MODEL 


1. Short wavelength sensor window configuration (a) 
and approximate thermal model (b) . 


It will be assumed that all fluxes are uniform over the window surface 


in order to calculate the radial distribution of window temperature. 

Given this distribution it is possible to calculate perturbed values 

of F which are accurate only if the perturbation is small compared to 
w 

F . The validity of this approach is obvious from the results obtained. 
3. Equilibrium Analysis 


The heat flow equation for a disc with uniform lateral power input 
is given by 

9 2 T 1 3T Cp 9T F N 

9 r 2 r 9r k 3t kx 

where 

T = T(r) = disc temperature at r 
k = thermal conductivity of the disc material 
C = specific heat of the disc 
p = density of the disc 

F n = net power flux per unit area into the disc. 


This equation assumes that the temperature gradient through the thick- 
ness of the disc is small compared to the radial gradients. 

The solution at equilibrium (8T(r)/3t = 0) is given by 
F n 3 2 

T(r) = T s + 4k^ (1 " r2/a2) 
where F n is given by 


F = F + F - 2 F . 
Ned w 


If we choose the parameter values 
x = 0.1 cm 
a = 1.0 cm 

k = 0.0138 W/cm°C (fused quartz) 



then 
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= 0.0181°K/(W/M 2 ) . (4) 

If we assume that the detector and window sinks are at the same 

temperature then F^ * F^. If we further assume that the flux incident 

from the earth is approximately given by F (r 300 W/M 2 ), for space and 

w 

earth views we have the following values of F q and F : 

F 

e 

space view 0 

earth view 300 W/M 2 

For this example the disc is isothermal at T g when it's at equilibrium. 

The steady state temperature during the space view is given by 

T (r) = T - 5 . 43°C (1 - r 2 /a 2 ). (5) 

s 

This is illustrated in Figure 2. The worst case value T(r=0) will 
be used to calculate the perturbed value F , i.e. 

270-5 4 4 

F « F (■ ) = 0.922 F . (6) 

w w 270 w 

Thus, for the assumed conditions the flux emitted by the window drops 
almost 8% between earth views and space views. This is equivalent to a 
20 W/M 2 change in background flux reaching the detector and would amount 
to a significant error if these equilibrium values were actually reached. 

4. Time Dependent Analysis 

The time dependence of the net flux input to the window can be 
approximated by a cosine function 

F.,(t) = F cos u)t (7) 

N o 

where 2 F q is the net flux difference between earth and space views and 
a) is the angular frequency of satellite rotation. It is more convenient 


IN 


-300 W/M 2 
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in the solution of (1) to rewrite (7) in the form 

F (t) - ± F e lmt + | F e - lmt ' 

N 2 o 2 o 

and to solve (1) for each term. Substitution of the trial solution 

T = T + IC(r) e ±1COt 
o 

into equation (1) yields the R equation 


( 8 ) 


(9) 


1 R - . 


-r— r R + 4- ia)a R = £ 

3 r 2 r r 

where we have used the notation 

F 

“ - Cp/k, E - . 

If we substitute 

7 - ( i • ^ 1/2 

Z = (+iwa) r 

and multiply (10) by r 2 we obtain an inhomogenous form of Bessel's 
equation 


( 10 ) 


(ID 


( 12 ) 


,2 3 ' 


9 „± 


z^ —y r~ + z -~r~ + z 2 r- •- e - . 

“i 4 oi + loaa 

which 'has the solution 

R(r) - = 5 - 7 — + A J (r (+iwa) ^^) . 

+ iu>a o 

Since there can be no time dependence at r = a (T = T g ) , A must have 


E Z" 


(13) 


(14) 


the value 
A = 


1/2 

±ima J (a(qFi^a) ) ’ 
o 


(15) 


Thus the final pair of R solutions is 


1/2 

+ E J Q (r(+ima) ) 

R(r) ” +iaia, [1 " J (aGPiu^l/Z)!- 


(16) 
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The values at r = o, the point of maximum temperature deviation, are 
given by 


RCo) 1 ==r^ — [1 - 1/J (a(+ima) 1/2 )]. 

+iwa o 

In order to evaluate the Bessel function it is useful to write 

J (a(-iu)a) ^ 2 ) = J (a\/<jja e 377 ^ 1 *) 
o o 

J (a(+itoa) 1/2 )= J (aVwa" i 37rl/l <) 
o o 

from which we can obtain expressions in terms of the real functions 

ber and bei , i.e. 
o o 

J Q (a\/a)a e 311 ^ 4 ) = ber^ ( 2 \/ toa) + i bei^ (a\Aoa) , 


J Q (a\/ukx" e 311 ^^^ 1 *) = ber Q (e\/u;a) - i bei Q (a\/aju) . 
Defining 

x = ber (aVuaj 
o 

y = bei Q (a\/oja) 


the time dependent solution r = o can be written as 

„ icut -imt „ imt -iwt 

s _ , E r e . e 1 E r e . e 

T(o,t) = T + — [ — r + 77 ] - — [ — ■ ■. — N + —r~, : — s 

o ma -1 +1 u>a -i(x+iy) +i(x-iy) 


] 


which can be reduced to 

9 T? /" X 1 / 2 

T(o,t) = T o + — V2 (1 - x 2 + y 2 ) cos(ut+4») 

where 

* - tan’ 1 + 1% . 

y 

For very low frequencies, i.e. for a Vu>a << 1 we find that 
x = 1-64 <^) 4 + ... 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


( 27 ) 
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, , ,aVwa N 2 , 

y * 16 ( — g — ) + 


(28) 


Thus as u ■> o 


x 1/2 sfi o 

(1 -^> + f ^ 


x ~ (x 2 +y 2 ) 


o. 


(29) 

(30) 


Inserting these limiting forms into (26) and (27) yields the low 
frequency expression 

g a 2 

T(o,t) = T o + — — cos rnt, m -> o. 

Comparison of (31) with (2) shows that, as expected, for very low 
frequencies the temperature variations are equal to the equilibrium 
values and no phase shift is present. 

However, for the expected operating conditions 
(u * 0.628 sec ^ ( 6 RPM) 

C = .753 W-sec/gm.°C 

3 

p = 2.203 gm/cm 
k = 0.0138 W/cm. °C 
a = 1 cm 
we f ind that 


(31) 


a\Aua = a \)4KJp/k = 8.689 
x = 39.44 
y = -50.05. 

In this case 


, x 1/2 


< l - -- 

(32) 

(|) = . 496tt ~ y 

(33) 
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and we can write (25) to good approximation as 

VTf 

T(0, t) = T + — — - cos (wt + tt/ 2) . 
o (x)C* p x 

o 

For F = 150 U/M and x = 0.1 cm (3) becomes 
o 

T(0, t) = T + 0.204°K cos (ait + tt/ 2). 
o 

If we assume that (35) applies to the average window temperature, a 
worst case assumption, then the background flux variation induced is 
given by 

4 V2*F 3 

F B (t) = ° T s COS (wt + * /2) 

which, for = 273°K, takes on the specific form 

F D (t) = 0.941 W/M 2 cos (uit + ti/2) . 

D 

This represents a peak to peak variation of 0.2% of the solar flux 
and 0.6% of the average outgoing shortwave flux from the earth. How- 
every, in actual operation these errors will have very small effects 
because of the tt/2 phase shift between window absorption and window 
temperature. Since data samples will be taken near points of maximum 
and minimum incident long wave flux inputs (i.e. at earth nadir and anti 
nadir), the window temperature deviation will be near minimum. This is 
illustrated in Figure 3 for a cosine variation of input flux. 

5 . Conclusions 

As a result of varying long wave flux inputs to the SW window 
its temperature varies sufficiently to produce variations in background 
flux levels incident on the detector. Significant features of this 


(34) 


(35) 


(36) 


(37) 


variation are as follows: 
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Figure 3. Relationship between window radiative heating and window 
temperature variation. Note the tt / 2 phase shift. 
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(1) the magnitude is independent of the window diameter. 

(2) the magnitude is inversely related to window thickness 
and rotational frequency. 

(3) for probable design parameters there is a phase shift 
nearly equal to tt/2 between incident long wave flux 
variations and background flux variations. 

(4) for probable design parameters background flux 
deviations should be negligible. 



APPENDIX F. FLUX ERRORS PRODUCED BY EMISSIVITY 
VARIATIONS WITH WAVELENGTH 


1. Basic Equations 


If a sensor with emissivity e(v) at wavenumber v is calibrated at 
temperature T , then the signal output in flux equivalent terms is given 


by 


n. 


(v)F(v) dv 


me as 


c(v)B(v,T c )dv 


crT 


( 1 ) 


where 


F(v) = Spectral incident flux 


B(v,T c ) = Spectral flux produced by a blackbody at temperature T c 
4 

oT^ = Total flux produced by a blackbody at temperature T^ . 

The ratio between the measured flux F and the exact total flux 

me as 


F = C F(v)dv 


( 2 ) 


is given by 


F e r 

meas F 


( 3 ) 


where 


. 

7 C = t /e(v)B(v, T c )dv]/aT c 4 
o 

oo oo 

Ej, = [ ^/*e(v)F(v)dv]/tJ*F(v)dv]. 


( 4 ) 

( 5 ) 


If F(v) has the same spectral dependence as B(v,T c ) , or if e(v) is independent 



F-2 


of wavelength then and the measured flux is the same as the true 
flux. However, since neither of these conditions are precisely net, the 
ratio expressed by equation (3) will differ from unity. 

2. Estimation of Errors 


Given a description of the surface emittance properties e(v) the flux 

error produced by this variation can be estimated by calculating e for 

F k 

each k, of a set of n spectral flux distributions F (v) 

k 

which represent typical cases for the earth. 

til 

The fractional error. in the k case is 
F -F £ F 

r k meas,k , k ... 

6 1 - — ( 6 ) 

k e 

c 

The bias error produced by the specific surface for the n cases is just 


— 1 n 1 1 n — 

6 = - Z 6=1-—- Z e- ] . 

" k-i k 7 " k=l F k 

c 

The standard deviation about the mean is given by 
n 1/2 

° ' [ 7 kli (4 k - 4 > 2 ! • 

Substituting (6) and (7) into (8) also yields the form 


(7) 


( 8 ) 


1 n i n — — 2 

a ' k=l [ 7 j = l £ j - £ k ! 1 

ne j j 


1/2 


(9) 


which can also be written as 

1 1 n — 2 1 n — 2 

0 ’ ^ '7 - [ 7 k£i E k J 1 


( 10 ) 


3. Test Conditions 


Three surface emissivity distributions were included in this study 


these are listed and described below: 
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(a) Parson's Black Lacquer (e^ at 180°C) ; data from Sydnor (1970) 

(b) 3M Black - 101-C10, baked (e^ at 77°K) ; data from Hall (1970) 

(c) 3M Black - 101-C10, baked (e v at 373°K) ; data from Hall (1970) 

Three distinct spectral distributions (Wark, 1962) were used in the study; 
their identifications are in terms of atmospheric conditions: 

(1) Cloudless Winter Arctic 

(2) Cloudless Summer Western US 

(3) Clouds - covered Tropics 

Figure 1 displays- the surface emissivities as a function of 
wavenumber and Figure 2 displays the spectral flux distributions 
for the three atmospheric conditions. 

4. Results 

The variation of with calibration temperature T^ is presented in 
Table 1. 

TABLE 1 . Average Emissivities for Blackbody Spectral Fluxes as a 
Function of Blackbody Temperature T^ . 

£ (a) e (b) £ (c) 

c c c 

0.915 .896 .935 

0.915 .898 .935 

0.915 .901 .934 

0.916 .903 .934 

0.916 .905 .934 

0.916 .906 .933 

0.917 .908 .933 

0.917 .909 .933 

0.917 .911 .932 

0.918 .912 


T c (°K) 

200 

210 

220 

230 

240 

250 

260 

270 

280 


290 


.932 
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The values of e are presented in Table 2 for each surface. Also 

K. 

included in this table are o and 5 values assuming a calibration temperature 
of 270°K. 

TABLE . 2 . Average Emissivities for Spectral Fluxes Typical of Earth 
Emissions . 

SPECTRAL FLUX CLASSIFICATION £ k^ b) E k ( ' C ' ) 


(1) 

Winter Arct. 


.918 

.903 

.933 


(2) 

West. US Summer 


.922 

.914 

.926 


(3) 

Cloudy Tropics 


.919 

.912 

.929 



1 3 7 

3 k=l £ k 

= 

0.920 

0.906 

0.929 



e 

c 

= 

0.917 

9.909 

0.933 



a 

= 

1.8 x lO" 3 

6.6 x 10~ 3 

3.1 x 

lO' 3 


1 

= 

-3.3 x 10' 3 

-0.7 x 10~ 3 

+3.9 x 

lO' 3 


lu all cases there is a bias error produced using a blackbody calibration 
procedure. The case of Parson's black demonstrates this most clearly. It's 
effective emissivity for absorbing earth emitted radiation is larger than its 
blackbody emissivity because the spectral distribution of earth fluxes always 
contain minima produced by water vapor and CO 2 absorption which serve to reduce 
the effect of the low emissivity of Parson's black at 660 cm The spectral 

emissivity of 3M black at 373°K shows an increase in this region and thus 
produces the opposite effect. In both cases bias errors of 0.3% - 0.4% 
are possible. 

The random errors varied from 0.2% -'0.7%. Although these might seem 
quite tolerable, it should be noted that these errors are actually correlated 
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with latitude, season, cloud cover, etc., as a result of the spectral 
flux distribution characteristics. Thus these "random" errors can result 
in bias errors in north-south gradients and other significant parameters 
as well. 








RADIANCE 
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WAVE NUMBER (CM**-1) 


FIGURE 2 : 


Radiance vs. Wavenumber 



APPENDIX G. FLUX ERRORS RESULTING FROM INCOMPLETE SPECTRAL COVERAGE AND 
IMPERFECT SEPARATION OF SOLAR AND LONG WAVE COMPONENTS 

Although It is possible to obtain spectral coverage adequate to measure 
total radiation leaving the earth, the radiation incident from the" sun, 
and thus the net radiative energy exchange, all experiments which attempt 
to separate reflected solar radiation from earth-emitted radiation must 
face a common problem: spectral separation techniques are not rigorous. 

The total radiation emitted by the earth can be rigorously separated from 
the total solar radiation reflected by the earth only if they lie in entirely 
separate spectral regions. On the planet earth this is not quite the case. 

In the spectral region from 4y to 5y there is a small, but not entirely 
negligible, overlap between the two components. Even a perfect long wave 
cutoff filter cannot make a complete separation. A poorly designed filter 
might produce significant errors. It is the purpose of this appendix to 
estimate how well this separation can be made in theory and how significantly 
the filter (window) cutoff wavelength affects the separation error. 

1. Approximating the Missing Flux 

The spectral flux from a blackbody at temperature T is given by 
. c v/T -1 

B v (T) = C^te - 1] (1) 

where T is the absolute temperature (°K) v is the wavenumber in cm ^ and 
the two radiation constants are 

= 1.1906 1 x 10 erg/cm^sec ^ 


C 2 = 1.43868 cm°K. 


( 2 ) 
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2 _ i 2 —i 

The units of B^(T) for these constants are erg/(cm -sec-(cm )) or mW/m -cm . 
Approximations for B which apply at large distances from its peak value are 


„ r 3'V /T 

B v» ,C l v e C 2 v/T » 1 


1 2 

B vS I(^)v 


C 2 v/T « 1. 


(3) 

(4) 


If we assume that the earth is a blackbody at temperature T^ then the total 


flux emitted over all wavelengths is given by 


CO 

1 

J O 


F(T e ) = | Bv(T E )dv. 


(5) 


If our measured flux is wavelength limited, i.e. it only includes wave- 
numbers between a lower cutoff v^ and an upper cutoff then the measured 

flux is given by 


V 2 °° 1 “ 

VV ' / B'XVdv - jB v (I E )dv - lB v (T E )dv -jB v (T E )dv (6) 

■v, J . v- 


If C 2 V 2^ >> an< * ^2 V 1 <<C ^ t ^ ien we can use t ^ ie approximations given by 
equations (3) and (4) to evaluate equation (6). The integral terms are 


given by 


Jf 


VV dv 5 f 


(7) 


V 2 

t 4 C 2 V 2 _1 _C 2 V 2 /T E 

/ B v (T E )dv » 7T Cl v 2 i — ) e 

E 


( 8 ) 


We can approximate equation (6) by the expression 

, . C 9 v. -1 , -C v 9 /T C v -1 

W ■ F( V - ” C 1 <■-¥■-> + ^ <4r> >• 

E E 


( 9 ) 
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-1 


If we take as an example v = 0 and v„ - 2500 cm (4y) and 


T„ = 273°K we find that 
E 


C 2 V 1 C 2 V 2 

44 = 0, 44 = 13.17 


( 10 ) 


and 


F (T ) - F(T ) = 0.212 W/M 
m E E 


( 11 ) 


We can calculate F(T„) using the relation 


/■ 

•'O 


^1 4 4 

VV dv =15 7^ T E = ffT E 


( 12 ) 


_5 2 —8 2 
where a = 5.669 x 10 mW/M = 5.669 x 10 W/M . 


Employing equation (12) , 


we find that 


F(T e ) = 314.9 W/M , and 


(13) 


VV- F( V . .. , n -4 

“ 6 - 73 x 10 


(14) 


- 1 , 


Thus a 4yi lower wavelength cutoff (v^ = 2500 cm ) leads to a missing flux 
which is 0.07% of F(T E ). 

If we use a 4p upper wavelength cutoff on reflected solar radiation 
(T s ~ 5900°K) then we find 


C v C v 

~ ~ ~ = .6096, 44 - = ». 


(15) 


Although C„v /T does not satisfy condition (3) the errors introduced by 
using this approximation (Equation (3)) are relatively small and positive, 
i.e., , the actual error will be slightly smaller than the value calculated 
using equation (10). 
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The fraction of missing flux of reflected (or incident) solar radiation 
is estimated to be 


C v 3 

-5 ( ■ ) 

F (T ) -F(T > v T * 

m s s _ s 

F(T ) ' k ~ 

S TT 


= -.0116. 


(16) 


Thus we find that approximately 1.2% of the solar flux is excluded using 
a 4y cutoff. In order to keep the omitted flux fraction below 0.1% we 
must have satisfy 


- ^2 V 1 ^ 4 

(-^H4 < .001 x ir /5 = 0.019482 

s 


which implies that 
v 1 < 1103 cm" 1 


(1 > 9. 06p) 


(17) 


(18) 


However, this would lead to the inclusion of large percentages of long 
wave flux from the earth. 

If we attempt to optimize the cutoff wavelength to produce equal per- 
centages of omitted flux in both spectral regions we find that must 
satisfy the condition. 


V i,V 

v l C, tn 3 ' ' 

2 s 

The dependence of v on T and the percentage of missing flux for each 
case are displayed in Table 1. 


(19) 


2. Estimating Cross Band Contamination 

No matter which cutoff wavelength is selected from Table 2, there will 
always be some observing conditions which result in missing fluxes of the 
order of 1%. In addition to the problem of not measuring part of what 
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should be measured, the spectral overlap also results in both SW and LW 
measurements containing flux contributions that they shouldn't be measuring. 

The SW band will accept some of the earth emitted radiation, and the LW 
band will accept some of the reflected solar radiation. 

Two examples of cross band contamination are presented in Table 2. 

Both cases are for a 4.57p cutoff wavelength. Case I approximates the 
radiation levels observed from a high albedo, high altitude cloud (cloud 
temperature 190°). Case II approximates the radiation levels present 
above a warm dark ocean (surface temperature 3 ; 300°K, albedo ~ 0.02). As 

Table 1. Values of and the corresponding A^ which result in omission of 

the same percent of flux for both long wave and short wave measurements.* 


h 

v(cm 1 ) 

Hv 0 

% MISSING FLUX** 

200 

1564.2 

6.39 

0.28% 

220 

1677.9 

5.96 

0.35% 

240 

1785.8 

5.60 

0.42% 

260 

1891.2 

5.29 

0.50% 

273 

1958.0 

5.11 

0.56% 

280 

1993.4 

5.02 

0.59% 

290 

2043.4 

4.89 

0.63% 

300 

2092.6 

4.78 

0 . 68 % 

320 

2189.1 

4.57 

0.78% 


* The SW measurement is assumed to include the range (0 -A^) anc * 

the LW measurement 0 - (A^ - °°) . 

** The percentages are defined as reflected outside the SW passband 

* total reflected solar , or^emitted radiation outside the LW. passband 


-j- the total emitted. 
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Table 2. Two Examples of Cross Band Contamination for a Wavelength Cutoff 


of 4.57 : 
Fluxes are 

Case I = 
in units 

bright cold cloud; Case II 
of W/M -2 . 

= dark warm ocean* 


CASE I 

CASE II 


CASE I 

CASE II 

TOTAL 

SOLAR 

700.0 

28.0 

TOTAL 

EMITTED 

73.88 

459.2 

SOLAR IN 
SW BAND 

694.5 

27.8 

EMITTED 
IN LW BAND 

73.88 

456.9 

EMITTED 
IN SW BAND 

.003 

2.26 

SOLAR IN 
LW BAND 

5.46 

0.2 

TOTAL FLUX 
IN SW BAND 

694.5 

30.1 

TOTAL FLUX 
IN LW BAND 

79.34 

457.1 

% DIFF BETW. 
SW IN BAND 
AND TOTAL 
SOLAR 

-0.8% 

+7.3% 

% DIFF BETW. 
LW IN BAND 
AND TOTAL 
EMITTED 

+7.4% 

-0.5% 


* The solar fluxes are based on wavelength independent reflectivity. 


indicated in the table contamination levels can be considerable both in 
absolute terms and as percentages. It should be noted that changing the 
cutoff wavelength will change the distribution of the contamination but 
will do little to reduce the contamination level. 


3. The Great White Hope 

The futility of attempting to achieve accurate spectral separation 
by means of filtering alone suggests a simpler method. Consider, for 
example, a radiometer with a SW passband from 0.2y (or O.Oy) to » where 
A 2 is short enough to exclude all significant earth emission. If this 
radiometer is calibrated against the sun so that it reads H g even though 
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Figure 1. Comparison of the near- inf rared solar spectrum with laboratory 
spectra of various atmospheric gases. (AFCRL, 1965) 



it misses all the solar radiation from to °°» anc * if the earth consists 
of a black (or dark grey) ball with perfectly white (or light grey) patches 
gliding over the surface (and no atmosphere), then, when the f orementioned 
radiometer views the earth, it will produce a reading which is exactly 
equal to the reflected solar radiation even though it does not accept all 
the reflected radiation from X ^ to In other words, the calibration pro- 
cess forces the radiometer to approximate the missing flux as if the average 
albedo from X^ to <® were the same as the average albedo from 0 to X^. 

For the idealized earth described above (the patches were meant to repre- 
sent clouds) the approximation works perfectly. For the real earth, with 
clouds made of water droplets and ice particles floating in an atmosphere 
sprinkled with water vapor and CO 2 , the approximation is highly questionable. 

As indicated in Figure 1 , there is significant atmospheric absorption 
in the overlap region from 4y to about 8y, mainly due to water vapor, CO^ 
and N 2 O. Thus solar radiation reflected at the bottom of the atmosphere 

leaves the top of the atmosphere with considerable attenuation in the 4y to 
8y interval violating the ideal of the white (grey) hope idea. The spectral 
reflectivity of clouds, which account for most of the reflected solar radiation, 
is even more important. Estimates of cloud reflectivity (see Figure 2, for 
example) indicate a very sharp drop in albedo beyond 3y due to the large IR 
absorbtion of liquid water drops in this region. It thus appears that reflected 
solar radiation in the overlap region is likely to be much smaller than would 
be deduced on the basis of grey (wavelength independent) reflectance. 

4. Choosing a Bandpass 

Based on the previous analysis the most reasonable way to separate 
reflected solar and emitted long wave components is to choose a cutoff 
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between 3y and 4y. This will insure that all significant emitted 
radiation is accepted by the long wave passband. Since most of the 
incident solar energy in the LW passband will not be reflected, the 
LW band will suffer little contamination and the SW passband will 
not miss any significant amount of reflected flux. This will require 
a correction factor, however, for using the solar calibration since 
approximately 1.5% of the solar flux will not be transmitted by the 
window. 

It should be noted that Figure 2, if extended to 6y would probably 
show some regions of increased albedo above 0.2. If the mean value is 
below 0.2 from 3.5y to about 6p the missing SW and LW contamination 
errors should still be within tolerable limits. 



\ (microns) 


Figure 2. The spectral albedo for middle layer clouds 
(after Novoseltsev, 1964). 



APPENDIX H. FLUX ERRORS RESULTING FROM NON- LAMBERTIAN 

ANGULAR RESPONSE 


1. Modeling Sensor Angular Response 


The ideal response of a plane flux sensor is that of a perfectly 

absorbing flat surface; i.e., the absorbed energy from a point source 

at an angle 0 away from the plane normal is proportional to cos0, the 

projected area of the surface normal to the source direction. 

We shall consider only one model for deviations from the ideal 

2 

cos0 response. This model is a linear combination of cos0 and cos 0, i. 

2 

R(0) = (1 -a) cos0 + acos 0, 0 £ a<l, 

where R(0) is the sensor relative angular response and a is a deviation 

parameter which can be varied between zero and one. 

The relative angular emissivity variation corresponding to equation 

(1) is just e(0) = e o R(0)/cos0, i.e. 

e(0) = e [1 -a (1 - cos0) ] . 
o 

As illustrated in Figure 2, the emissivity thus varies between at 0 = 

and e (1 - a) at 0 = 90°. 
o 

2. Modeling Angular Distribution of Radiation from the Earth 


Since our sensor will be in orbit around the earth the angular 
diameter of the earth will be tt/y, where Y > 1 and depends on orbital 
height according to the equation 


9 R 2 1/2 -1 

Y = [1 - 7 sin" 1 [1 -(^+h) ] ] . 


where R g = 6370 km (radius of the earth) and h is the satellite altitude 
For the value of h = 380 n.mi. (704 km) we find that y = 1.4014. 
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This corresponds to an earth angular extent of -64.22° $ 0 s + 64.22°. 

We shall consider three different models for earth radiance angular 
distribution 1(0), corresponding to uniform brightness, limb darkening, 
and limb brightening. These are defined as follows: 


IjO) = I ± Q , -n/2 < y0 $ tt/2 

I 2 (0) = I 2 q cos 2 y0, —tt/2 $ y0 $ tt/2 

I 3 (0) = I 3 0 sin 2 y0, -tt/2 S Y 0 $ tt/2. 

The constants I. n , I„ n , I- n are chosen to satisfy the normalization 

■1 j w 4 j U J) U 

condition 


Tf/2y 

F = 2tt J 1(0) cos0sin0d0 


where F is the plane flux incident on the detector. In order for all 
angular distributions to yield the same plane flux F (equation (7)) we 
find that the coefficients must satisfy 

F r . 2 tt .-1 

1 1. 0 = 7 [sln 2Y 1 

T _ J r . 2 tt 1 2 it .-1 

1 2.0 2 ir [sln 2Y + , .2 COS 2y ] 


t « F r • 2 tt 

I 3, 0 = 2 * [sin 2Y 


1-Y 
1 

1-Y 


2 TT .-1 
2 cos — ] . 


(4) 

(5) 

( 6 ) 


(7) 


( 8 ) 

(9) 

( 10 ) 


3. Sensor Calibration 


We shall assume that the sensor is calibrated with isotropic radiation 

limited to +tt/ 2Y in 0. The incident flux in this configuration is given by 
tt / 2 Y 

F. = I / 2irsin0cos0d0 = tt I sin 2 , 
icy c 2 y 

0 

where I is the angle independent radiance used in the calibration. The 


(ID 


absorbed flux is given by 
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tt/2y 


F = I / £ [(1-a) + acosG] 2rsin0cos0d0 

a c J o 


which reduces to the form 


( 12 ) 


F a = trl^ e Q [(l-a)sin 2 + -|a (1-cos 3 ) ] . (13) 

The calibration factor K which is subsequently used to derive incident 
flux from absorbed flux is thus given by 


F 2 

K = (_A) = sin (tt/2y ) 

F a cal £ q [ ( l-a)sin 2 (fr/2Y) + -|a (1-cos 3 -|— ) ] 

where the subscript "cal" on the flux ratio denotes the ratio existing 

for the calibration condition. For angular intensity distributions other 

than isotropic, K will not be equal to the ratio of incident and absorbed 

fluxes, although, in practice, sensor output will be treated as if it were 

true. In other words, we shall define the "measured" flux as 

tt/2y 

F = K F = K e / I (0) R(0) 2itsin0d0, 
meas a oj 

. o 

where in this case F and 1(0) denote absorbed fluxes and correspondent 

a 

incident intensities which might be encountered in orbit. 


(14) 


(15) 


4. Estimating Flux Errors 

Measured fluxes, i.e. those calculated according to equation (15), 
will be determined for 1^(0) and 1^(0) • The difference between corresponding 
measured fluxes F£ and F^ and the incident flux F are then the measurement 
errors. (Since 1^(0) is independent of 0, the calibration procedure assures 
zero error for this case) . 

The equations for F^ and F^ are: 
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Tr/2y 


F = K e I 
2 o 2 


/ 2 2 

cos Y0 [(l-a)cos© + acos 0]2Trsin©dO, 


tt/2y 


F 3 K £ o ^3,0 sin^Y© [(l-a)cos© + acos 2 0] 2irsin0d0. 


If we define two define integrals 


j/2y 2 


A = jf* cos Y© cos0sin0dy 


tt/2y 


/ 2-2- 


B = J cos Y© cos 0sin0d0 , 
o 


then Fp and F^ can be written as 


F„ = 2itK e I 0 [ (1-a) A + aB] 
2 o 2,0 


F„ = 2ttK e I. [(1-a) sin 2 ^ A) + aB] 

3 o 3,U 2 2y 


where 


A _ 1 / . 2,-rr x , 2 ,tt x 1 x 

A = 4 < sin (27) + cos < 27 >- 2 } 


‘Y !_• 


Y 


B - i COS 2 (J^> + f ( 1-cos J (^)) + * cos‘(|i) ? 


3 2,3n* 1 

4 cos <47' 


Since we can write the intensity constants in the form 

I = -* A' 1 

1 2.0 2 ir A 

T _ F ,1 . 2 TT -1 

1 3.0 2tt ''I 8111 2y^ “ A) 

equations (20) and (21) can be rewritten in the form 


F 0 = K e F [(1-a) + a 
2 o A 


B 


1 . 2 n 
"2 s in (jy) ~ A 


]• 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


F_ = K e F [(1-a) + a 
3 o 


( 27 ) 
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Evaluating A, B, and K for y = 1.4014, we find 


B = + . 147848 
A .15367 


.96212 


B 


= +.147848 

1 . 2.tt . . .25177 

2 Sin " A 


= .58723 


£ q K = (1 - 0.3882a) 1 . 


Thus the fluxes take on the specific values 

1-. 03788a 
2 1-. 24546a 


F„ = F 


1-. 41277a 


3 1-. 24546a * 


As a function of a we find the foilwing percentage errors in and F^ 

(LIMB DARKENING) (LIMB BRIGHTENING) 



Fo-F 

F -F 


100 x 

100 x 

a 

F 

F 

.01 

+ .2% 

- .2% 

.02 

+ .4% 

- .3% 

.05 

+ 1.1% 

- .8% 

.10 

+ 2.1% 

- 1.7% 

.20 

+ 4.4 % 

- 3.5% 

.40 

+ 9.2% 

- 7.4% 

.60 

+14.6% 

-11.8% 

.80 

+20.7% 

-16.7% 

1.00 

+27.5 % 

-22.2% 


Since typical limb darkening and brightening is much less than 100% the 
errors just tabulated should be considered as very conservative. 


(28) 

(29) 

(30) 

(31) 

(32) 



APPENDIX I. ESTIMATION OF FACR ANGULAR RESPONSE 
ERRORS DUE TO NONUNIFORM RADIATIVE COUPLING 
BETWEEN CAVITY AND SINK 

For the parameterization of errors due to nonuniform radiative 
coupling between the cavity and its surrounding heat sink we will model 
the cavity as a linear plate which is divided into two regions of 
distinctly different coupling. As indicated in Figure 1, the plate is 
divided at the midpoint into Region 1 and Region 2. Significant 
parameters required to describe the model are: 
s = the half length of the plate 
X = the thickness of the plate 
k = the thermal conductivity of the plate 
z = the distance along the plate 

e = the emissivity of the top surface of the plate 

^1*^2 = t * ie ra diative power transfer between the plate and the sink 
per unit area per °K for regions 1 and 2 respectively. 

T^(z),T 2 (z) = temperature distributions in region 1 and 2 respectively 

F^, F 2 = incident radiative fluxes in region 1 and 2 respectively 

(these do not include the radiation received from the sink) 

The transfer coefficients k^ and k 2 are mathematically defined as follows: 

k 1 = 4 E;L aT 3 (1) 

k 2 = 4e 2 aT 3 (2) 

where T is the mean temperature of the entire plate and e^, e 2 are effective 

emissivities for radiation transfer between plate and sink. The solution 

for the temperature distribution along the plate is determined to be 

„ eF 

m/\ . r a i z — 2a is -a,z, , 1 /*^\ 

T-^z) = A 1 [e 1 + e 1 e 1 ] + — — , (3) 
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m / \ » r 3.9^ i 4"2a 2 s a 2 z, 

T 0 (z) = A [e z + e z e * ] + 


eF 

2 

V ’ 
k 2 


where the constants are given by 

F„ F a 

A.. = e(-: — - - — )(1 - e aiS ) -1 [coth a,a H coth a s] -1 

1 k 2 ^1 1 2 

F F 1 a a 

A 2 = j^-)— (1 “ e a2S )" 1 [ coth a 2 s + -~c°t:h a 2 s] _ 

a l ■ 4 )V2 • a " d a 2 ‘ • 

The mean cavity temperature T is found by integrating T^(z) and T 
the length of the "cavity" and dividing by 2s. The result is 


T = 


a l 2 F 2 F 1 

(1 ^ 
e a 2 R 2 1 

2a 2 s 2 coth a,s coth a„s 
1 [ — + 1 


a l S 


a 2 S 


If we define the temperature servo origin so that T 
(8) implies a fixed relationship between F-^ and F 2> 

F 2 = -F 1 (l + -|hk 2 )(l +|hk 1 )" 1 

where the parameter h is given by 

s 2 

h = -j^[ f (a^s) + f(a 2 s)], where 
f(0 = coth Z ~ . 

For £ << 1 (the usual case) we can make the approximation 

fU) = 3 , 


is fixed at 0 

i.e. 


in which case (9) can be simplified to the form 


, k 7 - k 

F 2 - - F lO+ 3 ( -k^ )s 1 


(4) 

,(5) 

1 > ( 6 ) 

(7) 

2 ( z) over 

( 8 ) 

° , then 

(9) 

( 10 ) 
(ID 

( 12 ) 


( 13 ) 
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Defining the radiative flux inputs to regions 1 and 2 as and r> and 

F^ as the electric power input per unit area (which^ is uniform over the 
cavity surface) we find that 


F = -F + F, , and 
1 e e l,r 


F 0 = — F + F 
2 e e 2, r 


(14) 

(15) 


Inserting (14) and (15) into equation (13) yields the following expression 
for the electrical power flux in terms of the incident radiative fluxes: 


F 

e 


- [ f (F l,r + F 2,r> + ! 


(k - k ) 

— T ^-s 2 F, 1 

kx l,r 


(16) 


If we consider two extreme cases with different incident flux distributions, 
i.e. 

case (1): F, = F ; F„ =0 

l,r o 2,r 

case (2): F- = 0 ; F„ = F 

l,r 2,r o 

then the peak to peak fractional difference in electrical power flux for 
the same total radiative input power is just 


f < 2) - F g 1) x (k - k ) 
e e 12 1? 

F " 6 e k^ S 


(17) 


If the mean cavity temperature has the absolute value T, then the error 
expression can also be written as 


f ( 2) _ F (l) , _ . 

6 _ - 6 - \ COT 3— t- 

F„ 3 kx 


(18) 


where is the difference in effective emissivities for radiation 


transfer between cavity and sink. These emissivities depend on the 
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cavity-sink geometry as well as the surface emissivities of the two 
components . 

The s 2 factor in equation (18) applies rigorously only for the linear 
model used in the derivation. There is considerable uncertainty in choosing 
an appropriate geometry factor for a three dimensional cavity which will 
probably be resolved only by measurement. For order of magnitude estimates 
we will choose s to satisfy 

4s 2 = A (19) 

s 

where A g is the surface area of the cavity, in which case s % /tt x the 
cavity radius. Inserting (19) into (18) yields the final result: 


F (2) - F (1) eoT 3 A 

— = ^ — = ( e - e ) 

F 6kx ^2 V* 


( 20 ) 


For typical parameter values (see Section VII) of e = 1, A g ~ 19 cm , 
X a 7 x 10 -3 cm, T = 350°K and k = 0.58 W/cm°K, equation (20) has the 
value 


eoT' 


6k X s 


A ~ 0.19 


( 21 ) 


For a 3% difference between and z^ this would imply a fraction power 
error of 0.5% between the two regions. Based on the results of 
Appendix H, this should result in angular integration errors not exceeding 
0 . 1 %. 



APPENDIX J. PARAMETERIZATION OF THE FACR ANGULAR RESPONSE ERRORS RESULTING 
FROM LEAD WIRE AND SUPPORT CONDUCTION 

As indicated in Section VI. 5, lead wire and support conduction can 
produce angular response errors. In general, radiative flux incident on 
the cavity near a support or lead wire contact point has a somewhat 
smaller cavity heating effect than the same amount of radiation incident 
on other parts of the cavity. In order to make an order of magnitude 
estimate of the resulting errors, and to determine how cavity parameters 
can be optimized to reduce this effect, a similar, but geometrically 
simpler configuration is treated. The cavity is replaced by a disc with 
a single support wire. The electrical servo power is calculated for 
radiative power incident in an annulus outside this region. The power 
difference for these two cases is finally related to an angular response 
error. 

The disc model for treating the effects of wire conduction is illustrated 
in Figure 1. Parameters used in this model are 
a = Radius of the wire contacting the disk 
£ = The wire length between cavity and the sink 

b = The radius separating region I and region II on the "cavity" 

c = The outer radius of the cavity 

T g = The sink temperature 

Tj(r) = The cavity temperature distribtion in region I 
Tj.^.(r) = The cavity temperature distribtuion in region II 
x = The cavity thickness 

k = The thermal conductivity of the cavity (and of the lead wire) 
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Figure 1. Cavity Disc Model for Estimating Effects of Wire 
Conduction. 


The differential equations for and are 

9 2 T 1 3T 

~ 7T + - = -4a T , 

„ 2 r 3r I 

3r 


’Si . 1 3T n . 

+ — 4 “u* 

3r 


where a^. and <x are defined by 


dj = £ Fj/ ( 4kx) , 


a i;[ = eF i;[ /(4kx), 


where the additional parameters are defined as 
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Emissivity of the cavity surface 


Fj = Net power flux in Region I 
= Net power flux in Region II 


The solutions to both differential equations have the same form, i.e. 


T(r) = + A 2 r + A^Jtn r. 


(5) 


Applying the boundary conditions of energy conservation, i.e. 


2„ akx !!l | .ML [T (a) - I ], 
r=a 


( 6 ) 


3T 3T 

IT I - — I * T i (b > ' T u <b >' 
r=b r=b 


(7) 


9T 


II 


3r 


= 0, 


r=c 


( 8 ) 


yields the following expression for temperature distributions 


T^r) = D 1 + D 2 (r/a) + D 3 £n(r/a) 


Tj-j ( r) = (r/a) + £n(r/a) 


(9) 

( 10 ) 


D„ = 


Do = 


Do = 


D, = 


D_ = 


D, = 


4y(C 2 a i;][ - b 2 (a II ~a I )) - j a ] .a 2 (4Y-l) 


“ a I a 

2C 2 ° t II “ 2t|2 ( a ii _a i^ 


b (a IX - a x ) [1-2 £n(b/a) ] 
2 

" a Il a 


2C 2 a 


II 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


y = lx/ (a 2 ) 


(17) 
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The mean "cavity" temperature, which is the quantity that is servo controlled, 
is determined from (9) and (10) as 


T = [2ir JrTj. (r) rdr + 2 tt rdr] / [>(c 2 -a 2 ) ] 


(18) 


Assuming that the servo fixes the mean temperature at T = 0 (the temperature 
origin can be chosen arbitrarily) , equation (18) implies a relationship be- 
tween a^. and <x^, namely 


~~ - 1 = c 4 [4y - | + 2£n f]/[b 2 c 2 [l-4y - 2£n £] + j b 4 ]. 

z a a z 


(19) 


This also implies a relationship between electrical power input and radiative 
power input through equations (3) and (4). Defining 

F g = Electrical power input per unit area of the disk, 

F ^ = Net radiative flux incident on region I, and 

F^ 11 = Net radiative flux incident on region II, 
we find that the net power fluxes F^ and F^ are given by 


F 


I 


F + F 1 
e r 


( 20 ) 


F 


II 


F 

e 



( 21 ) 


Then the L.H.S. of equation (19) can be written as 


- 1 = 


II 


F +F 
e r 

F +F 


II 


-1 = 


F X -F 11 
r r 


F +F 


II 


e r e r 

In order to estimate the power error we shall consider two different 
distributions of radiative power which yield the same total radiative 
power incident on the cavity: 


( 22 ) 
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CASE 1: (Radiation incident in Region I) 


F r I = P r /(irb 2 ) 


f ti = 0 

r 


F = P o J{vc ) 
e e , 1 


CASE 2 : (Radiation incident in Region II) 


F 1 = 0 
r 


IT 2 2 

r = V (w(c " b }) 


F e -P e>2 /(.c) 

Denoting the R.H.S. of (19) by the symbol W, we find that the electrical 

powers for the two cases are 
2 

P e 1 " P r (23) 

6,1 Vrt> r 


= -P 


[1+W] 


e ’ 2 r W(c 2 -b 2 ) 

Since they should both be ideally equal to -P^, the fractional error 
between the two cases is just 


(24) 


_ ?e » 1 ? e,2 (1+W) . r 

P " W l t 2 + 2 2 J * 

r b c -b 

2 

For y >> 1 (typically y is of the order of 10 ) W can be approximated as 


2 2 
W * - h [1 - ^ ( | (1 + \) - £n f)]. 


( 26 ) 


in which case equation (25) can be simplified to the form 
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P -P 4^. + 2^ " £n b 

e,l e,2 c_ 

p k 2 

r 2y(l - 


(27) 


-2 

Recalling that y = £xa and allowing for wire conductivity k to differ 

w 

from the cavity conductivity k we can rewrite (27) in terms of two factors, 
i. e. 


P ,-P 0 k a 

e»l e » 2 _ _w 

P. £kx 

r 


(28) 


where S is a geometrical factor which, for the disc model, is given by 

2 2 

s = [ - i (1 + \) + £n £i a - V" 1 * 

c c 


(29) 


2 2 

If we require b and c to be in the same ratio as the aperture area to 

2 2 

the cavity area, a typical value is c /b % 16. In which case S has 
the value 

S ~ 1.2, c 2 /b 2 JS 16. 


(30) 


For the specific cavity parameters discussed in Section VII, the error 
estimate for a single lead can be calculated. For 


k = k = .426 J/°K 
w 


a = .0038 cm 
l = 1 cm 

x ~ a y = .006 cm 
we find that 





P 


(.0038 cm) 2 
1 cm x . 006 cm 


x 1.2 = 0.29 x 10 


-2 


(31) 


r 

Since angular response deviations of this magnitude typically resulted in 
integrated flux errors somewhat less than 0.1% the lead wire conduction 
does not present a significant problem. 



APPENDIX K. THERMAL MODEL OF THERMOPILE DETECTORS 


As indicated in Figure 1 the basic components of a thermopile radiation 

sensor are (1) a radiation receiver which absorbs and emits radiation, (2) a 

thermal resistor which conducts heat between the receiver and the heat sink, 

(3) a heat sink, and (4) a differential temperature sensor (a thermopile) 

which measures the temperature difference across the thermal ' resistor . 

Analysis of this model makes use of the following notation: 

A = area of the receiver 
r 

P^ = power absorbed by the receiver 
? r = power radiated by the receiver 
P^ = power conducted by the thermal resistor 
T^ = absolute temperature of the receiver ’ 

T g = absolute temperature of the sink 
= receiver emissivity at wavelength A 
K = total thermal conductance of thermal resistor 
V = voltage output of the differential thermopile 
F = incident spectral flux 

ttB^(T) = spectral flux emitted by a blackbody at temperature T 


Under steady-state conditions (power storage taking place only in the 
heat sink) , conservation of energy requires 


P = P - P , 
c a r 


where 


P 

a 


P 

r 


A r/ E »V-’ 

oo 

V I e AW dA * 


( 1 ) 


( 2 ) 
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Figure 1. Thermal model of a thermopile sensor. The 
absorbed power is P , the radiated power P , and the 
conducted power P . 



Thus (1) may be rewritten as follows: 

oo 

P C ■ A r/ E X tF X - ’VV IdX - 

If e is a weak and slowly varying function of A, and if F and B. (T ) have 
A A A r 

approximately the same wavelength dependence, then we may approximate (3) 

by the following expression: 

P = A e (F - aT 4 ) , 
c r r 


where 



OO 

aT r 4 = /*0 (T )dA, and 
“o 


e = 



r 

A more general expression is obtained by allowing £ to vary with wavelength 

A 

and F to have a significantly different spectral dependence than B (T ). 

In this case, if F^ varies in amplitude but retains its basic spectral 

characters (e.g. reflected solar radiation), it is useful to express 

equation (3) in the form 

P = A (e. F - e oT 4 ) 
c r 1 r 

where we define F and e as before and the parameter e^ by the following 
equation : 


e i ■ Av x/ 

0 

An even more general case allows for F to contain both a component dissimilar 
to B^(T ), which we will refer to as the short wave component (e.g. incident 

A IT 

and reflected solar radiation), and a component similar to B (T ), which we 

A r 
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will refer to as the long wave component (e.g. infrared flux emitted by 
the earth) . In this case we find 

P c = A r (e SW F SW + £ LW (F LW " aT r 4)) 

where e is the average emissivity for solar radiation and e is the 

o W J_iW 

average emissivity for blackbody radiation at the receiver temperature. 
The temperature difference required to conduct the power is just P^/K; 


( 10 ) 


thus we have the steady state relation 


K(T r -T s ) - A r (e SB F s|J + - 0T/)) 


(ID 


Since, for most high speed room temperature thermopiles, (T -T^) << T g , 
we may. make the approximation 


aT 4 = aT 4 + 4crT 3 (T -T ) 
r s s r s 


in which case equation (11) may be rewritten in the form 
( T r _T S ) = k l [e SW F SW + e LW (F LW~° T s 4) ] 


where k^ is defined by 


k = A [K + 4A e T „oT 3 } 1 . 

1 r r LW s 

It should be noted that for thermopiles with rapid response the denominator 

of equation (14) will be dominated by the conductivity term, i.e. it will be 

found that 

K » 4A e oT 3 . 
r LW s 

The .voltage output of the sensor is proportional to the temperature 

difference across the thermal resistor, i.e. 

V = a (T - T ), 
r s 

where the proportionality factor a will depend on the thermopile materials, 
the number of junctions, etc., and will have a weak dependence on heat sink 
temperature. A similarly weak temperature dependence is usually present in 


( 12 ) 


(13) 


(14) 


(13) 


( 16 ) 
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K as well. The combined effects of these temperature dependencies is 
typically of the order of 0.3% per °C, and will be entirely negligible 
as a result of temperature control of the heat sink. Substituting (16) 
into (14) yields 


V k 2 ^ E SW F SW + e Lt/ F LW aT s ^ 


where = ak^. 


(17) 



APPENDIX L. UNCERTAINTIES IN SOLAR INPUT RESULTING FROM ATTITUDE ERRORS 


The solar input to the wide FOV sensors is proportional to the cosine 
of the angle y between the vector to the sun s and the sensor pointing 
vector f (See Figure 1) . If the actual angle is y and the attitude 
uncertainty is 0, we shall assume that the estimated pointing vector r can 
lie anywhere on a cone of half angle 0 centered about 2 the actual pointing 
direction. In this case the estimated solar input factor depends on <j> 
according to the equation 


cos0' = cos0cosy - sin0sinycos<f>, 

where 0' is the angle between the estimated vector ? and the solar 
vector 3. 

The square of the error is then 

2 2 
(cos0' - cosy ) = (cosy (l-cos0) + sinysin0cos<J>) . 

The average value over all possible <j>'s is 

2 2 2 1 . 2.2 

<(cos0' - cosy) > = cos y(l-2cos0 + cos 0) + y sin ysin 0. 

<t> * 

Expanding cos0 and sin0 for small values of 0 yields 

cos© = 1 0 2 + 7 , 0 4 - ... 

2 4 ! 

sin0 = 0- -j, 0^ + . . . 

Inserting these expansions into (3) yields 

<(cos0' - cosy) ^>,5: cos^y [-^0 1 + ^ sin y [© - 1/30 ] 


Since 


— CO 


2 n 2 
-a 0 

e d0= 


r( 



n+1 

a 


we find that 


<0 


4 > 

0 


ref) 



(1) 

( 2 ) 

(3) 

( 4 ) 

(5) 

( 6 ) 
(7) 


4 

a it 


( 8 ) 
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where 

2 1/2 

°o ' (<0 V (9) 

assuming a normal distribution of G values about 0=0. Defining 

o 2 

a = (cosG' - cosy) , . (10) 

COS V ' (j) , 0 

we find that 

o ^ = 3/4 a 4 cos 2 y + 1/2 sin 2 y(a 2 - a A 4 ) (11) 

cos 0 b b 

The fractional standard deviation of the solar input factor cos( ') is 

then given by 

(^> ‘ 13/4 a 6 + 1/2(O 0 2 ‘ a 0 4 Han 2 l) 1/2 (12) 

The results of evaluating equation (12) as a function of y are presented 
in the following table for = 0.5° (.008727 radians). 
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Y ANGLE 


a 

cos 


/cosy 


0° 

6.60 

X 

io" 5 

5° 

8.11 

X 

io' 5 

10° 

2.03 

X 

io" 4 

15° 

4.48 

X 

io" 4 

20° 

8.20 

X 

10" 4 

25° 

1.34 

X 

IO -3 

30° 

2.06 

X 

10" 3 

35° 

3.03 

X 

io" 3 

0 

O 

4.34 

X 

io -3 

45° 

6.17 

X 

io” 3 

50° 

8.76 

X 

10" 3 

55° 

1.26 

X 

io" 2 

60° 

1.85 

X 

10’ 2 

65° 

2.84 

X 

io" 2 


% ERROR IN ESTIMATING SOLAR INPUT 

.007 

.008 

.020 

.045 

.082 

.134 

.206 

.303 

.434 

.62 

.88 

1.26 

1.85 


2.84 



APPENDIX M. DESCRIPTION OF FACR TEST MODEL AND RESULTS 


In order to verify the basic theory of the fast active cavity radio- 
meter (FACR) which was presented in Section VI, and to measure the geo- 
metrical coefficients applicable to the angular non-uniformity parameters, 
a test model of a FACR was constructed and tested. In order to avoid 
parts procurement problems and time consuming design and fabrication, 
available materials and simplified design were the rule. Although this 
resulted in a non-optimum design compared to that described in Section VII, 
it did result in a suitable instrument for testing the performance predic- 
tions of the theory. A description of the test model and the test results 
are presented in the following subsections. 

1. Description of FACR Test Model Construction 

The physical configuration of the main components of the FACR test 
model are indicated in Figures 1 and 2. The wire wound active cavity 
is mounted to an adjustable mounting ring by three stainless steel 
support wires. The mounting ring is attached to the heat sink (base) 
by a three point spring suspension and three alignment screws. These 
screws are adjusted to obtain the required spacing and parallelism 
between the cavity aperture and the area defining aperture mounted 
on the cover portion of the heat sink. The area defining aperture is 
subsequently adjusted to be concentric with the cavity aperture. In 
its final configuration, the cavity is surrounded by a heat sink of 
similar shape so that radiative transfer coefficients between the 
active cavity and the sink cavity will not vary significantly over the 
surface of the active cavity. 
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FIGURE 1. Physical Configuration of Assembled Fast Active Cavity Test Model. Cover 
base sink cavity, and interior and exterior surfaces of the active cavity 
with 3M Black (101-C10) . 




REVISIONS ^ i 

LTrI DESCRIPTION I DATE I APPROVED 


M-3 





M-4 


The active cavity is constructed in two pieces: an open-ended half 

which forms the cavity aperture, and closed-ended half which mates to 
the first section at the maximum cavity diameter. Each half is formed 
by winding wire non-inductively on a polyethelene mandril. The wire 
turns are cemented together by a very thin coating of Krylon lacquer. 

The two halves are subsequently coated on their interior surfaces with 
3M-Black (101-C10) and cured. They are then cemented together and the 
exterior of the assembled cavity is painted with 3M-Black. When the 
exterior coating is cured, the cavity becomes quite rigid and durable. 

The two piece heat sink is made of aluminum and vacuum sealed so 
that the volume within the heat sink cavity can be evacuated. The 
heat sink cavity is also painted with 3M-Black (101-C10) to maximize 
the radiative coupling with the active cavity. The straight cylindrical 
section on the base portion of the heat sink is used for winding the 
heater wire used in temperature control of the heat sink. It is also 
used for mounting the high current resistor in the cavity leg of the 
bridge servo used to control the active cavity temperature. The other 
bridge resistors and the preamp electronics are mounted on the back 
face of the heat sink base. 

Although the preferred window is a thin hemispherical dome, the 
actual window available was a thick quartz disc. As a result, angular 
response testing must separately account for the angular response 
effects of the window. The reason a window was used at all was that 
it makes laboratory testing much more convenient. 
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2. FACR Test Model Analysis 

Significant physical and geometrical parameters of the FACR test 
model are presented in Table 1. From these it is possible to derive 
predicted performance characteristics of the test model. 

In order to calculate the electrical power required to maintain 
the cavity at T q when no incident radiation is received, we must first 
estimate the effective emissivity for radiative exchange between the 
cavity and the sink. This can be approximated by equation (3) Section 
VII. For 


1 = / 2*630 \ 2 = 

2 1.905 J 


73 


( 1 ) 


and e, = e_ = e =0.94, we find that 
1 2 _ s ^ 


'eff 




= e 1 + (1-e ) A 


1 /A 2_ 


= 0.90 


( 2 ) 


The radiative power transferred directly from cavity to sink is thus 
given by 

T A Al 

= 1.636 W 


P = e rrOk 
1 eff 


‘ 4 4l 

, T - T 
1 o s 


(3) 


where A^ for the case of a quartz window in place, is taken to be the 

2 

sum of the cavity exterior surface area (29.11 cm ) and the cavity 

2 

aperture area (1.81 cm ). The portion of this power which is actually 

transferred to the window can be estimated by using equation (3) with 
2 

A^ = 1.267 cm . Additional power is conduced via leads and support 
wires. For each term, the power conducted has the form 


P. = A.k. (T -T )/£; 
i li os 


(4) 


, th 


where k^ is the thermal conductivity of the i wire and A^ and are 
the cross sectional area and length respectively. Assuming 


k =3.80 w/(cm - °c) 
cu 


(5) 


k =0.45 w/ (cm - °c) 
ss 


( 6 ) 
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Table 1. Fast Active Cavity Test 

Model Parameter Values 


PARAMETER 

MEASURED OR ESTIMATED VALUE 


Cavity Mass Distribution 


_ 

Copper (#36 T ’ire) 
wire insulation + Krylon 
3M-Black (interior) 
3M-Black (exterior) 

2.20 gm 
0.20 gm 
0.24 gm 
0.15 gm 

- 

TOTAL 

2.79 gm 


Cavity Wire Diameter 



Bare Copper (#36) 
with insulation 

0.0127 cm (.005”) 
0.0152 cm (.006") 


Cavity Resistance @ 20°C 

@ 100°C 

26.55 £7 
34.62 a 


Cavity Major Radius 

1.63 cm (0.640") 


Cavity Aperture Radius 

0.76 cm (0.300") 


Cavity Exterior Surface Area 

29.11 cm 2 


Defining Aperture Radius 

Area' 

Distance Between Cavity 

0.635 cm (0.250") 
1.267 cm 2 


Aperture and Defining Aperture 

0.051 cm (.020") 


Uniform Angular FOV 

+ 68.2° 

- 

Active Cavity Operating Temp.(T Q ) 

100°C 

- 

Sink Operating Temp.(T g ) 

35°C 


Diffuse Hemispherical Emissivity 

of 3M- Black (e ) 
s 

0.94 (Hall, 1970) 


Sink Cavity Major Radius 

1.905 cm (0.750") 


Cavity support wire diameter 

0.025 cm (0.010") 


Support wire length (£ ) 

s 

0.51 cm (0.20") 


Lead wire length (£ g ) 

0.53 cm (0.21") 
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we thus find that each lead wire conducts 59 mW and each support 
wire conducts 28 mW of power. The total electrical power P° is 
thus 1.838 W. These results are summarized in Table 2. 


Table 2. Electrical Power Budget for Zero Radiative 
Input (T = 100°C, T~ = 35“C) 


Radiative Transfer 

cavity to sink (direct) 
cavity to window 
SUBTOTAL 

Conductive Transfer 

lead wires (2) 
support wires (3) 

SUBTOTAL 

Total Electrical Power Required (P°) 

6 


1 . 569 w 
0.067 w 

1.636 w 


0.118 w 
0.084 w 

0.202 w 
1.838 w 


The radiative power received by the cavity from an incident flux 

equivalent to one solar constant H g is given by 

P = e A t H where t is the window transmittance. (7) 

r c a w s w 

2 

We can obtain an approximate value by assuming H = 1357 W/M , t = 0.92, 

S W 

2 

£ c = 0.993 (Section VII), and A = 1.267 cm (Table 1). The result is 

P^ = 157.1 mW/solar constant = 0.1157 mW/ (WM ^ (8) 

The voltage across the cavity is determined from the electrical power 
requirement using the relationship 

V = (R°P°) 1/2 (9) 

c ' 1 e' 

O 

where R^ = 34.62 is given in Table 3. Typical voltage levels are 
summarized in Table 3. 
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Table 3. Typical Values of Test Cavity Power and Voltage 

Space View Solar View 

-2 

(zero incident flux) (1357 wm incident flux) 


Electrical Power to 

Cavity (P ) 1.636 W 1.479 W 

Voltage across 

Cavity (V ) 7.5258 V 7.1556 V 

-2 

AV c per WM change 

in incident flux 0.27 mV 0.28 mV 


According to equation (16) Section VI, the FACR time constant is given 
by 

C 

T K+r ’ 

where C is the total heat capacity of the cavity, K is the effective 
conductance between the cavity and the sink, and F is the servo gain 
parameter. According to equation (8) of Section VI, K can be calculated 
from the power transfers listed in Table 2 as follows: 

4 —1 -1 

K= 1.636W x 4 T (T - T ) + 0.202W(T - T ) = 35.89 mW/o . 

o o s os K 


( 10 ) 


(ID 


The total heat capacity of the cavity is equal to the sum of the specific 
heat capacity - mass products of the components. Although the specific 
heat of copper is well known (0.393 w-sec/gm°c @ 100°C) , the same is not 
true for the other cavity components. We will assume a value of 0.902 
w.sec/gm°C for the other components since this is the value for carbon 
@ 100°C, the major constituent of 3M-Black. The result of this assumption 
is 

C = 2.20 x .393 + .59 x .902 = 1.397 w-sec/°C . (12) 

Since nearly 40% of the total heat capacity is due to the non-metallic 
components, the uncertainties in their specific heats makes (12) only 
an approximate result. 
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The natural time constant of the cavity (no servo action present) 
is given by 

T = = 38.9 sec. (13) 

In order to achieve a t of J 40 ms T must be given by 

r s: 10 3 X K * 36 w/ o„ . (14) 

K. 

3. Servo Electronics 

The circuit diagrams of the DC servo constructed for use with the 

test model are shown in Figures 3 thru 5. The separate portions of the circuit 

given in these figures are the preamplifier, the power amplifier and the 

analog divide and square root circuitry. 

The preamplifier is an instrumentation amplifier which is appropriate 

here because of its good common mode rejection. The MONO OP-07 operational 

amplifiers used in this circuit were chosen for their low frequency noise 

performance. The voltage gain of the preamplifer set up for testing was 

3 

measured to be 2.65 x 10 . 

The power amplifier portion shown in Figure 4 also performs an inversion 
and a level shift. The inversion simply gives the signal the proper sign for 
use with a positive supply voltage. The level shift is performed as a 
convenient technique to keep the voltage for the analog operation chips in 
their operating ranges. The shift which acts as a voltage gain w is accom- 
plished with the transistors shown in the circuit diagram. The magnitude 
of w is given by the resistance ratio R^/l^ where and R^ are defined 

in Figure 4. The value of w is related to z used in the previous 

3 2 

description of the DC servo by w = z . The reason that w does not equal 
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.0047pf 



Figure 3. Preamplifier Circuit Diagram. The ground labelled s is the signal 
ground which is connected to the base of the bridge circuit. The 
voltage difference V^-V^ is the imbalance voltage of the Wheatstone 

bridge containing the cavity. 



Figure 4. Power Amplifier. The hashed line ground symbol is power ground 
and the ground labelled s is signal ground. The voltage V, is 
-the voltage supplied to the bridge. 
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Figure 5. Analog Divide and Square Root Circuit Diagram. The transfer function 
of the Analog Devices AD 533 integrated circuits are given for both 
the divider and the square rooter. The ground labelled s is the 
signal ground. 
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z is that the preamplifier output voltage is not divided by the total 
bridge voltage V., as shown in the block diagram of Section VI, but is 

D 

divided by V„W 

The analog operations for which the circuit is shown in Figure 5 
are performed using two multiplier, divider, squarer, square rooter 
integrated circuits (Analog Devices A D533). The transfer functions 
used for this application are shown for both chips in the figure. They 
were set up using the standard instructions for their use as a divider 
and a square rooter. 

4 . Test Results 

A number of aspects of the predicted FACR performance have been 
verified using the test model FACR. The following characteristics have 
been examined . 

1. Time response 

2. Dependence of the time constant on servo gain 

3. Cavity power consumption. 

The tests of characteristics 1 and 2 were primarily aimed at verifying 
the operation of the FACR servo. The determination of cavity power 
consumption tests the predictions of the radiative and conductive power 
exchange between the cavity and the heat sink. 

The time response of the sensor to chopped incident flux is shown 
in Figure 6 for two different servo gain choices. Semi-log plots of the 
response were made for several different time constants. They showed 
the response to be expotential to within the uncertainty of reading 
the strip chart records except for the initial 10-20 ms of the 



:j HEWLETT-PACKARD Err: 



EIGURE6. Fast ACR Test Model Response to- Chopped -Radiation 
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response. This deviation is caused by the 60 Hz filtering in the 

preamplifier and earlier tests performed without this filtering show 

that it can be eliminated. The time response tests demonstrate that 

the desired time constant of 48 ms is easily within reach of the FACR. 

The results of measuring the FACR time constant for several different 

servo gains are shown in Figure 7. The time constants were measured from 

semilog plots of the strip chart records. The servo gain T was calculated 

from its theoretical expression which requires knowledge of the cavity 

resistance and temperature coefficient of resistance a along with the 

2 

servo parameter gz y. The cavity resistance was measured precisely 

for three different temperatures before testing began. Therefore, a was 

known, and measurements of the cavity resistance performed for each test 

by measuring the bridge current and voltage also provided the cavity 

temperature. The preamplifier gain g was measured during the servo 

2 

construction and z y was determined for each test from measurements of 

the bridge voltage output of the preamplifier. The results of Figure 7 

indicate that, as the theory predicted, the time constant varies linearly 

with T 1 for time constants greater than 50 ms. Again we see that the 

60 Hz filtering in the preamplifier affects the time response and causes 

deviations from the theoretical behavior for short time constants. The 

cavity heat capacity implied by these results is 1.7 W-sec/°K. 

The cavity power consumption was determined from measurements of the 

bridge voltage and current. Table 4 shows the comparison of the test 

results p meas an ^ the theoretical predictions p ca ^ c * The measured cavity 

temperature T and heat sink temperature T given in the table were used 
° s 
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FIGURE 7. Dependence of the Fast ACR on Servo Gain 
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TABLE 4. 



H 

O 

/-S 

0 

S-*' 

T (°K) 
s 

P calc <W) 

P (W) 

me as 

Pressure (microns) 

94.5 

32.4 

1.70 

1.92 

.013 

94.5 

32.4 

1.70 

1.90 

.013 

94.2 

32.3 

1.69 

1.91 

.013 

94.3 

32.3 

1.69 

1.88 

.012 

60.8 

32.3 

.680 

.681 

.011 

60.9 

32.2 

.682 

.681 

.011 

61.1 

32.2 

.687 

.677 

.010 
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to calculate P .. (the measurement of T was discussed in the previous 

C al C O 

paragraph and T was determined using a thermistor recessed in the sink) . 
s 

The measured power consumption was very strongly dependent on the vacuum 
gauge pressure reading (see last column of Table 4), for pressures greater 
than 0.010 microns. Since the gauge was located at the diffusion pump, the 
actual pressure inside the radiometer itself was not accurately known. 
However, when the vacuum gauge readings reached about 0.010 microns or 
smaller, there was good agreement between the measured and the calculated 
power consumption. 

Unfortunately, the initial FACR test phase was not completed, both for 
lack of time and money. The most obvious omission is testing of the angular 
response. Plans were made to verify the angular response characteristics 
of the FACR in a vacuum chamber to obviate the need for determining the 
angular response of the thick window on the test model. Other additional 
testing and design optimization of the FACR is necessary to completely 
demonstrate its capabilities. However, the success of the initial tests 
described in this report make us optimistic that future testing will also 
show agreement with the FACR theory presented here. 



APPENDIX N. AC SENSE - DC HEAT CAVITY SERVO DESIGN FOR THE FACR 


The DC serve designed to control the FACR cavity temperature was 
tested using a model FACR (see Appendix M) and was found to behave as 
predicted by the theory presented in Section VI. As a parallel effort 
to the DC servo testing a study was performed to determine whether other 
types of servo systems would offer advantages over the DC servo. The 
study was restricted to servo systems capable of producing a predictable, 
single time constant response to step function changes in the irradiance 
incident on the cavity. This was the major consideration of the DC servo 
design and as shown in Section VI amounts to requiring that the power 
supplied to the cavity vary linearly with changes in cavity temperature. 

The alternatives considered included a totally AC system, a pulsed 
system and an AC sense - DC heat system. The AC/DC combination was 
found to be superior to the other alternatives and was pursued in some 
detail. 

1. Basic Design and Operation of the AC Sense - DC Heat Servo 

The electronics for the AC sense - DC heat servo is illustrated with two 
levels of block diagrams in Figure 1. Like the DC design a Wheatstone 
bridge with the wire wound cavity as one resistance element (R^) is used 
to sense temperature changes of the cavity. However, here the bridge is 
excited with a fixed AC voltage V^,. The bridge imbalance voltage, which 
is linearly related to the cavity temperature, is then amplified using a 
transformer - preamplifier combination. The amplified voltage is 
synchronously demodulated and then filtered to restrict the broad band 
noise and to attenuate undesirable signal components generated in the 
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rectification process. The voltage emerging from the filter is of 
course still linearly related to the cavity temperature. Since to 
provide a single time constant response, the power supplied to the 
cavity must be linearly related to the cavity temperature, is next 
square rooted. A power amplifier sums the output of the square 
rooter and the AC excitation voltage* and powers the bridge. The DC 
bridge power is isolated to the R^^ b e S °f tbe bridge by placing 
capacitors in the other leg of the bridge. This allows the 
resistances R^ and to be small enough that the bridge does not 
load the transformer appreciably and avoids dissipating large amounts 
of power in R^ and R^. 

Two major advantages of the AC sense - DC heat servo are apparent 
at this point. First, an improvement of the signal to noise ratio over 
that of the DC servo can be obtained by coupling the AC signal from the 
bridge to the servo preamplifier via a transformer. As a result, whereas 
the preamplifier is the dominant source of noise in the DC servo, the 
preamplifier contribution to the noise of the AC/DC servo is very small. 
The system noise is essentially reduced to the Johnson noise of the 
bridge resistors. A noise analysis is given in Section 3 of this 
Appendix. The second advantage is the simplification of the analog 
operations required. By heating the cavity with a DC voltage and 
sensing its temperature with an AC voltage the need for the analog 
divide operation of the DC servo is eliminated. 
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2 . Analysis of the Servo Operation 


In this subsection the previous description of the servo operation 
is translated info servo operating equations similar to those given for 
the DC servo in Section VI. The notation of Figure 1 will be used. 

The bridge imbalance voltage can be written in terms of the AC 
excitation voltage V AC as follows, 

R 2 R 4 

V ( t ) = V - V = V r — - — 1 

ab^ ' a b V AC l I^(t)+K 2 R 3 +R 4 

(R 2 R 3 -R°R 4 ) aR°R 2 T(t) ^ 

‘ V AC [ (R“« 2 )(R 3 + R 4 )- (R . ,2 1 


( 1 ) 

( 2 ) 


where Equation (24) of Section VI has been used and where , as in Section 

VI, R° and a are the cavity resistance and temperature coefficient of 

resistance at the space view operating temperature T q , and T(t) is the 

deviation of the cavity temperature from T^. 

The voltage V m (t) following amplification, rectification and filtering 

is defined in terms of V , as follows, 

ab 

V (t) = n A m V (t) (3) 

m c ab 

where 

n = transformer turns ratio 
A c = preamplifier voltage gain 
m = demodulator gain 

It was assumed here that the transformer input impedance is sufficiently 

large that the transformer is not loaded by the bridge resistances. 

Finally, in terms of the voltage gain of the square rooter z and 

the voltage gain of the power amplifier A , the DC voltage supplied to 

1/2 

the bridge is given by V nr (t) = A z V (t) . Substituting V (t) from 
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Equation (3) and V ^(t) from Equation (2) yields, 


Vp C (t) [A z nm A c V AC ] ^(t)/^ 


(4) 


2 2 


(r 2 r 3 -r°r 4 ) qr°r 2 


V DC (t) [A p z nm A c V AC ] [ (R°+R 2 ) (R 3 +R a ) (R°+R 2 ) 2] 

Equation (5) can be recognized as Equation (25) of Section VI for the 

2 

DC servo where the first bracketed quantity above replaces gz y. The 
DC power to the bridge is again defined by Equation (27) of Section 
VI, yielding 


(5) 


P e C (t> 

= P ° Dc 
e 

- r t ( t) , 

(6) 

where 




p oDC = 
e 

2 2 ( R 3 P r 4 ) 

(z A p nm A c V AC ] (r ^ +r ^ r o ( 1+ g)3 

(7) 


e 3 

a[z 2 A 2 nm A V Ar ,] 
p c AC 

(8) 

r - 

(1+B) 4 

*1 

p 

r*/r 2 . 


(9) 


The major difference between the results for the AC/DC servo and the DC 
servo is the relationship of Equation (4) which gives the DC voltage 
in terms of the bridge imbalance voltage. The results of this analysis 
are summarized in Table 1 which can be compared with the DC servo 
analysis of Table VI. 2. 


3. Noise Analysis 

The important sources of noise in the AC/DC servo system are the bridge 
resistances and the preamplifier. The following calculation show that the 
Johnson noise of the bridge dominates. 
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Table 1. ; Summary of AC Sense - DC Heat Servo Relationships 


( 1 ) 

( 2 ) 


e = 


E l ,R 2 


r = 


( 1 + 6 ) 


a[z 2 A p 2 nm AW k( ,] 
R? 


where. 


V = AC sense voltage applied to the bridge (constant) 

n = Transformer turns ratio 

m = Demodulator Gain 

A p = Preamplifier voltage gain 

A p = Power amplifier voltage gain 

z =:, Gain of analog square rooter 


(3) 

V ab< C > 

(4) 

p*DC 

e 


DC 

(5) 

P T (t) 


DC 

(6) 

P T (t) 

(7) 

v DC (t) 


(l+l/B) 2 R° P° C (t) V AC 
[z 2 A p 2 nm A c V AC ] 

[z 2 A p 2 nm A c V AC J (R^gR^ (R^r^V+P)" V 

Total DC electrical power dissipated in cavity leg 
of the bridge 


(1+1/6) Pg (t) 



(l+l/S)[R“P^(t)] 


1/2 
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The expression for the mean square Johnson noise voltage of the 
bridge has been derived by P. Freymuth (1968). 

V 2 = 4 k T R Af [(1+6) (l+n/6) + AT/T ] (1+B)" 2 
n p- s 

where 

k = Boltzmann's constant 
T g = Sink temperature 

AT = The difference between the cavity temperature and T g 

Af = Equivalent noise bandwidth 

^ 1^2 ^ 3^4 

R = (r x +r 2 ) + (r 3 +r 4 ) 

6 = R 1 /R 2 

n = (R 1 +R 2 )/(R 3 +R 4 ) 

Using the following parameter values from Section VII: 

T = 320° 
s 

AT = 25° 

R x = R 3 = 520 Q 

R„ = R = 460 0 
2 4 

and an equivalent noise bandwidth of 52 Hz (this corresponds to a single 
section filter time constant of 4.8 ms), gives an RMS noise voltage of 
20 nV. 

The amplifier contribution to the noise is found by dividing its 
RMS noise voltage by the turns ratio of the transformer. For a typical 
low noise amplifier with a noise bandwidth of 52 Hz, the RMS noise 
voltage is 315 nV. Using a turns ratio of 50 the preamplifier noise 
referred to the bridge is 6.3 nV. 

Combining the bridge resistance noise with the amplifier noise in 
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RSS fashion yields a total noise voltage of 21 nV. This result was used 
in Section VII to calculate the signal to noise ratio for the FACR 
design. 



